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DEVELOPMENT  AND  APPLICATION  OF  A 
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By 
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December  1982 


Chairman:  Dr.  John  R.  Eyler 
Major  Department:  Chemistry 

A  dual  microcomputer  based  data  acquisition  system  for  an  ion 
cyclotron  resonance  mass  spectrometer  with  a  marginal  oscillator 
detector  is  described.    The  system  uses  an  Apple  II  microcomputer  acting 
as  the  master  processor  and  a  KIM-1  microcomputer  serving  as  the  slave. 
The  two  microcomputers  communicate  through  a  common  S-100  bus  and  its 
control  logic.    A  unique  experiment  interface  panel,  accessible  by 
either  microcomputer,  is  also  described.    The  hardware  and  software  for 
this  system  is  presented  and  discussed  in  detail.    A  single 
microcdmputer  system  based  on  an  Apple  II  microcomputer  controlling  a 
capacitance  bridge  detector  is  also  presented. 

The  dual  microcomputer  system  was  used  to  study  the  ion  chemistry 
of  CH^I,  CF^I,  and  a  mixture  of  the  two.     In  the  pure  CH^I  experiments, 
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reactions  of  the  electron  impact  fragments,  I*  and  CH3,  with  the 
neutral  parent  molecule  occured  at  the  collision  rate.    Experiments  with 
CF^I  showed  reactions  of  the  fragment  ions,  I''"  and  CF^,  occur  at  rates 
slower  than  the  corresponding  reactions  in  the  CH^I  system.    Reactions  of 
CF^I*  with  the  neutral  parent  lead  to  CF^ICF^,  CF3I2,  and  (CF3l)2. 
The  observed  reactions  in  a  mixture  of  CF3I  and  CH3I  indicate  charge 
transfer  reactions  dominate  the  reaction  chemistry.    Ionic  species 
containing  hydrogen  and  fluorine  were  not  seen;  however,  neutral  cross 
products  are  speculated.    These  experiments  demonstrated  the  effective 
use  of  a  marginal  oscillator  detector  at  high  masses  (m/e  >  300). 

The  single  microcomputer  system  was  used  to  study  two  isomers  of 
C^H^.    Rate  constants  for  the  more  reactive  of  the  two  isomers  were 
determined  for  reactions  with  diacetylene,  norbornadiene ,  styrene ,  and 
other  molecules.    Some  product  ions  observed  from  these  reactions  have 
been  reported  in  sooting  flames.    The  less  reactive  of  the  two  isomers 
was  investigated  with  respect  to  its  thermochemistry.    The  proton 
affinity  of  the  C^H^  neutral  forming  this  ion  was  experimentally 
determined  to  be  227.9  kcal/mole.    The  reactive  isomer  is  most  likely  a 
linear  structure,  while  the  unreactive  is  a  cyclopropenyl  cation  with  an 
ethylene  side  chain. 
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CHAPTER  1 
INTRODUCTION 


Ion  cyclotron  resonance  (icr)  mass  spectrometry  began  its  useful 
life  as  a  tool  for  use  by  the  physical  chemist  in  studying  ion-molecule 
reactions.    Due  to  the  limited  resolution  and  mass  range  of  earlier 
instruments  its  use  by  the  analytical  chemist  was  not  widespread. 
However,  following  significant  advances  in  detection  theory,  elec- 
tronics, and  digital  computers,  the  use  of  this  instrument  is  expanding 
rapidly.    The  following  text  describes  the  development  of  the  instrument 
to  the  present,  examines  the  role  of  computers  in  instrumentation,  and 
introduces  the  chapters  to  follow. 

Development  of  ICR  Spectrometry 
Work  in  icr  mass  spectrometry  has  been  extensive  enough  to  warrant 
the  publication  of  a  few  reviews  on  the  subject.    A  comprehensive  survey 
of  the  work  in  the  field  up  to  1973  has  been  presented  by  Hartmann  [1], 
while  Wanczek  recently  reviewed  the  instrumentation  and  theory  from  then 
to  the  present  [2].    In  1971  Beauchamp  presented  a  description  of  the 
applications  of  icr  spectrometry  to  physical  chemistry  [3].  These 
applications  include  determination  of  ion-molecule  collision 
frequencies,  reaction  rate  constants,  reaction  pathways  and  product 
distributions,  energy  dependence  of  rate  constants,  and  thermochemical 
properties  such  as  gas  phase  acidities  and  proton  affinities.    Most  of 
the  chemical  studies  to  date  using  icr  spectrometry  involve  one  or  more 
of  these  areas. 
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Earlier  versions  of  the  icr  spectrometer  involved  a  drift  cell  with 
three  regions  (ion  source,  analyzer,  and  ion  collector)  for  ion 
monitoring.    An  illustration  of  the  cell  and  an  explanation  of  how  it 
worked  has  been  given  by  Beauchamp  [3].    The  cell  required  constant 
production  and  detection  of  the  ions  as  they  drifted  from  the  source  to 
the  ion  collector.    Ions  absorbed  energy  from  a  radio  frequency  (r.f.) 
field  causing  the  resonant  ion  to  collide  with  one  of  the  cell  walls 
where  the  resulting  ion  current  was  measured  as  in  the  original 
omegatron  [4].    The  highest  performance,  however,  is  expected  from  a 
detector  which  does  not  eject  resonant  ions  from  the  ion  cell  [5]. 

The  use  of  a  r.f.  detector  which  monitored  the  power  absorbed  by 
resonant  ions  without  removing  them  from  the  cell  improved  performance 
[6].    The  operation  of  the  icr  was  further  enhanced  by  use  of  modulation 
techniques  which  varied  one  of  the  parameters  affecting  power 
absorption.    This  allowed  phase  sensitive  detectors,  monitoring  the 
modulation  frequency,  to  improve  the  signal  to  noise  ratio.    Types  of 
modulation  included  pulsing  electron  energy  above  and  below  the 
ionization  potential  of  the  ion  [7]  and  pulsing  the  electron  current  by 
controlling  the  grid  of  a  triode  gun,  thus  changing  the  ion  population 
itself  [8].    Other  techniques  involved  the  changing  of  the  cyclotron 
frequency  of  the  ion  by  pulsing  the  trapping  potential  [9]  or  magnetic 
field  [10]. 

Maximum  resolution  for  the  drift  cell  is  still  limited,  however,  by 
the  residence  time  of  the  ion  in  the  analyzer  region.    This  problem  was 
solved  by  the  longer  residence  times  made  possible  by  trapped  ion  cells 
[11].    The  original  application  of  this  cell  generated  ions  by  a  pulsed 
electron  beam,  then  probed  them  some  time  later  by  pulsed  magnetic  field 


modulation  while  detecting  their  power  absorption  with  a  marginal 
oscillator  detector.    In  the  mode  of  operation  commonly  used  with  these 
cells,  the  ions  are  probed  by  a  pulsed  r.f.  detector.    A  more  detailed 
explanation  of  the  behavior  of  ions  in  a  trapped  ion  cell  is  given  in 
Chapter  2 . 

After  detection,  a  quench  pulse,  collapsing  one  of  the  trapping 
potentials,  cleared  the  cell  of  all  ions  before  the  next  pulse  cycle. 
Recent  developments  have  used  the  unquenched  mode  of  operation  where  a 
steady  state  abundance  of  ions,  formed  by  continuous  electron  impact  and 
lost  by  wall  collisions,  makes  possible  the  study  of  slow  equilibrium 
reactions  [12].    The  trapped  ion  cell  and  pulsed  operation  allowed 
better  definition  of  reaction  time.    As  before,  the  ion  signal  was 
modulated  by  pulsing  the  trapping  potentials  [13],  but  eventually  this 
gave  way  to  a  pulsed  marginal  oscillator  (m.o.)  detector  [14]. 

The  marginal  oscillator  is  a  tuned  detector;  thus  only  one  ion  can 
be  seen  at  a  time.    Mass  sweeps  are  usually  accomplished  by  a  slow 
magnetic  field  scan  which  brings  the  ions  into  resonance  at  the 
frequency  of  the  m.o.  detector.    Broad  band  detector  circuits  capable  of 
fast  frequency  sweeps  and  simultaneous  observation  of  many  ions 
eventually  were  constructed. 

One  type  of  broad  band  detection  scheme  which  evolved  from  the 
trapped  ion  cell,  pulsed  timing  techniques,  and  digital  computers  is 
Fourier  transform  icr  (fticr).    The  principle  of  extracting  frequency 
domain  spectra  from  time  domain  signals  has  long  been  understood,  but  it 
wasn't  until  the  development  the  fast  Fourier  transform  algorithm  by 
Cooley  and  Tukey  that  such  routine  analysis  became  practical  [15]. 
Analogies  and  similarities  in  instrumentation  and  theory  between  fticr 


and  Fourier  transform  nuclear  magnetic  resonance  have  been  noted  [16]. 
Data  systems  for  the  later  and  some  specialized  electronics  led  to  the 
development  of  the  first  fticr  [17]. 

Since  the  introduction  of  fticr,  theories  of  its  operation  and 
limitations  have  flourished.    The  effects  of  digitizing  errors  and 
truncation  of  the  free  induction  decay  have  been  examined  [18],  and 
problems  with  peak  height  errors  have  been  corrected  by  zero  filling  of 
the  digitized  transient  [19]  and  by  accounting  for  excitation  level  ■ 
differences  produced  by  frequency  swept  excitation  [20].  Convolution 
based  techniques  for  phase  correcting  the  calculated  frequency  domain 
spectra  of  fticr  signals  have  been  proposed  [21]  and  verified  in 
experiments  [22].    The  line  shapes  at  low  pressures  have  also  been 
expressed  mathematically  [23]. 

Mass  resolution,  which  has  always  been  a  limitaion  of  icr 
spectrometry  in  the  past,  can  reach  95  percent  of  its  optimal  value  in 
fticr  experiments  by  acquiring  data  during  the  first  three  half  lives  of 
a  transient  signal  [24].    The  maximum  resolution  of  fticr  experiments 
can  be  made  arbitrarily  high  by  going  to  arbitrarily  low  pressures 
[21].    The  highest  recorded  resolution  using  a  1.2  T  conventional 
electromagnet  is  760,000  for  benzene  (m/e  78)  at  about  5  X  10"^  torr 
[25].    Instruments  using  a  4.7  T  superconducting  magnet  have  trapped 
ions  for  up  to  13-5  hours  and  obtained  a  resolution  of  1,500,000  for  m/e 
166  at  similar  pressures  [26], 

Wilkins  and  Gross  have  surveyed  the  use  of  Fourier  transform  icr 
spectrometry  as  an  analytical  technique  and  attribute  its  advance  in 
this  field  to  its  high  resolution,  high  sensitivity,  and  sufficiently 
accurate  mass  measurement  [27].    Using  a  three  parameter  empirical 


5 

formula,  mass  accuracies  of  a  few  parts  per  million  have  been  obtained 
over  limited  mass  ranges  [28].    Fticr  has  been  used  as  a  detector  for 
gas  chromatography  [29],  chemical  ionization  [30],  and  to  study 
fragments  from  collision  induced  dissociation  [31  ]• 

Other  broad  band  techniques  surfaced  after  the  invention  of  fticr. 
One  such  method,  rapid  scan  icr,  requires  the  use  of  a  computer  to 
digitize  and  analyse  the  transient  signal  obtained  by  sweeping  across 
the  frequencies  of  ions  of  interest  [32].    The  rapid  scan  and  fticr 
spectrometers  are  close  enough  in  performance  to  allow  the  functions  of 
the  former  to  be  accomplished  by  the  latter  [331.    Another  broad  band 
detector,  a  capacitance  bridge,  provides  the  convenience  of  frequency 
sweep  to  obtain  mass  spectra  but  does  not  require  the  use  of  a  computer 
to  obtain  spectra  [3^,35].    A  smilar  technique  had  been  introduced  in 
1965  by  Wobshall  but  did  not  gain  wide  use  until  after  the  development 
of  fticr  [36] . 

Computers 

The  digital  electronic  computer  has  evolved  over  the  years  from  a 
system  made  up  of  thousands  of  vacuum  tubes ,  to  thousands  of  discrete 
transistors,  to  a  few  silicon  chips  with  as  many  as  100,000  transistors 
on  a  quarter  inch  square.    With  each  new  generation  of  computers  came  a 
decrease  in  size,  an  increase  in  reliability,  an  increase  in  operation 
speed,  and  a  decrease  in  cost  for  the  same  amount  of  computing  power. 
Due  to  the  speed  at  vrfiich  they  could  perform  complex  and  tedious 
calculations,  computers  were  initially  used  for  data  analysis.  Later, 
they  were  used  to  perform  experiments  in  real  time  where  the  events 
monitored  were  too  fast  to  record  by  hand  or  too  agonizingly  slow  to 
justify  the  use  of  one's  valuable  time.    The  computer  is  noted  for  the 


speed  with  which  it  can  perform  operations,  the  precision  at  which  it 
can  reproduce  its  results,  and  its  ability  to  remember  and  faithfully 
execute  any  task  it  has  been  taught  without  complaining  or  becoming 
bored.    Its  only  problem  is  a  tendency  to  be  temperamental  and  do  what 
one  says  rather  than  what  one  means. 

The  early  years  of  computerized  experiments  were  dominated  by 
minicomputers  based  on  discrete  transistors.    The  development  of  medium 
and  large  scale  integrated  circuits  has  resulted  in  the  microcomputer:  a 
machine  with  sufficient  power  for  most  problems  but  at  an  affordable 
cost.    The  use  of  microcomputers  in  laboratory  experiments  is  not  new  at 
the  University  of  Florida.    One  instrument  used  a  microcomputer  for 
control  and  data  acquisition  of  photoacoustic  experiments  [37],  while 
our  own  lab  has  collected  photodissociation  dat^i  from  our  icr  [38]  and 
performed  fticr  experiments  [39,40].    Even  the  fast  kinetic  technique  of 
pulse  radiolysis  has  been  adapted  to  control  by  a  microcomputer  based 
system  [41].    Microcomputers  have  also  been  used  to  control  a  quadrupole 
[42],  magnetic  sector  [43],  and  spark  gap  mass  spectrometers  [44]. 

Preview  of  Chapters  to  Follow 
More  details  about  the  theory  of  icr  ion  motion  and  detection, 
analysis  of  kinetic  data,  and  signal  to  noise  considerations  are  given 
in  Chapter  2.    The  microcomputer  systems  designed  by  the  author  to 
perform  icr  experiments  are  described  in  Chapter  3  and  4.    Chapter  3 
furnishes  information  about  the  hardware  assembled  to  acquire  and 
process  the  icr  data.    Chapter  4  and  the  appendicies  supply  operating 
instructions  and  a  description  of  the  software  used  to  drive  the 
hardware  and  perform  data  analysis . 


Any  new  instrument  used  for  scientific  investigations  must  first 
prove  itself  before  its  results  can  be  regarded  as  a  true  representation 
of  the  event  being  monitored.    Unfortunately,  for  icr  studies  there  is 
no  a  priori  way  to  determine  all  the  factors  which  relate  the  electrical 
response  of  the  detector  and  measured  pressure  to  the  concentration  of 
ions  in  the  cell.    By  performing  an  experiment  whose  reaction  rate  is 
well  known,  however,  the  instrument  can  be  corrected  for  various 
errors.    The  procedure  used  to  verify  the  performance  of  the  icr  data 
system  described  in  this  work  is  given  in  Chapter  5. 

The  remaining  two  chapters  deal  with  applications  of  the  computer 
controlled  icr  to  chemical  problems.    The  C^H^  ion  has  been  shown  to 
have  more  than  one  isomer  in  the  gas  phase  [45]  as  well  as  being  an  ion 
found  in  sooting  flames.    Chapter  6  demonstrates  the  use  of  a  broad  band 
icr  detector  in  experiments  designed  to  identify  one  of  the  C^H^ 
isomers  and  to  measure  the  reaction  rate  of  the  other  with  various 
neutral  species  of  potential  importance  in  flame  systems.    Chapter  7 
details  the  investigation  of  ionic  processes  in  CH^I,  CF^I ,  and  a 
mixture  of  the  two.    These  systems  have  been  studied  before  in  relation 
to  their  radiation  chemistry  and/or  ion  chemistry. 


CHAPTER  2 
THEORY 


Ion  Motion 

An  ion  in  crossed  electric  and  magnetic  fields,  as  in  Figure  2.1, 
is  governed  by  the  Lorentz  equation: 

F  =  q(E  +  vxB) 

where  f  is  the  force  exerted  on  the  ion,  q  is  the  ion's  charge,  v  is  its 
velocity,  x  signifies  the  vector  cross  product,  and  ^  and  S  are  the 
electric  and  magnetic  field  vectors,  respectively  [46].    The  magnetic 
field  causes  the  ion  to  move  in  a  circular  orbit  in  a  plane  perpen- 
dicular to  the  field  with  a  characteristic  angular  frequency,  (Dq, 

0)  (2.1) 

for  an  ion  of  mass  m.    The  electric  field  experienced  by  an  ion  in  a 
trapped  ion  cell  with  flat  rectangular  sides  (see  Figure  2.1)  has  been 
numerically  calculated  to  correspond  closely  to  a  three  dimensional 
quadrupole  field  [46].    Experiments  with  hyperbolically  shaped  walls  in 
an  icr  cell  show  little  improvement  over  flat  plates,  thus  confirming 
the  validity  of  the  quadrupole  approximation  for  electric  fields 
[5,47]. 

The  ion  oscillates  back  and  forth  between  the  two  trapping  plates 
at  a  frequency,  Uj,  given  by 
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Here,  a  is  the  distance  between  the  plate  pairs,  and  V-j.  and  Vq 
are  the  trapping  and  side  plate  potentials,  respectively  [46].  The 
observed  cyclotron  frequency  of  the  ion,  is 

,    2       2  -1/2 

but  since  (jJj<<  cj^  for  practical  experiments,  cj  =       to  a  first 
approximation.    Even  smaller  frequency  shifts,  on  the  order  of  a  part 
per  million,  due  to  Coriolis  coupling  of  ion  cyclotron  motion  to 
rotational  degrees  of  freedom,  have  been  postulated  [48] . 

If  the  ion  is  probed  by  a  radio  frequency  electric  field  in  a  plane 
perpendicular  to  the  magnetic  field,  it  can  absorb  energy  from  the 
field  for  frequencies  near        The  ions  are  accelerated  to  an  average 
kinetic  energy,  T,  of 

^         Sm  2  0 

where  t  is  the  duration  of  the  r.f.  pulse  of  amplitude  E,  and  Vg  is  the 

initial  velocity  of  the  ion  [13].    The  first  term  is  the  energy  absorbed 

by  the  ion  from  the  r.f.  field,  while  the  second  term  corresponds  to  the 

initial  kinetic  energy  of  the  ion.    Due  to  differences  in  the  phase  of 

the  oscillation  of  various  ions  of  the  same  charge  to  mass  ratio  a 

particular  ion's  kinetic  energy  may  vary  up  to    ±qEvQt/4  [49]. 

As  the  ion  accelerates,  the  size  of  its  cyclotron  orbit  increases. 

The  radius,  r,  of  the  orbit  can  be  calculated  from 

nl/2 


r  = 


2T 

2 


The  maximum  orbital  radius  attainable  is  limited  by  the  physical 
dimensions  of  the  cell  and  is  equal  to  one  half  the  separation  between 
the  upper  and  lower  plates. 
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Energy  can  also  be  imparted  to  the  ions  from  a  r.f.  electric  field 
applied  to  the  trapping  plates  at  a  frequency  given  by  Equation  2.2. 
This  principle  has  been  applied  to  obtain  mass  spectra  from  an  icr  and 
is  used  in  the  ejection  of  electrons  from  the  cell  when  negative  ions 
are  being  trapped  [9]. 

Ion  Detection 

The  motion  of  an  ion  has  been  modelled  as  a  rotating  electric 
monopole  between  the  plates  of  the  icr  cell  [50].    However,  in  icr 
experiments  one  does  not  deal  with  just  one  ion  but  rather  a  macroscopic 
ensemble  of  ions  in  vrtiich  case  the  macroscopic  rotating  electric 
polarization  must  be  considered.    As  the  ions  are  excited  their  initial 
random  phases  are  driven  into  coherence.    This  coherent  packet  of  ions 
then  induces  an  oscillating  image  current  across  the  detector. 

A  generalized  circuit  diagram  for  icr  detectors  is  shown  in  Figure 
2.2A  and  applies  for  both  tuned  and  broad  band  detectors.    In  a  tuned 
circuit,  the  m.o.  detector  for  example,  a  parallel  tank  circuit  made  up 
of  a  capacitor,  inductor,  and  resistor  are  driven  by  a  constant  current 
source.    This  produces  an  instantaneous  voltage  change,  AV(t),  across 
the  tank  due  to  power  absorption  by  the  ions  which  is  given  by  the 
expression 

2 

ATT/.^        Nq  EtR  . 
AV(t)  =  —  sin  cot 

2md 

where  R  is  the  value  of  the  resistor,  N  is  the  total  number  of  ions  of 
mass  m,  and  all  other  variables  are  as  previously  defined  [51].  From 
the  above  expression  one  can  see  that  the  signal  voltage  is  inversely 
proportional  to  the  ion  mass.    Amano  has  shown  by  a  different  model  that 
the  signal  induced  in  the  detector  by  the  motion  of  the  ion  is  exactly 
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A) 


B) 


gure  2.2.     A)  General  circuit  diagram  for  tuned  and  broad  band  icr 

detectors.     A  radio  frequency  oscillator  drives  a  capaci- 
tor, inductor,  and  resistor.     A  change  in  the  number  of 
resonant  ions  in  the  icr  cell  causes  a  change  in  the 
amplified  output  signal.     B)  Circuit  diagram  for  the  paral- 
lel resistor  model  for  icr  power  absorption  detectors. 
Resonant  ions  cause  a  drop  in  the  amplified  output  signal 
by  shunting  power  to  ground. 
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out  of  phase  with  the  driving  voltage.    Thus  the  oscillation  voltage 
across  the  tank  decreases  [52]. 

Another  theory  applied  to  conventional  icr's  which  use  power 
absorption  detectors  is  illustrated  in  Figure  2.2B  [50].    The  power 
absorbed  by  the  ions  from  the  r.f .  generator  is  modelled  by  a  resistor 
which  shunts  the  icr  cell.    When  power  absorption  is  not  detected  the 
shunt  resistance  is  infinity.    The  m.o.  detector  has  been  modelled  this 
way  and  yields  the  same  detector  response  as  predicted  by  the  other 
theories . 

Comisarow's  rotating  electric  monopole  theory  yields  two  different 
cases  for  detectors  of  the  type  shown  in  Figure  2.2A.    If  the  circuit  is 
predominately  resistive  (or  tuned)  then  the  signal  voltage  behaves  as  in 
the  m.o.  detector:  the  voltage  is  inversely  proportional  to  ion  mass  and 
in  phase  with  the  current;  thus  the  ions  absorb  power.    On  the  other 
hand  if  the  detector  is  predominately  capacitive,  as  in  broad  band 
detectors,  then  the  signal  voltage  is  independent  of  ion  mass  and  the 
current  leads  the  voltage  by  90  degrees.    Thus  no  power  is  lost  in  the 
capacitor,  but  only  through  the  resistive  components  of  the  circuit. 
Mclver  et  al.  have  also  shown  that  the  equations  of  motion  for  the  ions 
are  coupled  to  the  circuit  equations  of  the  detector  and  therefore 
dependent  upon  the  circuit  used  [53]. 

The  resolving  power,  Rp,  of  a  mass  spectrometer  is  defined  as  the 
mass  of  the  ion  divided  by  the  width  of  the  ion  signal,    m,  in  terms  of 
mass  (commonly  full  width  at  half  peak  height)  as  shown  in  the  following 
equation . 
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For  icr  spectrometers  this  quantity  is  numerically  equal  to  the 
resolution  defined  in  terms  of  frequency  [23].    Thus  for  a  particular 
ion  mass  the  ultimate  resolution  is  determined  by  the  width  of  the  naass 
peak. 

At  low  pressures,  the  line  shapes  for  trapped  ion  and  drift  cells 
using  a  m.o.  detector  are  the  same  and  are  given  by  the  expression 

where  Acj^qcj  is  the  line  width  in  terms  of  frequency  and  T  is  the 
data  acquisition  time  in  seconds  [23].  When  using  MKS  units  the 
resolution  is  equal  to 


m 


qBT 


Am  5.566m 

Fourier  transform  icr's  have  a  somewhat  narrower  line  width  which  is 
given  by  the  expression 

A         -  3.971 

Even  though  equations  2.4  and  2.5  differ  by  only  a  factor  of  1.5,  the 
data  acquisition  period  for  the  m.o  detector  is  on  the  order  of  a  few 
milliseconds,  while  the  fticr  detector  can  acquire  data  for  hundreds  of 
milliseconds,  thus  greatly  increasing  the  resolution  attainable  by  fticr 

High  pressures  in  the  icr  cell  cause  the  mass  peaks  to  broaden  due 
to  collisions.    The  effect  of  this  pressure  broadening  takes  the  same 
form  in  drift  and  trapped  ion  cells  as  it  does  in  the  fticr  experiment 
and  has  been  expressed  analyticly  [24].     When  the  collisions  are 
nonreactive  then  the  peaks  have  a  Lorentzian  distribution  with  the  half 
line  width  at  half  peak  height  equal  to  the  elastic  collision  frequency 
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[49].    Even  when  reactive  collisions  are  taking  place  the  ion  signals 
are  still  Lorentzian  but  the  expression  for  the  peak  width  is  more 
complicated. 

Signal  and  Noise 

While  attempting  to  monitor  the  signal  from  the  detector  one  must 
always  contend  with  the  noise  that  comes  with  it.    Malmstadt  et  al.  have 
discussed  the  sources  of  noise  found  in  any  detection  circuit  [54]. 
Some  sources  of  noise  are  inherent  characteristics  of  the  electronic 
components  used  vrtiile  others  are  picked  up  from  the  environment,  and  are 
therefore  avoidable.    Johnson,  or  thermal,  noise  is  produced  by  random 
motion  of  electrons  in  resistive  elements  of  the  circuit.    It  can  be 
reduced  by  using  lower  value  resistors  wherever  possible  or  by  lowering 
the  operating  temperature  of  the  circuit.    Thermal  noise  is  totally 
absent  at  zero  Kelvin,  but  this  is  hardly  a  practical  operating 
condition.    Random  motion  of  electrons  across  a  junction  is  called  shot 
noise,  produced  at  junctions  of  gates  in  transistors  and  such  is  a 
problem  one  must  live  with;  at  electrical  contacts  like  solder  joints, 
it  is  a  problem  vdiich  can  reduced  by  good  soldering  technique. 

Excess  noise  is  considered  to  be  anything  above  and  beyond  that  of 
Johnson  and  shot.    There  are  two  classes  of  this,  the  first  being 
interference  noise.    Included  here  are  60  Hz  and  harmonics  from  power 
supplies  and  power  cords,  and  radio  frequency  interference  from  electric 
motors,  switches,  and  local  broadcast  stations.     The  second  is  flicker 
noise,  usually  referred  to  as  drift.    It  can  be  caused  by  power  supply 
voltage  fluctuations,  changes  in  component  values,  or  temperature 
drifts.    Flicker  noise  is  partially  due  to  natural  aging  of  components 
while  other  sources  can  be  reduced  by  placing  the  electronics  in  a 
stable  environment. 
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One  final  source  of  noise,  quantizing  noise,  is  introduced  when  one 
attempts  to  quantify  the  analog  output  of  the  detector.    The  magnitude 
of  this  noise  is  dependent  upon  the  instruments  used  and  eyesight  of  the 
investigator  when  quantification  is  done  by  hand  from  a  plot  of  the 
data.    When  the  data  are  digitized  by  an  analog  to  digital  converter  the 
finite  resolution  of  the  A/D  defines  the  amount  of  digitizing  noise 
which  can  be  quantified  as    ±1/2  bit.    Often  the  analog  signal  from  the 
detector  is  such  that  quantizing  noise  is  insignificant  relative  to 
other  noises . 

If  the  sources  of  noise  are  random  and  independent  of  each  other 
the  noise  adds  quadratically .    If  not,  then  cross  terms  must  be  included 
when  calculating  the  total  noise  present. 

The  signal  to  noise  ratio  of  icr  signals  has  been  given  much 
attention  since  the  development  of  fticr.    Marshall  has  shown  that  S/N 
in  fticr  is  improved  for  short  data  acquisition  times;  however, 
resolution  is  enhanced  by  long  times  [16].    White  et  al.  have  shown  the 
simultaneous  increase  in  S/N  and  resolution  resulting  from  a  reduction 
in  collisions  as  the  pressure  is  decreased  [25]. 

In  icr  signal  modeling  the  only  source  of  noise  vdiich  has  been 
included  in  S/N  calculations  is  thermal  noise.    Comisarow  has  derived  an 
expression  for  the  S/N  for  RC  parallel  circuit  icr  detectors  as  being 

S    ^  Nq^rBR-*-^^ 

^         2ma  /2k  TAf 
B 

vrtiere  T  is  the  temperature  of  resistor  R,  Af  is  the  detection  bandwidth, 
kg  is  the  Boltzmann  constant,  and  all  other  variables  are  as 
previously  defined  [50].    If  two  ions  of  different  masses  are  detected 
under  the  same  conditions  except  resonance  frequency  (but  the  same 
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magnetic  field  strength),  then  the  relative  signal  to  noise  ratio  is 

(S/N),  m 

^  ^  (2.6). 


(S/N)^ 

Thus,  if  one  doubles  the  mass  of  the  ion  being  investigated,  the  S/N 
drops  by  a  factor  of  two  relative  to  the  lighter  ion.    This  illustrates 
the  difficulty  of  monitoring  heavy  ions  in  an  icr. 

One  need  not  abandon  all  hope  for  studying  heavy  ions  if  one  has 
access  to  a  computer.    In  this  case  signal  averaging  may  be  performed  to 
increase  the  signal  to  noise  ratio.    However,  there  are  certain 
requirements  which  must  be  satisfied  for  accurate  averaging  [55] .  These 
requirements  are: 

1 .  It  must  be  possible  to  take  repetitive  scans  such  that 
exactly  the  same  regions  of  the  transient  signal  are  being 
scanned. 

2.  Data  must  be  acquired  at  a  fixed  sample  time  interval  and 
initiated  by  a  trigger  pulse  to  insure  that  the  information 
goes  into  the  same  channel  as  in  previous  scans. 

3.  The  signal  must  have  only  random  noise  centered  around  zero 
volts. 

4.  The  analog  to  digital  converter  must  represent  zero  volts 
as  zero  counts  or  be  corrected  in  software. 

If  these  conditions  are  met  then  the  signal,  S,  grows  linearly  with  the 
number  of  scans,  Nq. 

s  =  ^S^O 

If  the  noise,  N,  is  random,  with  no  coherent  contributions,  then  its 
growth  is  proportional  to  the  square  root  of  the  number  of  scans. 


N  =  k^N 
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This  leads  to  a  signal  to  noise  ratio  given  by 

S    _  ^1/2 


N  0 

where  the  above  k's  are  constants  of  proportionality. 

The  conditions  mentioned  above  can  easily  be  accomplished  by  a 
computer  controlled  icr.    Thus  one  can  obtain  a  signal  to  noise 
enhancement  which  is  propotional  to  the  square  root  of  the  number  of 
averaged  passes. 

Ton  Kinetics 

Experimental  conditions  in  an  icr  spectrometer  are  such  that  the 
concentration  of  the  ions  is  orders  of  magnitude  less  than  the 
concentration  of  the  neutral  gas.    Thus  when  the  ions  are  reacting  with 
the  neutral  gas  the  observed  ion  decay  signal  follows  pseudo-first  order 
kinetics.    In  this  case  the  intensity,  I,  of  the  ion  at  any  time,  t,  is 
given  by  the  expression 

I  =  I^e 

where  Iq  is  the  initial  intensity  at  time  zero  and  k  is  the  ion  decay 
rate.    Since  the  intensity  follows  an  exponential  decay,  one  way  the 
rate  can  be  determined  is  to  plot  the  natural  logarithm  of  the  ion 
intensity  versus  time.    The  plot  will  give  a  straight  line  with  a  slope 
of  -k.    Another  means  of  determining  the  decay  rate  is  to  calculate  the 
natural  logarithm  of  I/Ig  for  various  point  pairs  where  Iq  is  now  any 
concentration  taken  at  a  time,  At,  before  I.    The  decay  rate  for  this 
method  is  given  by 
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The  decay  rates  from  above  are  converted  to  a  rate  constant  by 
dividing  by  the  neutral  gas  density.    The  very  low  gas  pressures  allow 
use  of  the  ideal  gas  law  for  this  as  follows: 

PN, 

n  A 

V  RT 

where  n/V  is  the  gas  density,  P  is  its  pressure,  R  the  gas  constant,  T 
the  absolute  temperature,  and       is  Avagadro's  number.    If  the  pressure 
is  in  torr,  T  is  in  Kelvin,  and  R  is  in  cm^- torr/K- mole ,  then  the 
resulting  density  is  in  molecules  per  cm^.    Dividing  the  decay  rate  by 
the  gas  density  gives  the  reaction  rate  constant  in  units  of 
cm^/molecule-sec . 

The  two  methods  outlined  above  for  determining  the  reaction  rate 
constant  give  slightly  different  values  in  practical  use.    The  slope 
from  the  semilog  plot  is  usually  calculated  by  a  linear  regression 
routine  for  the  best  straight  line  fit  through  the  noisy  data.  The 
ratio  method  is  usually  repeated  for  several  point  pairs  separated  by  a 
constant  interval  of  time,  with  the  results  then  averaged.    The  first 
method  is  not  sensitive  to  the  order  of  the  noise  in  the  data  while  the 
second  one  is.    This  accounts  for  the  slightly  different  values  obtained 
when  using  the  two  methods  on  the  same  data.    The  use  of  linear 
regression  calculations  is  preferred .and  gives  better  results. 


CHAPTER  3 
SYSTEM  HARDWARE 


IntroduGtioq 

Two  microcomputer  controlled  data  acquisition  systems  are  described 
in  this  chapter.    The  first  to  be  presented  is  a  single  microcomputer 
setup  designed  for  the  National  Bureau  of  Standards  (NBS)  and  used  for 
the  C^H^  studies  described  in  Chapter  6.    Then  a  dual  microcomputer 
system,  used  for  the  CH^I  +  CF^I  experiments  described  in  Chapter  7,  and 
a  special  experiment  interface  are  discussed  in  detail.    These  two 
microcomputer  systems  were  designed  with  the  idea  of  future  expansion  in 
mind;  therefore,  comments  concerning  expansion  capabilities  will  be  made 
when  appropriate. 

Single  Microcomputer  System 
A  microcomputer  based  data  acquisition  and  analysis  system  for  an 
existing  icr  mass  spectrometer  [56]  with  a  new  capacitance  bridge 
detector  (CBD)  [3^]  was  installed  at  the  National  Bureau  of  Standards  by 
the  author.    The  microcomputer  system  is  based  upon  a  single  Apple  II 
(Apple)  [57],  which  utilizes  mainly  commercial  peripherals  and 
controllers,  and  carries  out  basic  experimental  control  and  data 
acquisition  functions  similar  to  the  dual  microcomputer  system  described 
later  in  this  chapter.  Due  to  the  unique  nature  of  the  instrument, 
however,  a  few  special  interfaces  were  built  to  couple  the  computer  with 
the  mass  spectrometer.    The  Apple-based  system  was  chosen  due  to  past 

proven  performance,  completeness  of  documentation,  and  low  cost. 

20 
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The  microcomputer  system,  shown  in  Figure  3.1,  is  intended 
primarily  for  use  as  a  pulsed  icr  mass  spectrometer  controller.  A 
12  inch  color  monitor  provides  video  display  and  graphics  while 
expansion  of  the  basic  computer  performance  is  furnished  by  the 
peripheral  connector  slots  on  the  microcomputer's  motherboard.  The 
software  is  designed  with  the  understanding  that  several  boards  must  be 
in  certain  slots.    A  list  of  the  slots  and  cards  expected  by  the 
software  along  with  the  function  performed  is  given  in  the  following 
paragraphs . 

Slot  1.    Mountain  Computer's  CPS  Multifunction  Card  [58]  provides  a 
parallel  port  for  driving  a  dot  matrix  printer.    The  printer  supplies 
hardcopy  of  experiment  parameters  and  spectra  while  a  24  hour 
clock/calender,  with  battery  backup,  aids  in  data  logging.    A  serial 
port  is  still  available  for  future  endeavors. 

Slot  3.    Interactive  Microware's  ADALAB  [59]  provides  an  A/D  and 
parallel  port  for  use  in  the  experiment.    The  12-bit  plus  sign  A/D 
digitizes  the  CBD  output  while  one  8-bit  parallel  port  delivers  the 
detector's  pulse  sequence.    An  unused  12-bit  plus  sign  D/A  and  one  8-bit 
parallel  port  are  available  for  expansion. 

Slot  4.    Solid  State  Music's  IEEE-488  [60]  controller  communicates 
with  devices  which  interface  through  an  IEEE-488  bus.    Presently  it 
controls  a  frequency  synthesizer  [61]  which  supplies  the  CBD's 
resonance  frequency.    The  board  is  capable  of  controlling  up  to  16 
different  devices  through  one  cable. 

Slot  6.    Apple  II  disk  drive  controller  with  two  floppy  disk  drives 
for  use  in  system  initialization  and  mass  storage  of  programs  and  data 
files. 
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Figure  3.1.     Block  diagram  of  the  single  microcomputer  controlled 
icr  spectrometer  at  the  National  Bureau  of  Standards. 
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Slot  7.    Contains  the  custom  made  spectrometer  interface  for 
driving  the  parallel  port  to  the  CBD  and  for  accepting  binary  coded 
decimal  input  from  an  ionization  gauge  controller.    This  card  is  the 
only  component  of  the  data  system  vrtiich  is  not  commercially  available. 

Not  required  by  the  software  but  available  for  future  development 
is  a  16  K  RAM  card  in  slot  0  [62].    This  expands  the  Apple's  memory  to 
64  K  of  RAM  and  allows  use  of  languages  other  than  floating  point  BASIC 
such  as  INTEGER  BASIC  and  PASCAL.    Two  remaining  slots,  2  and  5,  are 
unoccupied  and  can  be  used  for  additional  devices. 

The  program  for  running  the  system  is  similar  to  that  described  in 
the  next  chapter.    Modifications  to  the  program  take  advantage  of  the 
synthesizer's  fl^equency  sweep  capability  for  obtaining  mass  spectra  and 
determining  product  ions  through  double  resonance  experiments.  Other 
software  changes  manipulate  the  distinct  hardware,  correct  for  clock 
frequency  differences,  and  compensate  for  single  processor  control.  The 
use  of  a  single  microcomputer  requires  that  it  performs  all  tasks  needed 
to  run  an  experiment.    Since  it  can  only  perform  one  task  at  a  time, 
e.g.  pulse  the  detector  or  analyze  the  data  ,  all  tasks  must  be  stopped 
periodically.    The  output  from  pulsed  detectors,  however,  have  a 
tendency  to  drift  when  they  are  not  being  pulsed  and  must  be  corrected 
in  software.    Surely,  a  dual  or  multi-processor  system  could  overcome 
this  problem. 


Dual  MiGroGomputer  System 
The  microcomputer  system  designed  for  our  laboratory  evolved 
through  time  from  a  single  KIM-1  microcomputer  [63]  and  KIMSI  S-100  bus 
[64]  to  a  dual  microcomputer  system  with  numerous  peripherals.  The 
complexity  of  the  system  was  governed  by  the  most  difficult  experiment 
it  was  to  perform:  Fourier  transform  mass  spectrometry.    Our  fticr 
experiments  and  hardware  have  been  reported  earlier  [39,40].    Later,  the 
fticr  computer  system  was  renovated  to  a  general  data  acquisition  and 
control  instrument  capable  of  several  different  experiments  by  modifying 
the  software  and  changing  a  few  connectors. 

A  block  diagram  of  the  system  in  its  latest  configuration  for 
controlling  a  marginal  oscillator  detector  is  shown  in  Figure  3.2.  The 
icr  spectrometer  is  the  same  used  by  an  earlier  colleague  but  modified 
for  computer  control  [38].    Two  microcomputers  are  arranged  in  a 
master-slave  configuration  with  an  Apple  acting  as  the  master  and  a  KIM 
serving  as  the  slave.    The  master  microcomputer  is  responsible  for 
interaction  with  the  user,  program  development  (for  both  itself  and  the 
KIM),  graphics,  storage  of  programs  and  data  on  disk,  and  data 
analysis.    The  slave  microcomputer  performs  the  routine  tasks  needed  to 
run  the  icr  experiment.    These  include  system  initialization,  pulse 
generation,  signal  averaging,  and  system  status  display.  Both 
microcomputer  systems  have  6502  microprocessors  as  their  central 
processing  units  (CPU's),  thus  allowing  the  machine  code  from  one  to  run 
on  the  other. 

Communication  between  the  master  and  the  slave  is  accomplished 
through  a  common  S-100  bus  and  its  control  logic.    The  Apple 
microcomputer,  being  the  source  of  commands,  determines  the  use  of  the 
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S-100  bus  and  halts  instruction  execution  of  the  KIM  when  the  bus  is 
desired.    All  data  and  programs  which  pass  between  the  two 
microcomputers  are  placed  in  the  common  8  K  memory  block  on  the  S-100 
bus.     A  waveform  recorder  [65]  for  fast  analog  to  digital  conversion,  an 
experiment  interface  with  analog  and  digital  I/O,  and  a  frequency 
synthesizer  [66]  for  the  system  clock  and  sine  wave  generation  are 
accessible  through  the  S-100  bus.    All  components  on  the  S-100  bus, 
except  for  non-common  memory,  and  the  experiment  interface  can  be  used 
by  either  processor. 

Apple  Microcmputer  System 

The  Apple  microcomputer  used  for  this  system  is  one  of  the  earlier 
versions  which  has  a  set  of  memory  select  headers  on  the  motherboard, 
otherwise  it  is  the  same  as  the  Apple  described  for  the  single 
microcomputer  system.    A  9  inch  black  and  white  TV  monitor  serves  for 
video  display  and  graphics.    In  the  peripheral  slots  are  the  16  K  RAM 
memory  card  (slot  0),  the  multifuntion  card  (slot  1),  and  the  disk  drive 
controller  with  one  floppy  disk  drive  (slot  6)jall  functioning  as 
previously  described.    A  fast  floating  point  processor  board  in  slot  2 
contains  the  hardware  necessary  for  speeding  up  the  calculation  of  the 
transcendental  functions  called  from  the  languages  APPLEFAST  3.3  or 
PASCAL  [67].    The  appropriate  software  must  be  loaded  into  the  16  K  RAM 
card  before  the  fast  functions  can  be  used.    Interfacing  with  the  S-100 
bus  and  icr  spectrometer  is  accomplished  through  custom  logic  built  on  a 
prototype  card  in  slot  5.    Slots  3,  4,  and  7  are  unoccupied  and  ready 
for  expansion. 

To  the  microcomputers  the  world  appears  as  a  series  of  unique 
addressesj  each  of  which  allows  access  to  a  particular  device  or  memory 
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location.    Addresses  of  interest  to  a  programmer,  which  deal  with  the 
custom  logic  of  this  system,  are  shown  in  Tables  3.1  and  3-2.  All 
addresses  and  numbers  used  in  this  section  are  in  hexadecimal.  The 
Apple  addresses  listed  are  determined  by  the  peripheral  slot  occupied  by 
the  custom  interface.    The  DEV  address  in  Table  3.1  for  slot  5  is  CODX, 
where  the  X  corresponds  to  a  number  from  0  to  F.    The  S-100  addresses  in 
Table  3-2  all  begin  with  an  I/O  page  number.    For  the  Apple  using  slot  5 
the  first  byte  of  the  address,  often  refered  to  as  the  page  number,  is 
C5  while  the  KIM  sees  the  devices  in  page  FO.    This  difference  is  due  to 
the  way  the  I/O  page  is  decoded  by  the  Apple  and  KIMSI  and  must  be  taken 
into  account  when  using  code  written  for  one  microcomputer  on  the  other. 

Note  that  a  read  or  write  instruction  to  a  particular  address  may 
control  two  different  devices.    Write  devices  are  those  which  accept 
information  from  the  CPU,  such  as  D/A's  and  parallel  output  ports,  and  are 
activated  by  a  write  instruction.    A  read  device  is  one  which  sends 
information  to  the  CPU,  an  A/D  or  parallel  input  port  for  examples.  The 
addresses  marked  as  "unused"  are  decoded  by  the  hardware  but  not  assigned 
to  any  particular  device.    These  addresses  are  ready  to  aid  in  any  future 
expansion.    Addresses  marked  as  "not  available"  have  not  been  fully 
decoded  and  will  require  extra  logic  for  expansion. 

Access  to  these  addressed  devices  can  be  achieved  through  BASIC, 
machine  language,  or  the  Apple  monitor.    To  send  a  reset  command  to  the 
KIM'S  CPU  one  would  execute  the  read  instruction  "PEEK  49362"  from  BASIC 
or  "LDA  $C0D2"  from  assembly  language.    Typing  "CODO"  followed  by  a 
carriage  return  into  the  Apple  monitor  will  also  initiate  the  above 
function.    Writing  to  the  same  addrress,  though,  will  execute  a  different 
command,  namely  request  the  S-100  memory.     The  Apple  may  send  a  trigger 
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Table  Interface  device  address  map  for  the  Apple  .  All 

addresses  are  in  hexadecimal.     "Unused"  signifies 
the  address  is  decoded  but  not  in  use.  "Not 
available"  designates  an  address  which  has  not  been 
totally  decoded. 


Apple 

Address  Rgad  Pevicg  Write  Device 

(DEV)O  Interface  status  S-100  request 

1  KIM  status  S-100  release 

2  KIM  reset  S-100  memory  request 

3  KIM  NMI  S-100  memory  disable 

4  Unused  Enable  interface  reset 

5  Unused  Disable  interface  reset 

6  Unused  Apple  memory  request 

7  Unused  Apple  memory  disable 

8  Unused  Not  available 

9  Unused  Not  available 
A  Unused  Not  available 
B  Unused  Not  available 

C  Unused  Not  available 

D  Unused  Not  available 

E  Unused  Not  available 

F  Unused  Not  available 
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Table  3.2.  Interface  device  address  map  for  the  S-100  bus.  All 
addresses  are  in  hexadecimal.    "Unused"  signifies 
the  address  is  decoded  but  not  in  use.  "Not 
available"  designates  an  address  which  has  not  been 
totally  decoded.         ■  •  ^       ■  ,. 


S-100 

'    ;  \  :  -> 

Address 

Read  Device 

Write  Device 

(1/0)00 

Serial  Port  A  Control 

Serial  Port  A  Status 

01 

Serial  Port  A  Data  In 

Serial  Port  A  Data  Out 

02 

Serial  Port  B  Control 

Serial  Port  B  Status 

03 

Serial  Port  B  Data  In 

Serial  Port  B  Data  Out 

DxUtUduxUIi  oonui^ox 

05 

Biomation  Status 

Parallel  Port  B 

10 

Open  relay 

X-D/A  output 

11 

Close  relay 

Y-D/A  output 

12 

nXjgii  v^xu^K.  uyuc 

13 

A/D  input 

Low  clock  byte 

1M 

Unused 

Trigger  A 

15 

Unused 

Trigger  B 

16 

Unused 

Ilnn  ^pd 

17 

Unused 

Unused 

18 

Not  available 

Unused 

19 

Not  available 

Unused 

1A 

No  t  ava  i  1  a  h  1  p 

IB 

Not  available 

Unused 

1C 

Not  available 

Unused 

ID 

Not  available 

Unused 

IE 

liKj  \^    d  V  dJ.  Xd  u  X  c 

unusea 

IF 

Not  available 

KIM  status  port 

20 

PIA-1A  data/direction 

PIA-1A  data/direction 

21 

PIA-1A  control 

PIA-1A  control 

22 

PIA-1B  data/direction 

PIA-1B  data/direction- 

23 

PIA-1B  control 

PIA-1B  control 

24 

PIA-2A  data/direction 

PIA-2A  data/direction 

25 

PIA-2A  control 

PIA-2A  control 

26 

PIA-2B  data/direction 

PIA-2B  data/direction 

27 

PIA-2B  control 

PIA-2B  control 

28 

PIA-3A  data/direction 

PIA-3A  data/direction 

29 

PIA-3A  control 

PIA-3A  control 

2A 

PIA-3B  data/direction 

PIA-3B  data/direction 

2B 

PIA-3B  control 

PIA-3B  control 
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pulse  (see  Table  3-2)  out  the  experiment  interface  panel  (Trigger  A)  by 
the  write  instruction  "POKE  50452,0"  from  BASIC,  "STA  $C514"  from  assembly 
language,  or  "C514:0"  from  the  Apple  monitor.    The  KIM  can  send  the  same 
trigger  pulse  by  "STA  $F014"  from  assembly  language  only.    The  Apple  or 
KIM  must  have  control  of  the  S-100  bus  before  the  addresses  in  Table  3.2 
can  be  used  by  the  corresponding  microcomputer.    More  detail  on  the 
function  and  use  of  the  devices  in  the  tables  follows  in  the  remainder  of 
this  chapter.  Chapter  4,  and  the  appendices. 

The  Apple  to  S-100  control  logic  shown  in  Figure  3.3  is  mounted  on  a 
hobby/prototyping  board.    All  Apple  interface  functions  and  addresses  in 
Table  3.1  are  generated  by  this  board.    The  eight  write  activated  devices 
are  decoded  by  the  74LS138.    Since  the  Apple  uses  the  device  select  line, 
DEVS,  the  board  may  be  operated  from  any  slot  to  generate  a  500  nsec 
active  low  pulse  an  the  appropriate  line.    Writing  to  the  S-100  request 
address,  DEVO,  pulses  the  SIR  line.    This  activates  the  S-100  bus  for  use 
by  the  Apple  and  places  the  KIM's  processor  in  a  hold  state.    Writing  to 
the  S-100  release  address,  DEVI,  pulses  the  sId  line  causing  the  S-100  bus 
to  return  to  the  KIM's  control  and  restarts  the  processor  where  it 
stopped . 

The  Apple  and  S-100  bus  have  a  common  8  K  of  RAM  addresses  from  AOOO 
to  BFFF.    Two  memories  exist  in  this  area,  one  in  the  Apple's  dynamic  RAM 
and  the  other  in  the  S-100's  static  RAM.    The  active  memory  is  selected  by 
writing  to  the  memory  request  and  disable  addresses.    Disabling  the  Apple 
8  K  memory  leaves  the  microcomputer  with  40  K  of  RAM  and  sets  the  Apple 
memory  enable  line,  AME,  to  a  logical  zero.    The  S-100  memory  enable  line, 
S1ME,  goes  to  a  logical  one  when  the  Apple  is  using  the  S-100  memory; 
however,  the  S-100  bus  must  be  requested  before  the  memory  is  available. 
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KSRl 
KSR2 


PHLDA 


Figure  3,3.     Schematic  diagram  of  the  Apple  to  S-100  control  logic  for 
device  address  decoding  and  status  ports. 
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Both  memories  may  be  disabled  from  use  by  the  Apple  but  a  network  of 
flip-flops,  AND,  and  NAND  gates  prevents  both  memories  from  being  used  at 
the  same  time.    This  logic  also  causes  one  memory  to  be  released  if  the 
other  is  requested. 

Writing  to  the  enable  and  disable  interface  addresses  DEV4  and  DEV5 
allows  the  selection  of  the  interface  reset  feature.    When  the  reset  is 
enabled,  pressing  the  reset  key  on  the  Apple  keyboard  will  force  the 
interface  to  a  configuration  where  the  S-100  memory  and  bus  is  released 
and  the  Apple  8  K  memory  is  activated.    Disabling  the  feature  lets  the 
Apple  reset  itself  without  changing  the  configuration  of  the  interface 
hardware.    This  is  very  useful  when  debugging  software  in  the  S-100  memory 
from  the  Apple. 

A  74154,  4  to  16  line  decoder,  generates  the  read  activated  device 
addresses.    Two  of  the  addresses  cause  a  reset  and  non-maskable  interrupt 
(NMI)  request  to  be  sent  to  the  KIM.    These  addresses  function  the  same  as 
if  someone  pushed  the  corresponding  buttons  on  the  KIM.    The  interface 
status  and  KIM  status  ports  are  supplied  by  the  74LS367  hex  tristate 
driver  chip.    The  interface  status  port  tells  the  configuration  of  the 
Apple  to  S-100  interface  through  the  four  most  significant  bits  of  the 
port.    Bit  7  represents  the  Apple  8  K  memory  status  and  is  a  one  when  the 
Apple  has  48  K  of  RAM  and  a  zero  when  there  is  only  40  K  of  Apple  memory. 
Bit  6  is  a  one  when  the  S-100  memory  is  requested;  a  zero  vrtien  it  is  not. 
Bit  5  indicates  the  state  of  the  reset  feature  and  is  a  one  when  it  is 
activated.    Bit  4  monitors  the  S-100's  processor  hold  available  line, 
PHLDA,  and  is  a  one  when  the  bus  is  ready  for  use  by  the  Apple.    The  four 
least  significant  bits  of  the  port  are  not  connected  and  have  no  meaning. 


33 

The  8  K  bytes  of  memory  which  overlap  between  the  Apple  and  S-100 
systems  are  decoded  by  the  logic  in  Figure  3.^.    The  three  most 
significant  bits  of  the  1 6-bit  address  bus  are  decoded  for  the  range  AOOO 
to  BFFF.    This  signal  is  then  used  to  disable  the  memory  decoder  lines 
F1-8  and  D1-8  from  the  Apple  memory  select  headers  without  destroying  the 
contents  of  the  memory.    The  Apple  memory  enable  line  gates  the  memory 
select  line  "on"  when  necessary.    The  BUSEN  line  activates  the  data  bus 
drivers  during  phase  two  of  the  system  clock  when  S-100  memory  or  devices 
are  used  by  the  Apple.    This  coordinates  the  data  bus  so  it  does  not 
interfere  with  the  memory  refresh  cycles. 

The  seven  74LS245's  in  Figure  3.5  drive  the  data  and  address  buses 
between  the  Apple  and  S-100  bus.    The  three  chips  on  the  left  of  the 
figure  are  mounted  with  the  rest  of  the  Apple's  control  logic  on  the 
prototype  board,  normally  in  slot  5,  v*iich  has  just  been  described.  The 
right  four  are  on  the  wire  wrap  interface  board  occupying  one  slot  of  the 
KIMSI  S-100  bus.    When  the  S-100  bus  is  requested  the  PHLDA  line  goes  low 
and  activates  the  l6-bit  address.    This  bus  operates  in  only  one 
direction,  from  the  Apple  to  the  S-100  hardware.    The  bi-directional  data 
bus  from  the  Apple  is  driven  to  the  S-100  system  where  it  is  converted  to 
a  data  in  and  data  out  bus. 

KIM  and  S-100  Systems 

The  KIM  is  the  major  user  of  the  S-100  bus  and  its  peripherals  for 
experiments.    The  signal-to-noise  ratio  of  a  digitized  signal  is  enhanced 
by  signal  averaging.    However,  this  requires  a  very  reproducible  timing 
sequence  for  each  signal  recorded  and  is  especially  critical  in  the  fticr 
experiment.    The  necessary  synchronization  for  the  experiment  is 
accomplished  by  a  8  MHz  system  master  clock  provided  by  the  frequency 
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synthesizer's  internal  clock.    This  signal  is  divided  by  a  clock  generator 
circuit  on  the  experiment  interface  which  then  supplies  a  1  MHz  clock 
directly  to  the  KIM's  CPU  phase  zero  input  from  which  all  computer  timing 
is  derived.    A  variable  sample  clock  rate  used  to  sample  the  waveform 
recorder  is  derived  from  the  master  clock  by  a  divide-by-N  circuit  on  the 
experiment  interface . 

Once  the  Apple  address  and  control  lines  are  driven  over  to  the  S-100 
system,  the  logic  in  Figure  3.6  converts  these  signals  to  those  needed  by 
the  KIM  and  S-100  bus.    Two  mono-stable  vibrators,  9602's,  stretch  out  the 
0.5  ysec  KIM  reset  and  MNI  pulses  from  the  Apple  into  2  usee  pulses.  The 
longer  pulse  is  not  important  for  the  edge  sensitive  NMI  signal  but  is 
imperative  for  the  level  sensitive  reset  line.    Open  collector  inverters 
drive  the  lines  in  a  wired-or  fashion  to  the  KIM's  CPU.        *  -  - 

The  S-100  request  pulse,  SIR,  causes  the  KIM's  ready  line  to  drop  to 
a  logical  zero,  thus  stopping  the  CPU  while  the  S-100  release  pulse,  SID, 
causes  it  to  restart.    The  pulses  are  latched  by  one  74LS74  flip-flop  and 
synchronized  with  the  KIM's  phase  two  clock  by  a  second  flip-flop  to 
generate  the  KIMSI  ready  and  hold  signals,  KHLD  and  KRDY.    This  insures 
that  the  transition  only  occurs  during  phase  one  of  the  system  clock  as 
required  by  the  6502  CPU  chip.    When  the  S-100  bus  is  free  of  the  KIM's 
control,  the  processor  hold  available  line,  PHLDA,  goes  high  and  is  driven 
by  a  pair  of  inverters  to  the  Apple  and  other  circuits. 

Additional  logic  gates  convert  the  Apple  control  signals  to  S-100 
type  signals.    The  most  widely  used  signals  are  the  memory  select  lines, 
MWRITE  and  SMEMR,  I/O  select  lines,  SINP  and  SOUT ,  and  read/write  lines, 
PBDIN,  Pm,  and  SWO.    Selection  of  the  proper  S-100  data  bus  driver,  data 


Figure  3.6.     Schematic  diagram  of  the  Apple  to  S-100  control  line 

decoding  logic.    The  upper  half  converts  the  Apple  inter 
face  signals  to  standard  S-100  signals.    The  center  cir- 
cuit coordinates  the  S-100  request  and  release  commnads 
with  the  KIM's  phase  two  clock.     The  circuits  at  the  hot 
torn  produce  the  KIM  reset  and  interrupt  signals. 
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in  enable  (DIEN)  and  data  out  enable  (DOEN),  is  determined  by  a  pair  of 
NAND  gates  which  combine  the  PHLDA  and  R/W  signals. 

The  Bionation  waveform  recorder  is  controlled  and  interrogated  for 
data  by  two  parallel  input  ports  and  one  output  port  listed  in 
Table  3.1.    The  control  port  uses  three  of  the  four  least  significant  bits 
to  arm,  request  output,  and  toggle  the  next  data  point  out  of  the 
recorder's  memory  using  bits  0,  2,  and  3,  respectively.    Bits  number  2,  4, 
5,  5,  and  7  are  not  used.    Bits  7  and  5  of  the  status  port  display  the 
data  ready  flag  and  record  flag,  respectively.    The  8-bit  data  are  read 
from  the  recorder's  memory  through  the  Biomation  data  port.  Data 
recording  is  started  by  arming  the  Biomation  then  triggering  it  from  the 
front  panel.    For  an  example  of  how  to  use  this  peripheral  check  the 
Biomation  operator's  manual  and  the  subroutines  INIT,  BREQ,  GETB,  and 
BIODAT  in  the  assembly  language  program  KIM  ICR-MO  TEXT  listed  in 
Appendix  B. 

Experiment  Interface 

Most  computer  systems  are  designed  for  controlling  only  a  limited 
number  of  instruments,  usually  one.    This  often  makes  the  alteration  of  an 
experiment  difficult  since  the  interfaces  were  designed  only  for  a 
particular  task.    In  many  cases  the  hardware  is  present  but  the 
interconnections  are  integrated  into  the  system  such  that  it  is  virtually 
impossible,  or  at  least  undesirable,  to  modify  it  further.    The  experiment 
interface  described  here  provides  easy  hook-up  to  many  devices  through 
standard  types  of  connectors,  while  expansion  capabilities  are  present  in 
the  form  of  unused  gates  and  decoded  addresses. 

A  block  diagram  of  the  interface  is  shown  in  Figure  3.7.    A  19  in. 
rack  mount  front  panel  houses  a  series  of  BNC  connectors,  banana  jacks, 
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buttons,  and  a  switch.    The  chassis  of  the  interface  holds  a  row  of  eight 
44-pin  edge  connectors.    Five  of  the  slots  are  in  use  leaving  three  ready 
for  expansion  with  the  power  lines  and  data  bus  connected.    Starting  on 
the  right  while  looking  from  the  back,  the  first  board  houses  several 
miscellaneous  circuits.    The  second  board  generates  the  time  base  for  the 
computer  and  waveform  recorder.    The  third  holds  two  8-bit  D/A's,  the 
fourth  supplies  the  8-bit  A/D,  and  the  fifth  contains  three  peripheral 
interface  adapters  (PIA's).    All  boards  are  inserted  such  that  their 
components  face  to  the  left  when  seen  from  the  rear,  except  for  for  the 
first  board  due  to  the  large  heat  sinks  on  the  voltage  regulators. 

The  interface  control  logic  and  bus  drivers  (see  Figures  3.8  and  3.9) 
convert  the  S-100  bus  signals  to  those  used  by  the  custom  experiment 
interface.    This  involves  the  integrating  of  the  two  S-100  mono-direction 
buses  into  a  single  bi-direction  bus.    The  S-100  I/O  control  lines,  SINP 
and  SCUT,  are  decoded  with  four  of  the  address  lines  to  generate  the 
active  low  signals  I/01X  and  I/02X.    The  four  least  significant  bits  of 
the  address  bus  and  five  S-100  timing  signals  are  driven  to  the  experiment 
interface  slots.    One  driver  bit  is  still  available  for  expansion  on  the 
74LS367. 

The  KIM  status  port  is  produced  by  a  74LS138  and  74LS74  shown  in 
Figure  3.10.    Eight  addresses  are  write  activated  by  the  decoder  but  only 
one  is  used,  I/OIF,  leaving  seven  available  for  expansion.    A  write  to  the 
flip-flops  latches  the  contents  of  bits  6  and  7  of  the  data  bus;  thus  four 
different  statuses  can  be  indicated.    A  read  of  this  address  has  no 
significance  and  returns  garbage . 

A  2  MHz  clock  runs  through  the  interface  cable  and  supplies  the 
corresponding  S-100  clock  not  provided  by  the  KIMS I .    For  operation  of  the 
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Figure  3.8.     Schematic  diagram  of  the  experiment  interface  data  bus 
enable  and  address  decoding  logic. 
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Figure  3.10.     Schematic  diagram  of  the  KIM  status  port  (upper  half) 
and  the  KIM  interrupt  request  generators  on  the 
experiment  interface. 
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PIA's  the  phase  two  clock,  02 »       selected  from  either  the  KIM  or 
Apple  depending  on  which  has  control  of  the  S-100  bus.    Two  momentary 
contact  switches  are  debounced  by  a  pair  of  NAND  gates.    One  is  connected 
to  the  KIM's  NMI  line  while  the  other  produces  an  interrupt  request.  The 
7404  inverter  acts  as  a  2  ysec  pulse  generator  to  shorten  the  relatively 
long  time  the  IRQ  button  is  pressed  by  hand.    On  the  experiment  interface 
front  panel  there  is  an  active  high  and  active  low  interrupt  request  bnc . 
A  5.1  K  pull  up  resistor  is  used  on  the  active  low  line  while  a  Darlington 
transistor  makes  an  active  high  line.    The  IRQ  lines  are  combined  in  a 
four  input  NAND  gate  to  form  one  pulse.    Both  the  NMI  and  IRQ  lines  are 
driven  by  7403  open  collector  gates  in  a  wired-or  fashion. 

Experiment  interface  slots 

Signals  from  the  experiment  interface  bus  driver  run  to  the  boards  on 
the  experiment  interface  main  frame .    The  first  board  supports  the 
interface  power  supplies,  device  decoders,  trigger  generator,  and  relay 
control.    The  circuits  for  these  functions  are  shown  in  Figure  3.11.  The 
interface  has  two  D.C.  power  supplies;  each  is  capable  of  furnishing 
1  amp  at  +5  v  (for  digital  circuits),  100  ma  at  +15  v,  and  100  ma  at 
-15  V  (for  operational  amplifiers).    One  supply,  on  the  first  board, 
handles  the  power  bus  for  the  interface  cards  while  the  other,  mounted  on 
the  chassis,  is  accessible  through  the  front  panel  for  prototype  designs. 

The  74LS138  device  decoders  generate  eight  read  and  eight  write 
activated  device  addresses  from  the  I/01X  line.    Six  of  the  write 
addresses  are  used,  leaving  two  for  expansion,  while  only  four  of  the  read 
addresses  are  in  use.    The  decoders  produce  a  0.5  ysec  active  low  pulse 
when  selected.    Two  of  the  write  generated  pulses  are  channelled  directly 
to  the  interface  front  panel  to  the  TRIGGER  A  and  TRIGGER  B  bnc 

connectors.    These  outputs  are  capable  of  driving  ten  TTL  loads. 
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Figure  3.11. 


Schematic  diagram  of  the  experiment  interface  device 
decoders  (upper  half)  and  the  relay  control  logic  (lower 
half) . 
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A  single  relay  with  two  parallel  contact  sets  is  accessible  through 
the  interface  front  panel.    The  contacts  are  isolated  from  the  interface 
chassis  and  can  control  up  to  120  v  A.C.  at  2  amps.    The  relay  logic  uses 
a  single  74LS74  flip-flop  to  regulate  the  power  to  the  relay  coil  through 
a  2N3904  transistor.    Reading  the  address  1/011  powers  the  relay  coils 
causing  the  normally  open  and  normally  closed  circuits  to  change  state. 
Reading  address  1/010  returns  the  contacts  to  their  normal  states. 

The  system  clocks  are  generated  by  the  logic  shown  in  Figure  3.12. 
The  8  MHz  master  clock  from  the  frequency  synthesizer  is  first  divided  by 
a  74LS191  counter  into  4,  2,  1,  and  0.5  MHz.    The  4  MHz  clock  drives  the 
four  count  down  counters  synchronously  to  produce  the  variable  time  base. 
Two  74LS373  8-bit  latches  hold  the  l6-bit  clock  divide  value.    If  the 
clock  is  loaded  with  zero,  then  the  dividing  stops  and  the  output  stays  a 
logical  zero.    The  clock  may  be  divided  by  any  preset  value  from  1  to 
65535  to  generate  a  4  MHz  to  61.036  Hz  active  low  pulse  lasting  125  nsec. 
This  pulse  was  intended  primarily  to  serve  as  an  external  clock  for  the 
waveform  recorder,  however,  due  to  limitations  of  the  recorder; not  all  the 
times  are  usable.    Reading  from  address  1/012  causes,  at  the  end  of  the 
read  cycle,  the  clock  to  reset  itself  and  start  counting  over  again.  A 
series  of  AND  gates  and  inverters  act  as  fast  look  ahead  logic,  insuring 
reliable  operation  at  4  MHz.    All  clock  signals  are  driven  from  the  card 
by  either  a  74LS367  line  driver  or  inverter.    The  inverted  signals  are  the 
500  KHz  to  the  A/D  and  the  2~MHz  to  the  S-100  bus.    All  the  non-inverted 
clocks  go  directly  to  their  corresponding  bnc  connector  on  the  front 
panel. 


Figure  3.12,     Schematic  diagram  of  the  experiment  interface  clock 
generator  logic. 
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The  two  D/A's  and  their  associated  circuits  are  shovm  in 
Figure  3.13.    The  1408  8-bit  converters  generate  a  current  which  is 
changed  to  a  voltage  by  an  operational  amplifier.    Two  pots  for  each  op 
amp  adjust  the  output  voltage  zero  and  range.    The  voltage  varies  linearly 
with  the  input  from  -5.00  v  for  00  to  4.96  v  for  FF.    The  output  is 
smoothed  by  a  200  KHz  low  pass  filter.    Two  8-bit  latches  hold  the  input 
value  so  that  the  output  voltage  holds  steady  between  conversions. 

Analog  to  digital  conversion  is  performed  by  a  Datel  ADC-EH8B2  fast 
8-bit  converter  shown  in  Figure  3-14.    The  converter  is  triggered  by  the 
500  kHz  clock  and  latched  with  the  A/D's  end  of  conversion  signal.  This 
keeps  the  digitized  signal  synchronized  with  the  KIM.    The  Apple  is  not 
synchronized  with  the  A/D  converter;  therefore,  one  should  exercise  care 
when  using  them  together.    Bipolar  analog  signals  from  -5-00  to  4.96  v 
change  the  digitized  output  from  00  to  FF,  respectively.    Input  voltage 
overload  protection  is  accomplished  by  a  pair  of  back  to  back  18  v  zener 
diodes.    The  input  impedance  of  the  A/D  converter  is  4.45  Kohm  and  should 
be  buffered  when  connected  to  any  high  impedance  voltage  source. 

Three  PIA's,  shown  in  Figure  3-15,  supply  48  bits  of  parallel  I/O  and 
12  control  lines.    The  output  of  PIA-2  runs  to  the  front  panel  parallel 
output  BNC's  to  produce  an  A  and  B  port  with  their  corresponding  control 
lines.    PIA's  1  and  3  go  to  a  pair  of  DB-25  connectors  on  the  back  of  the 
chassis.    The  bits  represent  a  maximum  of  one  TTL  load  when  used  as 
inputs.    As  outputs  the  A  port  can  drive  one  TTL  load  v*iile  the  B  port  can 
directly  drive  the  base  of  a  transistor  switch.    The  PIA  interrupt  lines 
are  connected  in  a  wired-or  fashion  to  a  single  interrupt  line;  therefore, 
polling  must  be  used  if  more  than  one  chip  is  interrupt  active.    The  6820 
and  6821 's  PIA's  are  latched  with  the  phase  two  clock  and  operate  in  the 
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Figure  3.15.     Schematic  diagram  of  the  three  peripheral  interface 
adapters  on  the  experiment  interface  panel. 
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I/02X  address  range.    For  programming  considerations  one  should  consult 
the  reference  manual  [68] .    Sample  PIA  initialization  procedures  are  given 
in  subroutine  INIT  of  the  KIM  ICR-MO  TEXT  assembly  language  program  listed 
in  Appendix  B. 

Connections  to  the  marginal  oscillator 

The  pulse  sequence  for  the  icr  spectrometer  is  channelled  through 
parallel  port  B  on  the  front  panel  of  the  experiment  interface.    Six  of 
the  eight  bits  are  presently  used,  leaving  two  available  for  expansion 
purposes  such  as  triggering  a  dye  laser  or  gating  a  CO2  laser.  The 
pulse  sequence  needed  by  the  icr  and  m.o.  detector  is  shown  in 
Figure  3.16  along  with  the  corresponding  bit  supplying  the  pulse.  The 
grid  pulse  is  sent  to  the  EXT.  GRID  bnc  connector  on  the  icr  control 
panel.    The  EXT. /INT.  switch  is  then  set  to  external  allowing  the  computer 
to  gate  the  electron  beam.    The         line  is  connected  to  the 
external  trigger  of  the         oscillator.    The  quench  pulse  is  spliced 
into  the  icr  cell  plate  potential  control  electronics  through  the  quench 
pulse  bnc  on  the  controller  front  panel  after  disconnecting  the 
controller's  internal  clock  from  the  circuit.    The  reset,  detect,  and 
sample  and  hold  (S/H)  pulses  are  sent  to  the  m.o.  detector  and  post 
detection  system  with  the  timing  shown. 

Once  the  m.o.  detector  output  is  connected  to  the  input  of  the 
waveform  recorder  the  system  is  ready  to  go,  provided  the  software  is  all 
written. 
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Figure  3.16.     Pulse  timing  sequence  for  a  typical  icr  experiment  using 
our  marginal  oscillator  detector  and  cell  potential  con- 
trol electronics.     The  designated  bits  of  PIA-2B  on  the 
right  control  the  corresponding  signals  on  the  left. 


CHAPTER  4 
SOFTWARE 

While  hardware  gives  a  computer  its  physical  capabilities  it  is  the 
software  which  bestows  versatility.    This  chapter  discusses  programs 
written  to  operate  the  icr  spectrometer  in  our  laboratory  with  the 
microcomputer  system  previously  described.    These  programs  represent 
only  one  possible  application  of  the  computer  system.    The  full 
capabilities  of  the  hardware  are  not  even  close  to  being  thoroughly 
taxed,  and  thus  will  require  additional  software  in  order  to  satisfy  the 
needs  of  an  innovative  scientist  or  curious  student. 

The  subsequent  description  outlines  in  some  detail  an  operational 
and  functional  picture  of  the  software.    Various  tasks  are  examined  such 
as  initialization  and  control  of  the  hardware  and  loading  of  BASIC  and 
machine  language  code,  along  with  data  acquisition,  storage,  and 
analysis.    The  coordination  of  the  master  and  slave  computers  through 
software  instructions  is  also  explained. 

Program  Initialization 
The  program  is  started  by  inserting  the  "MO  Program  Disk"  and 
typing  the  lines 

<reset> 

C0D7:0 

6<ctrl-P> 

where  the  <  >  signifies  pressing  a  particular  key  on  the  keyboard. 
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otherwise  a  series  of  keys  is  pressed.    All  entry  lines  except  reset  end 
with  a  RETURN  key.    This  sequence  of  lines  resets  the  machine  to  monitor 
control,  configures  the  Apple  to  40  K  of  RAM  memory,  and  boots  the 
disk.    Under  control  of  the  "HELLO"  program  the  language  APPLEFAST  is 
loaded  into  the  language  card,  if  present,  and  an  EXEC  file  containing 
the  commands 

POKE  103,1 
POKE  104,64 
POKE  16384,0 
NEW 

POKE  232,0 

POKE  233,10 

BLOAD  ICR  SHAPE, A$AOO 

BLOAD  APPLE  ICR-MO  CODE 

RUN  ICR-MO 

is  started.    The  first  four  set  the  system  pointers  such  that  the  BASIC 
program  loads  in  above  page  one  of  the  high  resolution  graphics.  This 
overrides  the  normal  load  and  frees  memory  below  the  graphics  for  the 
Apple  machine  language  code  and  shape  table.    The  next  three  commands 
identify  to  BASIC  the  location  of  the  shape  table,  a  plus  sign  used  as  a 
data  display  cursor,  and  load  it  into  memory.    The  last  two  commands 
load  the  Apple  machine  code  for  communicating  with  the  KIM  and  then 
begin  loading  and  executing  the  BASIC  control  program,  ICR-MO. 

Functional  Description  of  Program 
The  operation  of  a  program  by  a  garden  variety  homo  sapien  does  not 
require  understanding  the  internal  workings  of  the  routine,  thus  a 
functional  description  of  the  program  is  sufficient  for  performing 
experiments  solely  to  obtain  chemical  information.    For  a  detailed 
description  of  how  the  programs  work  one  should  see  the  appendices. 
Appendix  A  explains  the  BASIC  program  used,  while  Appendix  B  covers  the 
assembly  language  code.    The  following  text  provides  a  brief  description 
of  the  commands  and  concepts  necessary  to  run  an  experiment. 
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Numberpfi  Command  Mode 

After  initializing  the  BASIC  program,  a  set  of  instructions  for 
initializing  the  KIM  appears  on  the  video  monitor.    These  instructions 
must  be  entered  by  hand  through  the  KIM's  keyboard.    They  disable  the 
processor's  interrupt  request  service  flag  through  its  status  port  and 
set  the  non-maskable  interrupt  request  vector  to  the  entry  point  of  the 
KIM's  experiment  code. 

Next  a  list  of  numbered  commands,  as  seen  in  Table  4.1,  appear  on 
the  video  screen  and  a  request  for  a  command  number  is  made.    Any  one  of 
the  ten  numbers  may  be  typed  at  this  time  in  order  to  perform  the 
specified  task.    An  invalid  entry  simply  causes  the  computer  to  ask  for 
another  command.    A  description  of  the  numbered  commands  follows, 
however,  not  in  the  order  given  in  the  table. 

The  EXIT  PROGRAM  (1)  command  provides  a  graceful  way  to  leave  the 
program  and  return  to  the  operating  system  for  control.  Inadvertent 
breaks  from  the  program  often  can  be  forgiven  by  typing  the  string  "GO 
TO  8000".    This  places  one  back  at  the  numbered  command  display  and 
entry  without  resetting  the  program  variables.    When  all  else  fails  the 
disk  should  be  rebooted. 

Occasional  hardware  glitches  in  the  KIM  and  S-100  system  can  be 
corrected  by  reloading  KIM's  program  with  the  LOAD  KIM  PROGRAM  (10) 
command.    This  restarts  the  same  routine  described  earlier  for 
initializing  the  KIM. 

The  MASS  CALCULATION  (6)  command  allows  determination  of  any  one  of 
the  three  icr  resonance  parameters  when  the  other  two  are  known.  The 
resonance  frequency,  magnetic  field  strength,  or  mass  to  charge  ratio 
are  determined,  neglecting  trapping  potential  effects,  by  Equation  2.1. 
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Table  4.1,    List  of  numbered  commands. 

Nvunbgr  Command 

1  EXIT  PROGRAM 

2  STORE  DATA  ON  DISK 

3  CATALOG  DISK 

4  UPDATE  EXPERIMENT  CONDITIONS 

5  CHANGE  PARAMETERS 

6  MASS  CALCULATION 

7  LOAD  PARAMETER  FILE 

8  SAVE  PARAMETER  FILE 

9  PRINT  DATA 

10  LOAD  KIM  PROGRAM 


60 

An  icr  experiment  is  set  up  and  executed  by  the  program  sequence 
invoked  by  the  command  CHANGE  PARAMETERS  (5).    At  this  time  the 
experiment  mode  of  operation  is  entered.    Due  to  its  significance,  this 
operation  mode  will  be  discussed  later. 

The  command  of  UPDATE  EXPERD^NT  CONDITIONS  (4)  allows  the 
recording  of  many  of  the  variables  not  controlled  by  the  computer  but 
useful  in  interpreting  the  data.    The  conditions  recorded  are  listed  in 
the  lower  half  of  Table  4.2  and  can  be  set  to  any  integer,  real,  or 
scientific  notation  value.    During  this  routine  four  key  commands  are 
operational  and  are  given  in  Table  4.3B.    These  key  commands  function  as 
explained  in  the  section  on  the  experiment  mode  of  operation. 

After  an  experiment  has  been  completed,  the  data  may  either  be 
saved  on  disk  or  listed  on  a  printer  by  the  appropriate  numbered 
commands  (2  and  9  respectively).    Data  is  stored  on  disk  in  a  text  file 
whose  format,  first  the  total  number  of  values,  then  a  list  of  X,Y 
pairs,  is  compatible  with  the  scientific  software  package  [69]. 
Depending  on  the  experiment  performed,  determined  by  the  scan  modes 
listed  in  Table  4.3C,  the  pairs  represent  different  quantities.  Their 
exact  significance  will  become  apparent  when  the  scan  modes  are 
described  later.    For  experiment  scan  modes  1  and  2  the  X,Y  pairs 
correspond  to  the  point  number  and  m.o.  output  level,  respectively, 
while  for  scan  modes  3  and  4  the  pairs  correspond  to  the  detect  delay 
time,  in  milliseconds,  and  the  m.o.  output  level.    After  saving  on  disk 
under  any  given  filename  the  program  prints,  if  desired,  the  filename, 
date,  time,  description,  experiment  parameters,  and  experiment 
conditions.    Printed  data,  again  as  X,Y  pairs,  is  preceded  by  the  same 
list  Just  described  but  with  the  filename  "NONE." 
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Table  4.2.  List  of  experiment  parameters  and  conditions 
along  with  typical  values. 


EXPERIMENT  PARAMETERS 


GRID  PULSE  WIDTH 

5 

MSEC 

DETECT  PULSE  DELAY 

250 

MSEC 

STEP 

10 

MSFC 

1  ilaJ  i-J  \y 

WIDTH 

3 

MSEC 

W2  PULSE  DELAY 

0 

MSEC 

WIDTH 

5 

MSEC 

QUENCH  PULSE  DELAY 

265 

MSEC 

WIDTH 

5 

MSEC 

DUTY  CYCLE  TIME 

270 

MSEC 

SCAN  MODE 

4 

NUMBER  OF  PASSES 

10 

EXPERIMENT  CONDITIONS 

ELECTRON  BEAM  : 

GRID  POTENTIAL 

-75 

VOLT 

FILAMENT  POTENTIAL 

-70 

VOLT 

COLLECTOR  POTENTIAL 

12 

VOLT 

COLLECTOR  CURRENT 

50 

NAMP 

CELL  POTENTIALS:  UPPER 

-1  .25 

VOLT 

LOWER 

-1  .25 

VOLT 

SIDE(F) 

1.5 

VOLT 

SIDE(C) 

1.5 

VOLT 

END1 

-1.25 

VOLT 

END2 

-1.25 

VOLT 

MAGNETIC  FIELD  STRENGTH 

12.1 

KG 

TEMPERATURE 

50 

C 

R.F,  LEVEL 

30 

MV 

W1  FREQUENCY 

87.4 

KHZ 

W2  FREQUENCY 

0 

KHZ 

PRESSURE 

1.5E-06 

TORR 
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Table  4.3.  List  of  keyboard  commands  and  scan  modes  for 
the  CHANGE  PARAMETERS  and  UPDATE  EXPERIMENT 
CONDITIONS  routines. 


A)  KEY^  FUNCTION 

ctrl-A    Abort 

ctrl-C    Cursor 

ctrl-E    Erase 

ctrl-G    Graphics 

ctrl-K    Rate  Constant 

ctrl-P    Pulse 

ctrl-Q    Quit 

ctrl-T    Text 

+    Increment  value 

  Decrement  value 

B)  KRY*^  FUNCTION 

 >    Advance  cursor 

<    Back-up  cursor 

RETURN    Change  value 

ESC    Escape  to  numbered  commands 


C)  Scan  Mndes^ 

1-  Continuous  scan 

2-  Signal  averaged  scan 

3-  Detect  delay  sweep 

Detect  delay  sweep  with  baseline  correction 


a)  For  use  in  CHANGE  PARAMETERS  routine  only. 

b)  For  use  in  CHANGE  PARAMETERS  and  UPDATE  EXPERIMENT 
CONDITIONS  routines. 
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The  experiment  parameters  and  conditions  can  be  saved  on  disk,  as  a 
text  file,  for  any  particular  set  up  desired  by  the  LOAD  PARAMETER  FILE 
(8)  command.    Any  disk  filename  may  be  specified  but  if  none  is  entered, 
by  just  pressing  the  RETURN  key,  then  the  default  filename,  EXP 
PARAMETERS,  is  used.    This  file  is  loaded  automatically  during  program 
initialization  to  put  the  program  into  a  non-zero  state.    The  LOAD 
PARAMETER  FILE  (7)  command  loads  any  parameter  file  generated  earlier  by 
the  save  command.    Again,  if  no  filename  is  given  the  default  value  is 
used. 

Parameter  filenames,  data  filenames,  or  any  information  from  the 
disk's  directory  may  be  examined  by  the  CATALOG  DISK  (3)  command.  This 
generates  a  catalog  of  the  active  disk  on  the  video  monitor. 

Experiment  Mode 

The  program  sequence  started  by  the  CHANGE  PARAMETERS  numbered 
command  provides  a  means  for  performing  pulsed  icr  data  acquisition 
along  vdth  some  data  analysis.    The  routine  starts  by  displaying  the 
list  of  experiment  parameters  shovm  in  the  top  half  of  Table  4.2.  A 
flashing  greater  than  sign  (>),  the  parameter  pointer,  appears  before 
one  of  the  numbers.    All  parameter  changes  occur  to  the  number  next  to 
the  pointer.    It  also  indicates  the  program  is  ready  for  the  single  key 
instructions  listed  in  Tables  4.3  A  and  B.    A  logical  key  stroke  for 
each  command  aids  in  their  memorization.    A  description  of  each  key  and 
its  corresponding  function  is  given  below. 

Abort  (ctrl-A).    Used  to  terminate  a  pulse  sequence  in  progress. 
This  provides  an  early  exit  from  a  pulse  sequence  while  leaving  the  data 
recorded  up  to  that  point  in  memory.    After  aborting  the  current  pulse 


sequence,  the  KIM  is  placed  in  a  continuous  operation  mode  even  though 
the  data  are  not  being  recorded  by  the  Apple. 

Cursor  (ctrl-C).    Allows  movement  of  a  cursor,  plus  sign,  across 
the  data  previously  recorded.    The  cursor  position  is  controlled  by  the 
two  hand  controllers:  paddle  zero  for  coarse  movement  and  paddle  one  for 
fine  movement.    Displayed  below  the  plot  are  the  cursor  position 
parameters:  point  number,  X- value ,  and  intensity. 

Erase  (ctrl-E).    Clears  the  high  resolution  screen,  regardless  of 
whether  the  graphics  is  currently  being  displayed  or  not. 

Graphics  (ctrl-G).    Switches  the  graphics  display  on.    Changes  may 
still  be  made  on  the  text  screen  and  parameters  but  they  are  not 
visible. 

Rate  constant  (ctrl-K).    Allows  determination  of  decay  rates  from 
data  previously  collected  during  a  detect  delay  sweep  (scan  modes  3  and 
4).    Calculations  are  performed  by  the  ratio  method  described  in  Chap- 
ter 2.    Specifying  the  first  and  last  data  points  to  be  used  in  the 
calculation  allows  evaluation  of  the  decay  constant  from  a  portion  of 
the  data.    Data  may  be  corrected  for  non-zero  decay  curve  tails  when 
determining  decay  constants  by  specifying  an  appropriate  offset  value. 

Pulse  (ctrl-P).    Initiates  a  pulse  sequence  with  the  designated 
parameters  and  switches  automatically  to  graphics  display.    Any  key 
command  may  be  entered  and  executed  during  a  pulse  sequence  except  for  a 
quit,  escape,  cursor,  rate  constant,  or  another  pulse  command.  When 
single  key  commands  are  given,  pulse  sequences  in  progress  are  delayed 
between  points  then  continued  with  the  updated  parameters,  if  any.  The 
parameter  pointer  stops  flashing  during  a  pulse  sequence. 
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Quit  (ctrl-Q).    Causes  an  exit  from  the  program  and  a  return  to  DOS 
keyboard  control.    This  key  operates  the  same  as  the  EXIT  PROGRAM 
command  (number  1  in  Table  4.1). 

Text  (ctrl-T).    Switches  the  screen  from  graphics  to  text  display. 
When  an  experiment  is  in  progress  the  text  remains  on  the  video  screen 
as  experiments  are  performed  and  plotted  on  the  high  resolution  graphics 
screen. 

Plus  (+) .    Increments  the  value  of  the  experiment  parameter 
currently  indicated  by  the  pointer.    When  the  upper  limit  is  reached  the 
parameter  remains  at  the  maximum  value.    The  maximum  value  is  determined 
by  the  limits  set  in  the  "experiment  parameter  limit  check"  routine  (see 
program  lines  4000  to  4l80  in  Appendix  A). 

Minus  (-).    Decrements  the  value  of  the  experiment  parameter 
currently  indicated  by  the  pointer.    When  the  lower  limit  is  reached  the 
parameter  remains  at  the  minimum  value.    The  lower  limit  is  determined 
the  same  way  as  described  in  the  previous  paragraph. 

Advance  ( — >) .    Shifts  the  parameter  pointer  down  one  line.  The 
pointer  returns  to  the  top  of  the  list  if  the  key  is  pressed  when  the 
pointer  is  at  the  bottom  of  the  list. 

Back-up  (< — ).    Shifts  the  parameter  pointer  up  one  line.  The 
pointer  returns  to  the  bottom  of  the  list  if  the  key  is  pressed  vhen  the 
pointer  is  at  the  top  of  the  list. 

Change  (Return) .    Modifies  the  value  next  to  the  parameter 
pointer.    Changes  are  checked  against  the  parameter's  upper  and  lower 
limits,  then  corrected  as  necessary. 

Escape  (ESC).    Causes  exit  from  the  experiment  mode  of  operation 
and  return  to  the  numbered  command  mode. 
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The  experiment  parameters  controlled  by  the  computer  are  the  pulse 
times,  duty  cycle,  scan  mode,  and  number  of  passes  (when  signal 
averaging).    Four  different  types  of  experiments  (scan  nodes)  can  be 
performed  by  the  computer  and  are  given  in  Table  4.3C.    The  experiment 
is  selected  by  the  scan  mode  value. 

Continuous  scan  (1).    The  pulse  sequence  repeats  with  all 
parameters  held  constant.    The  detect  pulse  delay  parameter  equals  the 
detect  delay  time.    The  maximum  number  of  points  on  the  screen  and  in 
memory  is  280.    The  screen  is  cleared  when  full  and  the  previous  data 
are  replaced  with  the  incoming  data.    This  scan  mode  is  most  commonly 
used  when  tuning  the  icr  for  a  particular  ion. 

Signal  averaged  scan  (2).    Experiments,  such  as  magnetic  field  and 
double  resonance  sweeps,  are  signal  averaged  for  the  specified  number  of 
passes  during  one  sweep,  then  displayed  on  the  screen.    The  screen  is 
filled  with  280  points  representing  280  experiments  averaged  N  times 
each  (where  N  is  the  number  of  passes  averaged)  then  the  pulse  execution 
is  terminated  and  data  saved  in  memory. 

'  "  *    Detect  delay  sweep  (3).  The  grid  pulse  is  enabled  and  the  m.o. 
output  signal  is  averaged  for  the  specified  number  of  passes  at  a  given 
detect  pulse  delay  time.    This  detect  pulse  delay  time  is  swept  from 
zero  to  the  maximum  detect  pulse  delay  time  displayed  on  the  screen  in 
increments  of  the  displayed  step  size. 

Detect  delay  sweep  with  baseline  correction  (4).    Ion  signal  and 
baseline  are  averaged,  in  turn,  for  the  specified  number  of  passes.  The 
detect  delay  is  swept  as  in  scan  mode  3-    The  baseline  is  linear  least 
squares  fitted,  then  subtracted  from  the  ion  signal. 
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The  grid  pulse  starts  at  time  zero  and  can  vary  in  width  from  0  to 
25  msec.    Its  value  may  be  incremented  by  0.5  msec  steps  within  the 
above  limits  using  the  +  and  -  key  commands. 

The  detect  pulse  is  controlled  by  three  times:  detect  width,  delay, 
and  step.    The  width  can  vary  from  1.5  to  25  msec  in  0.5  msec  steps. 
The  detect  delay  and  step  size  have  two  different  meanings  depending  on 
the  scan  mode.    If  a  continuous  or  signal  averaged  scan,  modes  1  and  2, 
is  selected,  then  the  screen  shows  the  delay  time  used  in  the 
experiment .    This  delay  time  can  vary  from  zero  to  the  quench  delay 
minus  the  detect  pulse  width.    The  step  value  is  used  as  the  amount  of 
increase  or  decrease  when  the  detect  delay  time  is  changed  by  the  +  or  - 
keys.    During  a  detect  delay  sweep,  modes  3  and  4,  the  actual  detect 
delay  time  varies  from  zero  to  the  detect  pulse  delay  time  displayed  on 
the  screen.    The  detect  pulse  delay  time  is  swept  between  the  above 
limits  by  incrementing  its  value  by  the  detect  pulse  step  size. 

The        pulse  width  and  delay  times  are  controlled  by  their 
corresponding  parameters.    The  values  are  not  bound  by  any  limits.  The 
quench  pulse  width  can  vary  from  0  to  10  msec  while  the  quench  pulse 
delay  must  be  less  than  or  equal  to  the  duty  cycle  time  minus  the  quench 
pulse  width.    The  duty  cycle  time  determines  the  maximum  time  for  pulse 
generation  and  must  be  greater  than  or  equal  to  the  quench  pulse  width 
plus  quench  pulse  delay.    The  maximum  duty  cycle  time  is  currently  set 
at  16  seconds. 

When  a  signal  averaging  scan  is  in  progress,  the  number  of  passes 
taken  for  each  point  is  set  by  the  corresponding  parameter  and  may  vary 
from  1  to  255. 


Pulse  Timing  Parameters 

An  8-bit  pulse  routine,  whose  operation  is  described  in  Appendix  B, 
provides  the  Apple  or  KIM  microcomputer  with  the  ability  to  produce 
accurate  and  variable  pulse  sequences  needed  for  various  experiments.  A 
single  8-bit  parallel  port  produces  eight  separate  pulses  which  are 
independent  of  each  other  but  change  sychronously. 

The  pulses  generated  are  active  high,  that  is,  a  logical  one  when 
"on"  and  a  logical  zero  when  "off".    Each  pulse  corresponds  to  a 
particular  bit  in  the  pulse  hardware  port.    Two  parameters  are  needed 
for  each  pulse  bit,  an  "on"  time  and  an  "off"  time;  therefore,  16 
parameters  are  needed  to  specify  pulses  for  all  eight  bits.  One 
additional  parameter,  the  duty  cycle  time,  is  needed  to  specify  the 
maximum  time  for  pulse  generation.    All  parameters  are  16  bits  in  length 
and  thus  require  double  precision  values  (two  8-bit  bytes). 

The  pulse  parameters  mentioned  above  are  in  units  of  0.5  msec; 
consequently,  the  shortest  pulse  lasts  0,5  msec,  the  longest  lasts 
32568  msec.    Any  pulse  size  between  these  limits  is  possible  but  only  in 
0.5  msec  increments. 

Several  different  signal  types  can  be  generated  by  this  routine.  A 
typical  pulse  is  when  the  bit  goes  high  then  low  some  time  later.  The 
pulse  may  go  high  immediately  then  low  as  in  Figure  4.1A,  or  first  delay 
and  then  go  "on"  and  "off"  as  in  Figure  4. IB.    For  these  pulses  the  "on" 
time  is  greater  than  or  equal  to  zero,  and  the  "off"  time  is  greater 
than  the  "on"  time  but  less  than  or  equal  to  the  duty  cycle  time.  This 
will  generate  a  pulse  equal  in  length  to  the  difference  between  "on"  and 
"off"  times.    The  pulse  is  generated  each  time  the  program  is  executed. 
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Figure  4.1.     Some  of  the  pulse  shapes  which  can  be  produced  by  the 
8-bit  pulse  routine.     The  bits  are  switched  between  a 
logical  1  and  0  by  setting  the  "on"  and  "off"  times 
to  the  appropriate  values;    A)  "on"  time  equal  to  zero 
and  the  "off"  time  less  than  the  duty  cycle;  B)  "on" 
time  greater  than  zero  and  the  "off"  time  less  than  the 
duty  cycle;  C)  "on"  time  greater  than  the  "off  time"; 
D)  "on"  time  greater  than  zero  and  the  "off"  time  greater 
than  the  duty  cycle;  E)  "on"  time  equal  to  zero  and  the 
"off"  time  greater  than  the  duty  cycle. 
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A  pulse  port  bit  can  be  always  "off"  as  in  Figure  4.1C.    This  can 
be  generated  in  three  different  ways:  "on"  time  equal  to  "off"  time, 
"off"  time  less  than  "on"  time,  or  "on"  time  greater  than  duty  cycle 
time.    The  preferred  way  is  to  set  the  "off"  time  equal  to  zero,  then 
the  "on"  time  can  be  any  value.  ,  v; 

The  pulse  port  bit  may  remain  "on"  at  end  of  the  duty  cycle  as  in 
Figure  4. ID.    This  is  accomplished  by  making  the  "on"  time  less  than  or 
equal  to  the  duty  cycle  time,  while  the  "off"  time  is  greater  than  the 
duty  cycle  time.    A  pulse  bit  can  be  kept  high  at  all  time,  as  in  Figure 
4. IE,  by  setting  the  "on"  time  equal  to  zero  and  the  "off"  time  greater 
than  the  duty  cycle.    After  the  first  pass  through  the  routine  the  bit 
will  go  high  and  remain  in  that  state. 

KIM  Status  DisDlav 

Since  a  computer  which  is  working  properly  often  looks  like  one 
which  is  not,  the  slave  microcomputer  is  instructed  to  flash  its  vital 
signs  periodically.    The  KIM  displays  its  status  in  three  bytes  on  the 
KIM'S  6-digit  LED  front  panel  display.    The  first  byte  on  the  left 
displays  the  scan  mode  status  of  the  computer.    There  are  three  possible 
values:    00  for  end  of  signal  averaged  scan,  01  for  continuous  scan,  and 
02  for  signal  averaging  in  progress.    The  display's  middle  byte  changes 
value  by  one  each  time  a  pass  through  the  KIM  ICR-MO  program  is 
completed.    When  a  continuous  scan  is  in  progress  the  display  value 
increments,  but  during  an  averaging  scan  the  value  begins  at  the  number 
of  passes  then  decreases  to  zero.    The  right  byte  displays  the  data  read 
from  the  8-bit  A/D  during  the  last  pass.    A  hexadecimal  00  is  the  lowest 
value  and  FF  is  the  highest. 
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Data  Analysis 

Most  of  the  data  are  analyzed  through  the  use  of  three  commercial 
scientific  routines  when  possible  [69].    A  curve  fitter  is  used  to 
perform  linear  least  squares  calculations,  while  a  scientific  plotter 
enables  the  combining  of  several  different  data  sets  for  visual 
comparison.    The  program  VIDICHART  is  particularly  useful  when  comparing 
mass  spectra  to  each  other.    For  specific  details  about  operating  these 
routines  one  should  consult  the  manuals.    Hard  copy  of  plotted  data  is 
obtained  from  a  dot  matrix  printer  controlled  by  a  graphics  screen  dump 
routine  [70]. 


CHAPTER  5 
INSTRUMENT  CALIBRATION 


Introduction 

Due  to  the  low  gas  densities  involved  in  icr  studies,  the  most 
difficult  quantity  to  measure,  if  accurate  ion-molecule  rate  constants 
are  to  be  determined,  is  pressure.  All  experiments  were  carried  out  with 
gas  pressures  in  the  10"^  to  10"^  torr  range,  thus  allowing  the  use  of 
an  ionization  gauge  for  pressure  measurement.    Since  this  transducer 
works  by  ionizing  the  gas  and  then  measuring  the  resulting  ion  current, 
the  gauge  sensitivity  depends  upon  gas  characteristics,  particularly  the 
cross-section  for  ionization.  Unfortunately,  the  sensitivity  is  not 
directly  related  to  some  known  quantity  such  as  the  ionization  potential 
or  polarizability  of  the  gas;  thus,  the  gauge  must  be  calibrated  for 
each  sample  used.    If  the  ion  trapping  cell  is  not  near  the  ionization 
gauge,  as  in  our  instrument,  then  the  geometry  of  the  vacuum  system  must 
also  be  accounted  for  in  the  evaluation  of  the  true  pressure  by 
determining  a  system  correction  factor. 

If  one  lets       and  S  be  the  pressure  reading  and  sensitivity  of  the 
ionization  gauge,  and  F  the  system  correction  factor,  then  the  pressure 
of  the  gas  at  the  cell,  P^ ,  is  given  by  the  following  equation: 

Pq  =  P^  »  S  »  F  . 

The  sensitivity  can  be  determined  by  calibrating  the  ionization  gauge 
against  a  pressure  gauge  whose  output  is  independent  of  the  species 
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used,  such  as  a  capacitance  manometer.    The  system  correction  factor  is 
evaluated  by  comparing  reaction  rates  from  our  icr  spectrometer  with 
well-established  rates  for  a  known  reaction  and  then  adjusting  the 
results  by  the  appropriate  constant. 

One  of  the  most  highly  studied  chemical  systems  in  ion-molecule 
reaction  kinetics  must  be  the  ions  produced  in  methane  gas.    The  only 
reactions  which  occur  at  an  appreciable  rate  (k  >  10"^^  cm^/sec)  are 
those  of  the  primary  ions,  CH^  and  CH]J,  with  the  neutral  parent  [71]: 

CH]|[  +  CHi^  >  CH^  +  CH^  k^ 

and 

CH^  +  CB^   >  C^a^  +  H2  k2. 

The  rate  constants  measured  by  a  variety  of  methods  (such  as  drift  tube, 
icr,  and  other  mass  spectrometers)  for  these  two  reactions  are  in  good 
agreement  resulting  in  average  values  of       =  11  +  IX  10"^^  and 
k2  =  9  ±  1  X  10"'^  cm-'/sec  [72].    If  only  icr  results  are  considered  one 
obtains  averages  of  k^  =  1 1 .5  ±  0 .4  X  10"^^  and  k2  =  9.9  ±  1.6  X  10~^° 
cm^/sec  (all  errors  reported  are  95%  confidence  limits)  [73]. 

Due  to  the  abundance  of  earlier  work  on  ion-molecule  reactions  in 
methane  gas,  the  rate  constants  for  reactions  in  this  system  were  used 
for  calibrating  the  instrument.    The  following  text  describes  the 
procedures  used  for  determining  ionization  gauge  sensitivity,  system 
correction  factor,  and  reaction  rate  constants.    The  ion-molecule 
reaction  rates  for  several  other  chemical  systems,  reported  in  the 
literature,  also  are  determined  to  illustrate  the  validity  of  the 
corrections  applied. 
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Experimental 

Gauge  Sensitivity  Experiments 

The  ionization  gauge  was  calibrated  against  a  single-sided 
capacitance  manometer  [7^]  in  the  mid  10"^  to  10"^  torr  range.    To  aid 
in  reading  the  manometer,  a  digital  voltmeter  was  connected  to  the 
analog  output  of  the  manometer  control  unit.    Experiments  were  performed 
by  first  pumping  out  the  vacuum  chamber  to  a  low  background  pressure, 
then  reading  both  pressure  gauges.    Next,  a  sample  gas  was  leaked  into 
the  vacuum  system  followed  by  the  reading  of  both  gauges  again  when  the 
pressure  had  stabilized.    This  procedure  was  repeated  for  several 
different  sample  pressures  in  the  range  specified  above. 

Ion  Decav  Determination 

Ion  decay  experiments  were  performed  in  the  10"^  to  lO"*^  torr 
pressure  range.    Ions  were  formed  by  electron  impact  with  30  to  70  eV 
el-ectrons  and  trapped  in  the  icr  cell,  then  the  signal  intensity  was 
monitored  as  a  function  of  the  detect  pulse  delay  time  by  the  data 
system  previously  described.    The  cyclotron  frequency  for  most 
experiments  was  set  for  ion  resonance  at  1.2  T,  except  for  methane  which 
was  run  at  0.8  T  due  to  the  upper  frequency  limit  of  the  marginal 
oscillator.    For  a  typical  experiment  the  cell  trapping  plates  were  at 
+1  V,  the  upper,  lower,  and  end  plates  were  at  -1  v,  the  grid  was  pulsed 
for  5  ms,  and  a  4  ms  detect  pulse  gated  the  30  mv  peak-to-peak  marginal 
oscillator  irradiating  frequency.     After  each  experiment  ions  were 
removed  from  the  cell  with  a  5  ms  quench  pulse.     Ion  signals  were 
averaged  from  10  to  25  times  at  each  detect  pulse  delay  time,  then 
stored  on  disk.    Gaseous  samples  in  glass  bulbs  were  routinely  purified 
by  repeated  freeze-pump-thaw  cycles  in  our  experiments. 
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Results 

Gauge  Sensitivity  Experiments 

The  appropriate  background  readings  were  subtracted  from  the  sample 
pressure  for  both  the  ionization  and  manometer  readings  to  give  pressure 
changes.    Linear  variations  were  assumed  when  correcting  for  drifts  in 
the  capacitance  manometer  background  measurements.    Plots  of  the 
capacitance  manometer  vs  ionization  gauge  pressure  readings  for  a  few 
selected  gases  are  shown  in  Figure  5.1.    Slopes  of  the  lines  were 
determined  by  linear  least  squares  analysis  and  are  listed  under 
sensitivity  in  Table  5.1.    Linear  correlation  coefficients  for  all  fits 
were  better  than  0.996,  thus  indicating  good  linear  behavior  over  the 
pressure  range.  , 

Methane  Experiments 

Three  typical  ion  decay  curves  for  the  reaction  of  the  methane 
parent  ion  with  neutral  methane  gas  are  shown  in  Figure  5.2.    The  curves 
were  normalized  to  the  same  extrapolated  intensity  at  time  equal  to 
zero.    Plots  of  the  natural  logarithm  of  ion  intensity  versus  time 
yielded  straight  lines  verifying  the  application  of  pseudo-first  order 
reaction  kinetics.    Slopes  of  the  lines  were  determined  by  linear  least 
square  analysis .    When  the  pressure  read  from  the  ionization  gauge  is 
corrected  for  the  gauge  sensitivity  using  curves  as  in  Figure  5.1  to 
obtain  the  values  in  Table  5.1,  and  one  assumes  a  neutral  gas 
temperature  of  323  K  [75],  then  the  reaction  rate  constant  can  be 
calculated  by  the  method  described  at  the  end  of  Chapter  2.  The 
evaluation  of  the  CH^  reaction  rate  at  six  different  pressures  gives 
an  average  of  18.9  ±  1  .7  X  10"''°  cm^/sec.    Similar  analysis  of  four 
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CH^  ion  decay  signals  in  methane  results  in  an  average  of  16.8  ±  1  .2  X 
10"''°  cm^/sec. 

These  two  reaction  rates  are  obviously  much  higher  than  the  rates 

determined  by  icr  studies  described  earlier,  thus  indicating  a 

significant  system  correction  factor.    The  ratio  of  the  CH]J  reaction 

rate  to  the  CH^  rate  determined  in  this  study  is  1.13  and  compares  well 

with  the  ratio  of  1.16  determined  from  an  average  of  icr  experiments 

[73].    A  relative  rate  of  1.2  from  the  two  reactions  determined  on 

another  icr  is  in  good  agreement  with  our  results  [82].    If  one  assumes 

the  average  icr  literature  results  as  being  the  true  value  then  one 

obtains  a  multiplicative  system  correction  factor  of  1.64  ±  0.20  for  the 

CH]J  reaction  and  1.70  ±  0.40  for  the  CH^  results  or  an  average  of 

1.67  ±  0.21  (error  limits  were  determined  by  use  of  the  CH^  results). 

Using  this  value  one  calculates  a  corrected  rate  constant,  k  of 

'  corr' 

=  11.4  ±  2.5  X  10'"'°  cm3/sec  and  k2  =  1 0 . 1  ±  2 . 0  X  10"'' °  cm^/sec. 
Reactions  in  acetylene,  ethylene,  benzene,  and  chlorocyclopropane 
were  also  investigated  for  comparison  with  previous  studies.  Pressures 
read  from  the  ionization  gauge  were  corrected  for  gauge  sensitivity  to 
the  gas  and  then  the  reaction  rate  constants  were  determined  in  a  manner 
analogous  to  that  used  in  the  methane  system.    Using  the  system 
correction  factor  obtained  from  the  methane  calibration  reactions  one 
gets  the  rate  constants  k         in  Table  5.1. 
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PiscugsiQn 

Several  measured  quantities  involved  in  the  determination  of 
reaction  rate  constants  introduce  error  into  the  results.    One  important 
factor  is  reaction  time,  but  this  quantity  is  computer  controlled  and 
accurate  to  within  a  part  per  million,  thus  not  a  significant  source  of 
error.    Another  factor  is  the  resolution  of  the  digitizer.    In  this 
system  the  8-bit  A/D  can  distinguish  changes  of  0.4?  in  signal 
intensity,  a  value  much  smaller  than  the  noise  in  the  observed  ion 
signal.    As  can  be  seen  from  the  plots  in  Figure  5.2,  a  few  percent 
precision  can  be  obtained  by  signal  averaging.    Pressure  measurement 
inaccuracies  are  by  far  the  most  significant  source  of  errors  due  in 
part  to  the  inability  to  measure  accurately  the  pressure  and  also  due  to 
pressure  fluctuations  caused  by  pumping  speed  changes,  leak  valve 
variations,  and  gas  resevoir  pressure  changes.    Another  possible  source 
of  error  is  the  sensitivity  of  the  ionization  gauge  measured  in  the 
10"^  to  10"^  torr  pressure  range  which  must  be  assumed  the  same  for  the 
10"^  torr  range. 

The  system  correction  factor  of  1.67  ±  0.21  seems  reasonable  in 
light  of  the  location  of  the  ionization  gauge,  closer  to  the  oil 
diffusion  pump  than  the  cell,  and  the  gas  sample  inlet  to  the  high 
vacuum  chamber  being  directly  below  the  cell.    This  system  correction 
factor  indicates  a  higher  pressure  at  the  cell  than  at  the  gauge,  which 
is  consistent  with  the  facts  just  presented. 

Average  rate  constants,  k^Q^j,,  for  the  reactions  investigated  in 
this  study  are  presented  in  Table  5.1  with  95$  confidence  limits 
calculated  by  propagation  of  errors  from  the  pressure  correction  factor, 
ion  signal  decay  rate,  gauge  sensitivity,  and  uncertainty  in  the  methane 
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literature  values.    The  errors  are  such  that  all  values  overlap  for  the 
methane  system  indicating  good  agreement  with  the  literature.    In  the 
acetylene  system  the  reaction  of  the  parent  ion  with  the  neutral  forms  a 
condensation  product  which  can  fragment  into  various  ions  and  neutrals: 

C2H2  +   ^  Products 

As  before,  the  rate  of  k^  =  10.1  ±  2.5  X  10"^^  cm^/sec  compares  well 
with  the  literature  values.    A  similar  case  exists  for  the  reaction  of 
C2H2  and  Z^^^  ions  in  ethylene 

C2H2  +   >  Products  ki^ 

and 

CgH^  +  Z<^^  >  Products  k^  . 

Only  a  few  rates  have  been  reported  for  k|^  while  considerable  attention 
has  been  given  to  the  determination  of  k^ .    Again,  as  can  be  seen  in 
Table  5.1,  the  results  from  this  study  are  consistent  with  the 
literature  values. 

There  are  two  possible  mechanisms  for  removal  of  ions  from  the 
cell.    The  first,  which  supplies  chemical  information,  is  the  reaction 
of  the  ions  with  the  neutral  background  gas  to  generate  another  ion  and 
often  another  neutral.    The  second  mechanism,  ions  diffusing  to  the  cell 
walls  where  they  are  neutralized,  complicates  the  analysis  of  the  data. 
One  can  experimentally  determine  the  ion  loss  rate  by  looking  at  a 
non-reactive  or  very  slowly  reacting  chemical  system  so  the  diffusive 
losses  from  the  cell  are  dominant. 

The  parent  ion  in  benzene  is  known  to  react  slowly  to  form  a 
condensation  product  by 
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^6^6     ^6^6  ^  ^12^12  ^6 

with  a  rate  constant  of  1.2  X  10"^^  cm^/sec  for  ions  formed  by  30  eV 

electron  impact  [79] .    The  rate  constant  determined  may  be  complicated 

due  to  the  reformation  of  CgHg  ions  by  charge  transfer  from  fragments 

[83].    The  electron  impact  studies  on  benzene  showed  that  57%  of  the 

ions  formed  are  the  parent  [79] .    The  only  ions  which  can  transfer 

charge  to  the  neutral  benzene  are  C^jH^,  C^H2,  and  C2H2  (ionization 

potentials  were  taken  from  Rosenstock    et  al.  [84])  but  they  initially 

comprise  only  17$  of  the  total  ion  population  at  the  pressures  in  these 

experiments.    Our  experiments  with  benzene  at  a  pressure  of  1.2  X 

10"^  torr,  Figure  5.3,  show  that  7k%  of  the  peak  parent  ion  signal 

remains  500  ms  after  formation.    Assuming  no  loss  from  the  trapping  cell 

and  correcting  for  the  ionization  gauge  and  sytem  correction  factor, 

then  the  reaction  rate  constant  for  the  ion  decay  is  1.4  ±  1.1  X 
—  11 

10       cm-'/sec.    This  compares  well  with  the  rates  reported  in  the 
literature. 

The  slowest  ion  decay  signal,  also  shown  in  Figure  5.3,  is  observed 
for  the  reaction  of  cyclic  C^H^  with  propargyl  chloride: 

C^H^  +  CH=CCH2C1   >  products  k^. 

The  calculated  reaction  rate  from  this  study  is  k^  =  2.8  ±  1.0  X  10"^^ 
cm^/sec.    A  previous  study  of  this  ion  reported  that  no  reaction 
occurs  with  propargyl  chloride  [81].    The  value  obtained  does,  however, 
set  an  upper  limit  on  the  ion  loss  rate  which  in  turn  establishes  the 
ion  decay  detection  limit  for  this  instrument.    Since  almost  all 
experiments  involve  observation  of  the  ion  for  times  less  than  500  ms, 
cell  loss  may  be  neglected. 
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Conclusion 

The  trapping  efficiency  of  the  ion  cell  is  sufficient  to  justify 
the  neglect  of  non-reactive  ion  loss  processes  when  determining  most 
rate  constants  (k's  >  10"^^).    Sensitivity  of  the  ion  gauge  to  gas 
composition  is  linear  in  the  pressure  range  of  interest  but  pressure 
measurements  are  by  far  the  most  inaccurate  quantity  measured  and  must 
be  improved  to  increase  the  accuracy  of  the  reaction  rates  obtained. 
Our  instrument  appears  to  be  seeing  a  representative  picture  of  the 
reactions  of  interest  as  indicated  by  the  agreement  of  results  from  this 
study  with  previous  results. 


CHAPTER  6 
C^H*  ION  STUDIES 

Introduction 

Soot  Formation 

The  combustion  of  hydrocarbons  in  the  presence  of  oxygen  is  an 
important  process  presently  being  given  considerable  attention.  This 
process  results  in  the  liberation  of  heat  and  the  formation  of  various 
byproducts.    Under  certain  circumstances,  however,  byproducts  such  as 
soot  may  be  undesirable.    The  mechanism  responsible  for  the  formation  of 
soot  in  fuel  rich  flames  is  not  very  well  known  but  has  been  the  subject 
of  much  investigation.    A  recent  review  by  Calcote  discusses  this 
problem  in  detail  [85]. 

During  combustion,  the  fuel  and  oxidizer  are  changed  into  various 
neutral  molecules,  free  radicals,  and  ions  by  pyrolysis  [86],  oxidation 
[87] ,  and  chemi-ionization  mechanisms  [88]  .    These  primary  species  are 
low  molecular  weight  and  short  chain  or  small  cyclic  molecules  and 
ions.    The  reaction  of  the  various  fragments  with  one  another  and  the 
initial  fuel  and  oxidizer  leads  to  formation  of  monocyclic  and 
polycyclic  aromatic  hydrocarbons,  PCAH's.    This  transformation  is 
referred  to  as  "nucleation"  and  is  one  of  the  least  understood  parts  of 
soot  formation.    Two  schools  of  thought  differ  on  the  major  cause  of 
soot  nucleation.    One  emphasizes  the  role  of  ions  [88],  while  the  other 
highlights  the  importance  of  neutral  free  radicals  [89] . 
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The  entire  process  of  soot  formation  proceeds  very  rapidly,  on  the 
order  of  milliseconds,  even  in  reduced  pressure  flames.  Kinetic 
arguments  favor  the  ion-molecule  reaction  mechanism,  due  in  part  to  the 
very  fast  reaction  rates  possible  for  some  ions.    Another  reason  is  that 
isomeric  rearrangement  of  ions  is  known  to  proceed  very  rapidly  [90]. 
This  is  very  important  in  explaining  the  formation  of  cyclic  aromatic 
canpounds  in  sooting  flames,  since  cyclization  is  expected  to  proceed 
only  slowly  for  neutral  free  radicals.    At  least  some  thermodynamic 
arguments  tend  to  favor  the  ion-molecule  mechanism  for  soot  nucleation 
and  also  cast  doubt  on  the  feasibility  of  a  neutral  reaction  mechanism 
to  lead  to  the  formation  of  very  large  cyclic  molecules  [85]. 

A  better  understanding  of  the  structures  and  reactivities  of  ions 
during  nucleation  of  small  molecules  to  soot  should  help  clarify  the 
mechanism  of  soot  formation.    One  ion  observed  in  sooting  flames,  which 
has  been  suggested  as  an  intermediate  ion  leading  to  soot,  is 
C^H^  [85].    Certainly  more  studies  dealing  with  the  structure  and 
reactivity  of  this  ion  would  help  in  the  understanding  of  soot  formation 
and  ion  chemstry  in  general. 

Thermochemistry  of  Criit  Ions 

The  number  of  atoms  in  the  ion  is  such  that  several  different 

linear  and  cyclic  isomers  are  possible.    These  isomers  cannot  be 
differentiated  directly  by  mass  spectrometry  and  thus  can  only  be 
inferred  from  properties  such  as  measured  heats  of  formation  and 
chemical  reactivities.    Evidence  has  been  reported  for  the  existence  of 
more  than  one  stable  isomer  but  conclusive  data  as  to  which  isomers 
exist  is  still  elusive  [45]. 
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The  appearance  potentials  of  C^H^  ions  from  a  large  number  of 
alkylbenzenes  and  chloro-alkylbenzenes  have  been  measured  [91]-  Results 
of  this  work  indicate  that  in  these  compounds  the  C^H^  generated  has  a 
heat  of  formation  less  than  or  equal  to  284  ±  5  kcal/mole.    An  earlier 
study  measured  the  appearance  potential  of  C^H^  from  HCsCCHsCHCH^ 
and  calculated  a  heat  of  formation  of  277  kcal/mole  for  the  ion  formed 
in  the  reaction 

C^H^   >  C^H^  +  H 

[92].    The  linear  ionic  structure  is  assumed  to  result  from  the 
fragmentation  of  this  ion,  thus  leading  to  a  heat  of  formation  of  about 
280  kcal/mole  for  linear  C^H^  isomers. 

Another  possible  isomer  of  the  C^H^  ion  is  a  five  membered  ring 
(cyclopentadiene  cation).    Using  the  most  recent  ionization  potential  of 
8.41  eV  for  the  neutral  parent  radical  [93]  and  accepting  a  value  of 
60.9  kcal/mole  for  the  heat  of  formation  of  the  neutral  [94]  leads  one 
to  AHj.  =  255  kcal/mole  for  the  cyclic  cation.    A  separate  study  reported 
the  appearance  energy  of  C^H^  from  C^E^Cti2  (3.55  eV)  resulting  in  a 
heat  of  formation  of  242  kcal/mole  [95].    This  lower  value  suggests  the 
existance  of  another,  more  stable,  isomer  for  the  CcHJ  ion. 

Experimental 

All  experiments  were  performed  at  NBS  using  the  single  micro- 
computer controlled  icr  mass  spectrometer  described  in  the  single 
microcomputer  secton  in  Chapter  3-    Norbornadiene  was  used  to  generate 
two  different  isomers  of  C^H^.    The  internal  energy  of  the  ions  was 
controlled  via  charge  transfer  reactions  to  norbornadiene  involving  CO 
(I. P.  =  14.0  eV  [84]).    A  typical  experiment  used  1  X  10"^  torr  of 
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norbornadiene  with  1  to  2  X  10"^  torr  of  a  reactant  gas  with  CO  then 
added  to  a  total  pressure  of  about  1.5  X  10      torr.    This  procedure 
insured  that  most  of  the  hydrocarbon  ions  were  produced  by  charge 
transfer  and  not  by  direct  electron  impact  of  norbornadiene,  and  also 
increased  the  concentration  of  C^H^  in  the  icr  cell  compared  to  what 
might  be  obtained  by  direct  electron  impact.    The  CO*  ions  were  produced 
by  electron  impact  on  CO  of  25  to  55  eV  electrons  gated  for  5  msec  with 
40  namp  of  current  seen  at  the  collecter  plate.    C^H^  ions  also  were 
prepared  by  charge  transfer  from  CO"*"  to  p-chlorotoluene  and  benzyl- 
chloride  under  the  conditions  stated  above. 

Proton  transfer  reactions  of  one  of  the  C^H^  isomers  were  observed 
under  conditions  similar  to  those  mentioned  above,  except  that  the  ions 
were  produced  from  charge  transfer  reactions  with  Ar*  (I.P.(Ar)  = 
15.6  eV).    All  experiments  were  performed  at  a  magnetic  field  strength 
of  13  kG  and  a  gas  temperature  of  61    C  as  measured  by  an  iron- 
constantan  thermocouple  mounted  to  the  icr  cell. 

Ion-molecule  reaction  pathways  were  determined  from  double 
resonance  experiments  where  the  reactant  ion  was  ejected  out  of  the  cell 
while  a  mass  sweep  was  performed  through  the  possible  product  ions.  The 
double  resonance  oscillator  was  gated  on,  then  off,  while  averaging 
signals  for  each  ion  at  a  particular  frequency.    The  averaged  signal 
from  the  oscillator  "off"  experiments  was  then  subtracted  from  the  "on" 
data.    Product  ions  appeared  as  negative  going  peaks  at  their 
characteristic  cyclotron  frequency. 
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Results 

Kinetics  and  Products 

The  two  different  isomers  of  C^H^  from  norbornadiene  are  evident 

in  the  decay  curve  shovm  in  Figure  6.1.    Table  6.1  lists  the  calculated 

rate  constants  and  reaction  products  observed  for  the  "reactive"  C^H^ 

isomer.    The  rate  constants  reported  were  determined  by  subtracting  the 

unreactive  ion  signal  from  the  decay  curve  then  running  a  linear  least 

squares  anaylsis  on  a  semi-log  plot  of  the  reactive  part  of  the  data. 

Due  to  the  fast  reaction  of  C^H^  with  norbornadiene,  rate  constants 
—  1 0 

less  than  1X10       cm-'/sec  with  the  added  reaction  gas  could  not  be 
determined  with  any  useful  precision. 

Proton  Affinities 

Four  proton  transfer  reactions  of  the  "unreactive"  C^H^  with 
different  bases  were  performed  to  bracket  the  proton  affinity  of  the 
C^Hj^  neutral  produced  in  the  reaction 

C5H5  +  BASE  — >  C^Hij  +  BASEH* 

The  results  of  the  proton  transfer  reactions  are  given  in  Table  5.2  and 
indicate  a  proton  affinity  of  227.9  ±  0.3  kcal/mole  for  the  C^Hj^  neutral 
[96]. 

Discussion 

Reaction  Products  and  Kinetics 

Isotope  studies  on  carbon-13  labeled  norbornadiene  indicate  the 
main  route  for  C^H*  formation  is  the  follov;ing  [97]. 

C^hJ   >  C3HJ  +  C2H2   >  +  H  +  C2H2 


Two  C^H^  structures  are  produced  from  norbornadiene  [45]  but  only  one, 
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Table  6.1.    Experimental  rate  constants  and  reaction  products  for  the 
"reactive"  C^H^  isomer  reacting  with  several  neutral 
molecules. 


Reaction  k  (cm^/sec)  /  10 


s4  + 

^2°2 

0.6 

 > 

<  0.1 

+ 

HCHCCECH 

2.6 

+ 

0.7^ 

+ 

H2C=C=CHCH2 

s<  ^ 

«2 

2.5 

+ 

1.0 

+ 

H2C=CHCH=CH2 

1.6 

+ 

h\ 

<  0.1 

J. 

Products 

12 . 1 

+ 

1.2 

(norbornadiene) 

^io< 
So«ii  + 

b 

S«2 

+ 

1.2 

+ 

0.2 

+ 

C,H_CH=CH^ 

Products 

0.9 

b 

+ 

0.3 

^13^13 

^13^12  ^ 

H 

+ 

^^18 

(isooctane) 

<  0.1 

+ 

5.0 

+ 

1.7 

+ 

C^H3CH2C1 

S»7  + 

C^H^Cl 

3.6 

+ 

1.3 

+ 

p-ClCgH^CH^ 

 > 

C3H3CI 

4.3 

+ 

0.6 

a)  95%  confidence  intervals  given  when  more  than  one  determination  of  the 
rate  constant  was  made. 

b)  Major  reaction  pathway. 
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Table  6.2.     Results  from  observations  of  C^H^  (unreactive  structure) 
formed  in  norbornadiene  with  a  few  additive  molecules. 


Additive 
Molecule 

Proton 
Affinity 
(kcal/mole) 

Ionization^ 
Potential 

(eV) 

Mode  of 
Reaction 

k 

exp 

(cm"^/sec) 

225.9 

7.82 

No 

Reaction 

<io-ii 

Azulene 

227.6 

7.4 

No 

Reaction 

<io-i^ 

(n-C2H^)2NH 

228.3 

7.84 

Proton 
Transfer 

=10-1° 

233.0 

7.5 

Proton 
Transfer 

>io-9 

a)    Proton  affinities  taken  from  Aue  and  Bowers  [98]. 


b)     Ionization  potentials  taken  from  Rosenstock  et  al.  [84]. 
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the  "reactive",  is  seen  when  aromatic  six-membered  rings  are  used  as 
precursor  molecules  [97] .    In  this  case  the  C^H^  is  formed  from  the 
benzyl  or  tropylium  CjHj  structure  by  the  reaction 

Our  results  agree  with  the  investigation  mentioned  above  in  that  only 
the  "reactive"  isomer  was  observed  when  p-chlorotoluene  and  benzyl- 
chlorotoluene  were  used  as  the  source  of  C^H^  (see  results  at  the 
bottom  of  Table  6.1). 

The  reaction  products  shown  in  Table  5.1  indicate  most  reactions 
proceed  by  condensation  followed  by  fragmentation.    Many  of  the  product 
ions  have  been  seen  in  flames  and  have  a  low  hydrogen  to  carbon  atom 
ratio  as  found  in  PACH's.    Reaction  of  C^H^  with  the  larger  neutrals 
indicates  that  ions  of  high  molecular  weight  can  be  formed  quickly  by 
ion-molecule  reactions.    Benzene  and  acetylene  do  not  react  at  an 
appreciable  rate  with  the  C^H^  ion,  yet  diacetylene  and  aromatics  with 
side  chains  do.    This  suggests  soot  nucleation  may  proceed  by  adding  a 
few  large  molecules  rather  than  through  addition  of  many  smaller  ones. 

As  can  be  seen  from  Table  6.1,  the  rates  of  many  of  the  reactions 
reported  in  this  study  appear  to  be  occuring  at  a  fraction  of  the 
collision  rate  (which  is  about  1  X  10~9  cm^/sec  for  most  reactions 
listed).    Little  previous  work  dealing  with  this  ion  is  reported  in  the 
literature.    The  only  literature  value  known  to  the  author  for  a 
reaction  similar  to  those  examined  in  this  study  is 

C5H;  +  C6H5CH3   >  C^^h5  +  CH4 

with  a  rate  constant  of  7  X  lO"''^  cm^/sec  [76].    This  rate  compares 
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reasonably  well  with  our  value  of  1.2  ±  0.2  X  10       cm-^/sec  determined 
for  partially  deuterated  toluene  reacting  with  C^H^  from  norbornadiene . 

The  initial  steps  of  soot  formation  have  been  computer  simulated 
and  demonstrate  that  large  ions  evolve  rapidly  enough  to  account  for  the 
appearance  of  soot  even  when  using  a  modest  reaction  rate  constant  of 
2  X  10"^^  cm^/sec  for  all  reactions  [88].    Our  rates  show  reactions  of 
C^H^  an  order  of  magnitude  faster  than  those  used  in  the  simulation 
while  producing  various  ions  found  in  flames. 

CcHg  Thermochemistry 

The  proton  affinity  of  C^H^  determined  in  our  experiments 
corresponds  to  the  energy  released  in  the  reaction 

C5HH  +  H+   >  C5H+ 

where  C^H^  is  the  unreactive  species.    Using  our  value  for  the  proton 
affinity  of  C^Hi^  and  a  heat  of  formation  of  365.7  kcal/mole  for  H"*"  [84] 
leads  to  a  difference  between  the  heats  of  formation  of  "unreactive" 
C^H^  and  of 

Ati^iC^ap  -  AHf.(C5Hi,)  =  137.8  kcal/mole. 

The  cannot  be  the  cyclopentadiene  cation  because  the  ionization 

potentials  are  such  that  this  species  would  transfer  charge  to  all  the 
bases  in  Table  6.1. 

Assuming  that  the  "unreactive"  C^H^  is  not  the  higher  energy 
linear  ion,  and  accepting  that  it  is  not  the  cyclopentadiene  ion,  then 
the  most  likely  structure  is  that  below  [99]. 

|.+^C-CH=CH, 
HC^  ^ 
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Since  McCrery  and  Preiser  formed  an  ion  with  a  heat  of  formation  of  242 
kcal/mole  [95],  and  the  next  most  stable  structure  (cyclopentadiene 
cation)  has  an  estimated  heat  of  formation  equal  to  255  kcal/mole,  we 
will  assume  that  they  were  forming  this  ion.    If  this  ion  transferred  a 
proton,  then  the  only  molecule  that  could  be  formed  without  rearrangment 
is  the  following. 

HCv 
ll>=C=CH 

The  estimated  heat  of  formation  of  this  molecule  is  about  130  kcal/mole 
(calculated  below)  leading  to 

[Mi^iC^np  -  AHf.(C5H^)]  =  112  kcal/mole  . 

If  rearrangement  occurs  during  proton  transfer,  a  square  planar 
structure  may  result;  however,  this  is  not  expected  due  to  the  amount  of 
strain  present  in  the  molecule.    Other  C^Hi^  molecules  which  might  result 
if  rearrangement  occurs  during  proton  transfer  are 

CH2=C=CHC=CH 
[AHfCC^Hp  -  AHf(C5Hj^)]  =  131  kcal/mole 

or 

HC=CCeCCH3 
[AHf-CC^H*)  -  AHf-CC^Hij)]  =  140  kcal/mole 

or 

H2C=C=C=C=CH2 
[AHf.(C5H^)  -  AHfCCgHij)]  =  150  kcal/mole 

None  of  the  heats  of  formation  for  these  neutrals  are  known  and  all  have 
been  estimated  by  using  the  technique  of  "macro-incrementation"  [100]. 


II   "C=C=CH      =      |i^C=CH        +      CH-=C=CH-      -      CH  =CH, 
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Instead  of  adding  up  the  contributions  of  the  bonds  as  in  Benson 
estimations  [101],  one  estimates  the  heat  of  formation  using  "groups" 
from  molecules  of  known  AH^,  choosing  molecules  having  the  same 
structural  features  (or  strain).    For  example, 

HC^  HC, 

II    C=C=CH      =  II 
HC^  HC 

130  (est.)  =    96.4  [102]    +      45.6  [103]    -      12.45  [103] 
Using  this  method,  an  estimate  of  the  heat  of  formation  of  the  radical 
corresponding  to  the  "unreactive"  ion  is  123  kcal/mole.  There- 

fore, we  can  estimate  an  ionization  potential  of  about  5.2  eV  for  this 
radical,  and  the  observation  that  the  ion  does  not  transfer  charge  to 
the  amines  with  ionization  potentials  of  lA  to  7,8  eV  is  consistent 
with  our  proton  affinity  data. 

Conclusion 

The  C^H^  ion  does  react  rapidly  with  some  of  the  species  found  in 
flames  to  produce  much  larger  ions  which  are  also  seen.  The  reactions 
are  fast  enough  to  account  for  the  rapid  generation  of  soot  in  fuel  rich 
flames.  These  ions,  however,  are  near  ground  state  in  energy  unlike  a 
significant  number  of  species  found  in  flames.  These  excited  ions  may 
be  expected  to  be  even  more  reactive  due  their  added  internal  energy. 

The  "unreactive"  C5H5  ion  is  most  likely  the  vinylcyclopropenyl 
structure  shown  earlier.    On  the  basis  of  thermodynamic  arguments,  it 
can  be  inferred  that  this  ion  deprotonates  to  give  a  linear  C^Hjj 
neutral.    The  "reactive"  C^H*  could  be  the  cyclopentadiene  ion  or  have 
a  linear  structure.    However,  more  experiments  are  necessary  to 
determine  its  structure  more  conclusively. 


CHAPTER  7 
ION  CHEMISTRY  OF  CHoI  AND  CF,I 


Introduction 

Ion  chemistry 

The  radiation  chemistry  of  gas  phase  CH^I  [104,105]  ,  CF^I 
[106,107,108],  and  mixtures  of  the  two  [109]  have  been  investigated  in 
recent  years.    In  the  CF^I  system,  the  initial  absorption  of  radiation  by 
CF^I  generates  the  ionic  species  CF^,  I*,  and  CF^I*  [105].    These  ions 
are  believed  to  produce  I2  via  the  I  radical  formed  by  charge  transfer 
of  I*  to  CF^I,  or  ion-molecule  reactions  which  generate  I2  or  I2 
[108].    Free  radical  scavenger  studies  support  the  role  of  ions  in  the 
formation  of  I2  from  the  radiolysis  of  CF^I  [107].    Iodine  seen  in  the 
radiolysis  of  CH^I  [104]  and  mixtures  of  CH^I  and  CF^I  [109],  may  be 
produced  through  a  similar  ionic  mechanism.    End  product  analysis  of 
irradiated  CF-^I  and  CH^I  mixtures  shows  a  variety  of  cross  product  species, 
CH^F  and  CF^H  for  examples,  and  there  has  been  speculation  that 
ion-molecule  reactions  might  be  involved  in  forming  some  of  them  [109]. 

The  ion  chemistry  of  CH^I  has  been  reported  previously 
[110,111],  as  well  as  deuterated  studies  which  probed  the  mechanism  of 
the  reactions  [112].    The  ion  chemistry  of  CF^I  has  also  been  published 
but  with  some  conflicting  results  [113,114].    Both  the  CH^I  and  CF^I 
work  show  the  production  of  I"^  as  a  fragment  from  ionization,  while  one 
of  the  CF^I  studies  reports  the  formation  of  I2  [113].    The  ion 
chemistry  of  the  mixture  has  not  been  reported  to  date,  yet  such  a  study 

would  shed  light  on  the  role  of  ions  in  the  radiolysis  of  this  system. 
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Detector  limitations 

Due  to  inherent  characteristics  of  the  icr  a.o.  detector  used  in 
these  studies  (see  Chapter  2  for  details),  high  mass  signals  are  weak 
compared  to  low  mass  signals  under  similar  conditions.    This  makes  the 
detection  of  ions  containing  more  than  one  iodine  difficult  at  best. 
The  detector  sensitivity  is  inversely  proportional  to  ion  mass;  thus  the 
S/N  ratio  at  m/e  265  (CF^ICFp  is  about  four  times  less  than  that  at 
m/e  69  (CFp  as  determined  by  use  of  Equation  2.6.    Signal  averaging  is 
necessary  to  pull  the  signal  out  of  the  noise  yet  one  must  average  16  of 
the  high  mass  experiments  just  to  reach  the  same  S/N  ratio  as  one  at  the 
low  mass.    The  S/N  problem  is  compounded  when  low  intensity  ions  are 
observed  in  the  high  mass  range. 

While  the  ions  are  losing  intensity  at  high  masses  one  must  also 
contend  with  a  loss  of  resolution.    A  zero  pressure  resolution  for  m/e 
69  at  12  kG  magnetic  field  strength  is  1500  (calculated  from  Equation 
2.3),  v*iereas    a  maximum  of  only  392  is  possible  for  m/e  265.  Resolu- 
tion is  not  a  problem  in  the  icr  study  of  CF^I  since  the  possible 
products  are  widely  separated  in  the  mass  spectrum.    However  for  the 
CH^I  system,  the  gain  or  loss  of  one  hydrogen  may  be  difficult  to  see. 
High  pressures  collisions  broaden  the  peaks,  further  degrading  the 
resolution  by  spreading  the  peak  out  and  decreasing  its  maximum 
intensity. 

Experimental 

Chemicals  were  purified  and  supplied  by  Mr.  Dawit  Teclemariam  and 
degassed  by  repeated  freeze-pump-thaw  cycles.    Ionization  gauge 
sensitivities  for  CH^I ,  CF3I ,  and  a  1:1  mixture  were  determined  in  the 
same  fashion  as  described  in  the  instrument  calibration  section  of 
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Chapter  5.    Values  of  1.51  ±  0.09,  2.19  ±  0.17,  and  1.82  ±  0.19 
respectively  were  obtained  (95%  confidence  intervals  given).  These 
sensitivities  are  such  that  the  1:1  mixture  value  equals,  within 
experimental  error,  half  the  sum  of  the  component  values.  All 
experiments  were  performed  on  the  dual  microcomputer  based  data  system 
previously  described  in  Chapters  3  and  4.     Marginal  oscillator  signals 
were  averaged  from  7  to  30  times  for  each  experiment.    Double  resonance 
experiments  were  performed  by  sweeping  across  the  frequency  of  the  ion 
of  interest  with  a  function  generator  and  averaging  the  output  of  the 
m.o.  detector  with  the  computer.    Chemically  coupled  ions  are  indicated 
by  a  negative-going  peak  in  the  spectrum  [115]. 

All  fragment  ions  were  formed  by  70  eV  electron  impact  on  reactant 
gases  in  the  3  X  10"*^  to  5  X  10"^  torr  pressure  range.    In  a  typical 
experiment  the  electron  beam  was  pulsed  for  5  msec  while  measuring 
currents  of  30  to  600  namp  at  the  collector  plate.    The  trapping,  side, 
and  end  plates  of  the  icr  cell  were  at  +1 ,  -1 ,  and  -1  v  respectively. 

Rggults 

Pure  CH-^I 

The  parent  ion,  CH^I*,  and  two  fragment  ions,  CH^  and  I*,  are 

produced  by  electron  impact  of  the  CH^I  neutral.    The  product  ions 

CH^ICHI  and  CH^I^  were  both  seen  but  the  condensation  dimer  (CH2l)2 

could  not  be  detected  over  the  pressure  range  used  in  these 

experiments.    Ion  decay  signals  were  analyzed  by  linear  least  squares 

analysis  of  semi-log  plots.    The  results  of  the  experiments  on  CH3I  are 

given  in  Table  7.1.     Figure  7.1  shows  the  presence  of  CH2I*  as  a  shoulder 

on  the  CH^I"^  mass  peak;  however,  due  to  the  limited  resolution  of  the 

instrument  at  this  mass,  quantitative  results  could  not  be  obtained  for 
the  ion . 
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Pure  CF-,1 

Ionization  by  electron  impact  of  trifluoromethyl  iodide  generated 
the  parent  ion  and  two  fragments:  I"*"  and  CF^.    These  species  reacted  to 
form  various  products:  CF2I*,         CF^ICF^,  CF2l2»  and  (CF^Dl  in 
order  of  decreasing  intensity.    The  latter  was  seen  at  very  high 
pressures  (5  X  10"^  torr)  and  high  collector  currents  (600  namp)  as 
a  very  broad  and  weak  peak.    The  results  of  the  experiments  are  given  in 
Table  7.2. 

CH-.I  ■>•  CF-^I  Mixture 

The  ion  chemistry  of  the  mixture  looks  very  much  like  the  ion 
chemistry  of  CH^I.    Fragment  ions  from  electron  impact  of  both  CH^I  and 
CF^I  are  seen  in  the  mass  spectrum  of  the  mixture.    The  product  ions 
CFgl"*",  CH^ICH^,         and  CH^ll    appear  in  the  mass  spectrum  of  the 
mixture.    The  results  of  reactions  characteristic  of  the  mixture  are 
given  in  Table  7.3. 

Discussion 

Pure  CH->T 

Only  a  limited  number  of  rate  constants  have  been  reported  in  the 
literature  for  the  CH^I  system.    One  reaction  studied  previously  is 

CH^I*  +  CH^I  — >  CH^ICH*  +  I 

resulting  in  measured  rate  contsants  of  0.12  X  10"^*^  [119]  and  0.06 
X  10"      cm^/sec  [110],    A  more  recent  value  of  0.237  X  10"^°  [120] 
compares  well  with  0.24  X  10"''°  cm^/sec  from  this  study.    The  previous 
studies  were  done  by  forming  the  CH^I"^  ion  with  10  to  11  eV  photons 
[119]  or  electrons  [110]  thus  generating  an  ion  with  little  excess 
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internal  energy;  yet  their  results  are  in  good  agreement  with  those  from 
our  experiments,  in  which  the  ions  were  formed  by  impact  with  70  eV 
electrons.    A  reaction  of  CH3lCH^  with  CH3I  was  observed  to  occur 
slowly  with  a  rate  constant  of  9  ±  2  X  lO'''^  cm^/sec.  Unfortunately, 
the  products  of  this  reaction  were  not  determined. 

The  CH3I2  detected  in  this  study  has  been  reported  in  only  one 
earlier  investigation  [ill].    This  is  due  in  part  to  the  limited  mass 
range  of  the  icr  of  Beauchamp  et  al.  [110].    The  photoionization 
experiments  concentrated  on  the  formation  of  the  association  product  and 
may  have  missed  this  low  intensity  fragment  from  the  association  process 
[119].    Double  resonance  experiments  verified  the  production  of  CH^ICH^ 
and  CH3II  from  the  parent  ion. 

The  1+  fragment  ion  reacts  (k  =  9.7±2.8X  10" ''Q  cm^/sec)  via 
charge  transfer  with  the  neutral  parent  molecule  at  the  collision  rate. 
This  is  expected  due  to  the  difference  in  ionization  potentials  of  the 
two  species;  I  =  10.451  eV  and  CH3I  =  9.54  eV  [84],    This  also  explains 
the  absence  of  I*,  seen  in  the  CF3I  system,  since  it  would  most  likely 
be  produced  by  the  reaction 

I-'  +  CK3I  >  1+  +  CH3  . 

This  reaction  channel  evidently  cannot  compete  with  the  fast  charge 
transfer  reaction. 

The  CH|  fragment  ion  is  believed  to  react  rapidly  via  the  reaction 
CH^  +  CH3I  >  CH4  +  CH2I-' 

as  suggested  earlier  [110].    However,  the  reported  ionization  potential 
of  9.83  eV  for  the  methyl  radical  [84]  indicates  charge  transfer  is  also 
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possible  between  the  fragment  ion  and  neutral  parent  molecule.    The  rate 
constant  obtained  (19.2  ±  1.2  X  10"     cm-'/sec)  has  no  literature  value 
for  comparison.    However,  the  reaction  appears  to  proceed  at  the 
collision  rate.    The  low  intensity  of  CH2I"*'  relative  to  CH-^I'*'  seen  in 
Figure  7.1  suggests  that  the  charge  transfer  reaction  is  dominant.  This 
is  in  contrast  to  isotope  studies  on  reactions  of  CD^  with  CH^I  which 
indicate  both  reactions  do  occur  but  the  reaction  producing  methane  is 
of  primary  importance  [112].    This  reaction,  which  could  explain  the 
formation  of  methane  in  the  radiolysis  of  CH^I ,  has  been  postulated  by 
Teclemariam  and  Hanrahan  [109]. 

Pure  CF-,1 

The  parent  ion  produced  from  CF^I  reacts  to  form  a  variety  of 
products.    Analogous  to  the  CH^I  system,  CF^ICF^  and  CF2I2  are 
produced  through  the  reactions 


and 


CF3I+  +  CF3I  ™>  CF^ICF^  +  I 


 >  CF3I2  +  CF3 


These  reactions  occur  at  the  same  or  slightly  faster  rate  than  in  the 
CH^I  system.    Earlier  work  did  not  report  the  formation  of  CF^ICF^  yet 
did  report  ions  of  heavier  mass  {CF-^l2  and  CF2I2CFP  [113].    A  low 
pressure  study  (1.7  X  10""^  torr  maximum)  on  CF^I  failed  to  identify  any 
ions  of  mass  greater  than  CF^I*  [114].    This  is  to  be  expected  due  to 
the  small  rate  constants  involved  for  many  of  the  reactions  forming  the 
heavy  product  ions.    The  formation  of  the  (CF^D^  dimer  is  observed  at 
5  X  10"^  torr,  unlike  the  corresponding  dimer  ion  in  pure  CH^I.  The 
association  reactions  in  methyl  iodide  show  an  increase  in  rate  when 
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deuterium  is  substituted  for  the  hydrogens  in  the  methyl  iodide  molecule 
[120].    This  isotope  effect  is  attributed  to  a  decrease  in  the 
vibrational  frequency  of  the  C-H  bond  and  the  resulting  increase  in  the 
density  of  vibrational  states.    When  fluorine  atoms  take  the  place  of 
hydrogens  in  methyl  iodide  one  would  expect  an  enormous  change  in 
vibrational  frequency  along  with  a  greatly  increased  density  of  states. 
This  makes  the  fluorinated  compound  better  suited  to  partition  the 
energy  of  association  internally,  thus  stabilizing  the  ion  without  the 
help  of  a  third  body. 

The  fragment  ion  I*  reacts  to  form  I2  as  demonstrated  by  icr 
double  resonance  experiments.    The  charge  transfer  reaction: 

I*  +  CF3I  ™>  CF3I+  +  I 

was  reported  earlier  [113]. 

Analogous  to  the  methyl  iodide  ion  chemistry  the  fragment  CF^ 
reacts  to  form  CF2I*  via  the  reaction: 

CF^  +  CF3I   >  CFn  +  CF2I*  . 

Previous  work  in  a  time  of  flight  mass  spectrometer  with  a  high  pressure 
source  did  not  detect  a  decay  mechanism  for  the  CF2I*  ion  [113].  The 
difference  between  this  observation  and  the  above  reaction  seen  in  an 
icr,  was  explained  in  terms  of  energy  differences  for  the  reacting 
ions.    It  was  suggested  that  the  ions  in  the  icr  retain  some  excess 
energy,  while  those  in  the  high  pressure  source  mass  spectrometer  were 
presumably  therraalized.    Our  experiments  do  show  a  slow  reaction  with  a 
rate  constant  of  1,4  X  10"''''  cm^/sec.    However,  any  excess  internal 
energy  in  the  ion  should  be  removed  through  the  numerous  non-reactive 
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collisions  suffered  before  reaction  occurs.    The  total  loss  of  CF^ 
cannot  be  accounted  for  by  the  above  reaction  alone,  yet  the  charge 
transfer  reaction 

CF3  +  CF3I  — >  CF^I*  +  CF3 

proposed  in  an  earlier  study  [113]  is  too  endothermic  to  be  seen  in  an 
icr  experiment.    This  is  due  to  the  difference  in  ionization  potentials 
of  CF3I  (10.23  eV)  and  CF^  (9.25  eV)  [84]. 

CH-jl  +  CF-,1  Mixture 

The  CH^I  ion  chemistry  appears  to  dominate  in  the  ion  chemistry 
of  the  mixture.    This  is  due  in  part  to  the  loss  of  CF^I*  through  charge 
transfer  in  the  reaction 

CF^I*  +  CH3I  -— >  CF3I  + 

which  is  thermodynamically  favored  due  to  the  difference  in  ionization 
potentials  between  CH^I  (9.54  eV)  and  CF^I  (10.23  eV) .    The  slower 
reactions  of  CF^I*  which  produce  the  species  CF^ll,  CF^ICF^  and 
(CF^Dl  cannot  compete  with  this  reaction  channel.    The  above  charge 
transfer  reaction  was  verified  by  icr  double  resonance  experiments. 

The  CH^  ion  generated  from  electron  impact  of  methyl  iodide  does 
not  react  with  the  neutral  CF^I  molecule  since  its  decay  rate  in  the  \ 
mixture  can  be  accounted  for  by  the  reactions  seen  in  the  pure  CH^I 
system.    The  CF^  ion,  however,  does  react  with  the  CH^I  molecule  with  a 
rate  constant  of  3.0  ±  1.2  X  10"''°  cm^/sec.    The  equation  for  this 
reaction  is  assumed  to  be 

CF^  +  CH3I  --->  CF3H  +  CH2l'^  (7.1). 
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Unfortunately,  no  previous  work  has  been  reported  for  comparison  of  rate 
constants  for  the  mixture. 

The  iodine  ions  produced  in  the  mixture  appear,  not  surprisingly, 
to  undergo  the  same  reactions  observed  in  the  two  pure  systems.  The 
production  of  l|  can  be  explained  by 

I*  +  CF3I  — >  I2  +  CF^ 

or  by 

CHgl*  +  CF3I   >  CH3CF2  +  I2  (7.2). 

There  seems  to  be  an  absence  of  cross  product  ions  containing  hydrogens 
and  fluorines  but  neutral  cross  product  species  do  seem  probable. 

CpHQlugion 

The  ion  chemistry  of  pure  CH^I  and  pure  CF^I  are  similar  in  many 
respects  including  the  production  of  CHi^,  CH2I''",  and  their  fluorinated 
analogs.    The  ion  chemistry  of  the  mixture  resembles  the  chemistry  of 
the  CH^I  system  due  to  charge  transfer  from  CF^I*  to  CH^I •    The  mixture 
does  not  show  the  presence  of  ionic  cross  product  species,  but  the  ion 
molecule  reactions  determined  in  this  study  do  indicate  the  formation  of 
neutral  cross  product  species  (see  reactions  7.1  and  7.2  above). 

A  marginal  oscillator  detector  on  an  icr  mass  spectrometer  is  not 
normally  noted  for  its  operation  in  the  high  mass  region.    However,  with 
the  signal  averaging  capabilities  of  a  computer  controlled  data 
acquisition  system  low  intensity  ions  up  to  mass  400  can  be  probed  for 
chemical  information. 


APPENDIX  A 
BASIC  PROGRAM 


The  subsequent  text  and  program  listing  describe  the  BASIC  program 
used  to  control  the  marginal  oscillator  detector  on  our  icr  spectro- 
meter.   Understanding  of  the  code  requires  a  good  knowledge  of  the 
APPLESOFT  II  BASIC  language  and  the  DOS  3-3  disk  operating  system.  The 
appropriate  manuals  should  be  consulted  for  the  function  and  syntax  of 
the  various  commands  [121]. 

Due  to  the  interpreting  mechanism  of  the  BASIC  language ,  program 
execution  time  is  shortened  by  placing  more  frequently  used  routines 
early  in  the  program  and  less  used  code  towards  the  end.    Memory  space 
is  conserved  and  execution  time  shortened  by  not  using  remark  statements 
in  the  program.    To  aid  in  deciphering  of  the  code,  Table  A.I  identifies 
numbered  sections  of  the  BASIC  routine  and  their  purpose. 

The  most  important  part  of  software  initialization  is  the  defining 
of  certain  variables  and  an  array  used  by  the  machine  code  programs 
described  in  Appendix  B.    The  first  four  variables  defined  must  be  the 
data,  1)%,  number  of  passes,  KIM  program  vector,  KV?,  and  KIM 

status,  KE^.    The  names  are  not  important  but  they  must  be  integers  in 
the  order  given.    Next,  the  first  array  defined  must  be  the  pulse 
parameter  integer  array  PP^(8,1),  nine  long  by  two  wide.    All  variables 
must  be  as  indicated  to  insure  proper  communication  between  the 
routines . 
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m 

Table  A.I.    List  of  ICR-MO  code  sections  along  with  line 
numbers  and  brief  description  of  the  routine. 


Line 

Number 

From 

la 

Routine  Description 

10 

30 

Pulse  parameter  set  up 

50 

54 

Data  plotting 

100 

1 10 

Get  data  and  key  command  check 

200 

250 

Continuous  scan  experiment 

260 

370 

Detect  delay  sweep  experiments 

400 

420 

Signal  average  scan  experiment 

800 

840 

String  input  and  evaluation 

900 

910 

Key  command  input 

920 

1080 

Print  experiment  parameters 

1300 

1320 

Print  experiment  numerical  values 

1500 

1680 

Key  command  interpreter 

1700 

1920 

Experiment  parameter  alteration 

2000 

2170 

Print  experiment  conditions 

2200 

2230 

Print  condition  numerical  values 

2405 

2490 

Condition  change  command  interpreter 

2500 

2680 

ExDeriment  condition  alteration 

4000 

4180 

Experiment  parameter  limit  check 

5000 

5070 

Store  data  on  disk 

5100 

5140 

Print  hardcopy  of  parameters  and  conditions 

5400 

5500 

Load  parameter  file  from  disk 

5505 

5550 

Save  parameter  file  on  disk 

5800 

5830 

Disk  error  handler 

5900 

5910 

Disk  catalog 

6000 

6200 

Mass  calculation 

6500 

6597 

Rate  constant  calculation 

7000 

7040 

Print  data  hardcopy 

7400 

7430 

Load  KIM'S  control  program 

7500 

7900 

Cursor  movement  routine 

8000 

9000 

Numbered  command  display  and  interpreter 

10000 

10040 

Initialize  arrays,  variables,  and  parameter 
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The  plot  data  routine  takes  X,Y  pairs,  checks  for  Y  value  under- 
and  overflows,  then  displays  them  as  points  on  the  high  resolution 
graphics  screen.    The  routine  which  activates  the  Apple  machine  code 
also  checks  for  key  commands  pressed  during  the  machine  code  execution. 
When  a  key  is  pressed  the  machine  code  routine  is  left  and  the  key 
command  interpreted  and  executed.    The  program  then  reenters  the  machine 
code,  waiting  for  data  from  the  KIM  or  another  key.    The  string  input 
and  evaluation  routine  accepts  characters  during  experiment  parameter  or 
condition  changes  and  returns  the  numerical  Vcilue  of  the  expression  in 
the  variable  PAR.     If  the  string  cannot  be  evaluated  the  error  flag,  ER, 
is  set  to  one,  otherwise  it  is  zero.    Disk  errors  handled  by  the 
designated  routine  result  in  the  display  of  the  error  number  and 
sounding  of  the  bell,  then  program  control  returns  to  the  numbered 
command  display  and  input  routine.    The  code  for  printing  hardcopy  of 
the  parameters  and  conditions  assumes  that  the  printer  interface  card  is 
in  slot  1 . 

The  following  is  a  listing  of  the  BASIC  program  ICR-MO. 


1    GOTO  10000 

10  PP$(2,0)  =  TD  »  TW:PP55(2,1)  =  (TD  +  DW)  »  TW:PPJ(1,0)  = 

PP$(2,0):PPJ(1  ,1)  =  PP$(1,0)  +  1:PP5S(3,0)  =  PP*(2, 1 )  :PP$(3, 1 ) 

=  PP?(2,1)  +  1 :  RETURN 
20  PP%m,0)  =  WD  »  TW:PP?(4,1)  =  (WD  +  WW)  »  TW:  RETURN 
30  PP?(7,0)  =  QD  »  TW:PP?(7,1)  =  (QD  +  QW)  »  TW:PP$(8,0)  =  DC 

*  TW:PP?(0,1)  =  GW  »  TW:  RETURN 
50  YV  =  159  -    INT  (YV  *  FE  +  .5):  IF  YV  >  159  THEN  YV  =  159: 

GOTO  54 
52    IF  YV  <  0  THEN  YV  =  0 
54    HPLOT  XV, YV:  RETURN 

100    CALL  2048:  IF  KE%  <    >  0  THEN  A  =    PEEK  (  -  16384)  -  128: 

POKE    -  16368,0:  GOSUB  1500:  GOTO  100 
110  RETURN 

200  TD  =  DD:  GOSUB  10:  HGR  :KV$  =    -  23912 
210    CALL  62450 :ME  =  0 :  FOR  I  =  0  TO  279 

240    GOSUB  100:DA(I,0)  =  D$:N  =  I  +  1 :XV  =  I:YV  =  D?:  GOSUB  50 
250    NEXT  :  GOTO  210 

260    POKE    -  16304,0:N  =  0:XS  =  279  /  DD:  IF  DD  /  DS  >  280 
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THEN  EX  =  0:  PRINT  BE$:  GOTO  900 
270  KV$  =    -  23854:TD  =  0:  GOSUB  10:XV  =  0:PPi(0,1)  =  0:Z  =  1: 

GOSUB  100:  FOR  TD  =  0  TO  DD  STEP  DS:  GOSUB  10:PP$(0,1)  =  GW 

TW:  GOSUB  100:DV  =  DJ:DA(N,Z)  =  DV  /  PA$ 
280  YV  =  DA(N,Z):  GOSUB  50 :Z  =  0:XV  =    INT  (XS  »  TD  +  .5):  IF 

MDJ  =  4  THEN  P?%{0,^)  =  0:  GOSUB  100:DV  =  D$:DA(N,Z)  =  DV  / 

PA$:YV  =  DA(N,Z):  GOSUB  50:Z  =  1: 
290  IF  TD  =  0  AND  Z  =  0  THEN  NEXT 
295  N  =  N  +  1 :  NEXT 

300  KV?  =    -  23912:PPil(0,1)  =  GW  »  TW:  GOSUB  lOOrDV  r 

Di:DA(N,Z)  =  DV  /  PA$:YV  =  DA(N,Z):  GOSUB  50 
305    IF  MD$  =  3  THEN  N  =  N  +  1 :  GOTO  370 

310  C  =  0:T  =  0:U  =  0:V  =  0 :W  =  0 :  FOR  I  =  1  TO  N:X  =  I  -  1 :Y 
=  DA(I,1):C  =  C  +  X:T  =  T  +  Y:U  =  U  +  X  »  X:V  =  V  +  Y  «  Y:W 
W  +  X  *  Y:  NEXT 

320  B  =  (N  »  W  -  C  »  T)  /  (N  »  U  -  C  »  C):AA  =  (T  -  B  »  C)  /  N 
325  XV  =  0:YV  =  AA:  GOSUB  50:YL  =  YV:XV  =  279:YV  =  AA  +  B  »  (N 

-  1):  GOSUB  50:  HPLOT    TO  O.YL 
360    FOR  I  =  0  TO  N  -  1:DA(I,0)  =  DA(I,0)  -  (AA  +  B  »  I):  NEXT 
370  EX  =  0:  GOTO  900 

400    POKE    -  16304,0:KV$  =    -  23854 :TD  =  DD:  GOSUB  10:PP$(0,1) 

=  GW  »  TW:  GOSUB  100 
410    FOR  I  =  0  TO  278:  GOSUB  100:DV  =  DJ:DA(I,0)  =  DV  /  PA$:YV 

=  DA (1,0): XV  =  I:  GOSUB  50:N  =  I  +  1:  NEXT 
420  KV%  =    -  23912:  GOSUB  100:DV  =  D$:DA(N,0)  =  DV  /  PA3;:YV  = 

DA(N,0):XV  =  N:  GOSUB  50 :N  =  N  +  1:EX  =0:  GOTO  900 
800  ER  =  0:B$  =    CHR$  (32) 

810    GET  A$:  IF  A$  <    >    CHR$  (13)  THEN  B$  =  B$  +  A$:  PRINT 

A$; :  GOTO  810 
820    IF    LEN  (B$)   <  2  THEN  ER  =  1 :  RETURN 
830  A$  =    MID$  (B$,2,1):  IF  A$  =  "0"  THEN  PAR  =  0:  RETURN 
840  PAR  =    VAL  (B$):  RETURN 
900    GOSUB  920 

910    GOSUB  1680:  VTAB  PV*:  GET  A$:A  =    ASC  (A$):  GOSUB  1500: 

GOTO  910 
920    TEXT  :  HOME 

930    PRINT  "MO  EXPERIMENT  PARAMETERS":  PRINT 
940    PRINT  "GRID  PULSE  WIDTH";  TAB(  27);GW;  TAB(  36); "MSEC" 
950    PRINT  "DETECT  PULSE  DELAY";  TAB(  27);DD;  TAB(  36); "MSEC" 
960    PRINT    TAB(  14);"STEP";  TAB(  27);DS;  TAB(  36);"MSEC" 
970    PRINT    TAB(  14);"WIDTH";  TAB(  27);DW;  TAB(  36);"MSEC" 
980    PRINT  "W2  PULSE  DELAY";  TAB(  27);WD;  TAB(  36); "MSEC" 
990    PRINT    TAB(  10); "WIDTH";  TAB(  27); WW;  TAB(  36); "MSEC" 
1000    PRINT  "QUENCH  PULSE  DELAY";  TAB(  27);QD;  TAB(  36); "MSEC" 
1010    PRINT    TAB(  14); "WIDTH";  TAB(  27);QW;  TAB(  36); "MSEC" 
1020    PRINT  "DUTY  CYCLE  TIME";  TAB(  27);DC;  TAB(  36);"MSEC" 
1060    PRINT  "SCAN  MODE";  TAB(  27);MD$ 

1070  PRINT  "NUMBER  OF  PASSES";  TAB(  27);PA5t 

1080  RETURN 

1300  POKE  33,9:  POKE  32,26:  POKE  34,2:  HOME 

1310  PRINT  GW:  PRINT  DD:  PRINT  DS :  PRINT  DW:  PRINT  WD:  PRINT 

WW:  PRINT  QD:  PRINT  QW:  PRINT  DC:  PRINT  MD$:  PRINT  PA$ 

1 320  RETURN 

1500  IF  EX  =  1  GOTO  1540 
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1505  IF  A  =  3  THEN    GOSUB  7500:  RETURN 

1510  IF  A  =  17  THEN    POP  :  GOTO  9000 

1515  IF  A  =  11  THEN    GOSUB  6500:  RETURN 

1520  IF  A  =  27  THEN    POP  :  GOTO  8000 

1530  IF  A  =  16  THEN    POP  :EX  =  1:  ON  MD^  GOTO  200,400,260,260 

1540  IF  A  =  1  AND  EX  <    >  0  THEN    POP  :  POP  :PP?(0,1)  =  GW  * 

TW:TD  =  DD:KV$  =    -  23912:  CALL  2048:EX  =  0:  GOTO  900 

1545  IF  A  =  18  THEN    GOSUB  100:  GOSUB  1300:  RETURN 

1570  IF  A  =  5  THEN    CALL  62450 :ME  =  0:  RETURN 

1580  IF  A  =  7  THEN    POKE    -  16304,0:  RETURN 

1590  IF  A  =  20  THEN    POKE    -  16303,0:  RETURN 

1610  IF  A  =  13  THEN    GOSUB  1700:  RETURN 

1620  IF  A  =  43  THEN  SD  =  1 :  GOSUB  4000:  GOSUB  1300:  RETURN 

1630  IF  A  =  45  THEN  SD  =    -  1:  GOSUB  4000:  GOSUB  1300:  RETURN 

1640  IF  A  r  8  THEN  ?V%  =  PV$  -  1 

1650  IF  A  =  21  THEN  PV$  =  PVJ  +  1 

1660  IF  PV^  <  3  THEN  PV?  =  13 

1670  IF  PV$  >  13  THEN  PV%  =  3 

1680  POKE  32,25:  POKE  33,1:  POKE  34,2:  HOME  :  VTAB  PV?:  PRINT 

">":  VTAB  PVJ:  RETURN 

1700  POKE  32,26:  POKE  33,9:  VTAB  HTAB  1:  GOSUB  800 

1710  IF  ER  =  1  THEN  RETURN 

1720  IF  PAR  <  0  THEN  RETURN 

1730  1?%  =  TW  »  PAR:PAR  =  TP?  /  TW 

1740  ON  PV$  -  2  GOTO  1780  , 1790, 1800, 1810, 20, 1830, 1840, 1850, 

1860,1900,1910 

1770  GOTO  1500 


1780  GW 

PAR 

GOTO 

1920 

1790 

DD 

PAR 

:  GOTO 

1920 

1800 

DS 

PAR 

GOTO 

1920 

1810 

DW 

PAR 

GOTO 

1920 

1820 

WD 

PAR 

GOTO 

1920 

1830 

WW 

PAR 

GOTO 

1920 

1840 

QD 

PAR 

GOTO 

1920 

1850 

QW 

PAR 

GOTO 

1920 

i860 

DC 

PAR 

GOTO 

1920 

1900 

IF 

EX  = 

0  THEN 

MD% 

1905    PRINT  BE$:  GOSUB  1300:  RETURN 
1910  PA?  =  PAR:  GOTO  1920 

1920  SD  =  0:  GOSUB  4000:  GOSUB  I3OO:  RETURN 
2000    TEXT  :  HOME 

2010  PRINT  "MO  EXPERMENT  CONDITIONS" 
2020    PRINT  :  PRINT  "ELECTRON  BEAM  :" 

2030    HTAB  5:  PRINT  "GRID  POTENTIAL";  TAB(  27);GV;  TAB( 
36); "VOLT" 

2040    HTAB  5:  PRINT  "FILAMENT  POTENTIAL";  TAB(  27);FV;  TAB( 
36); "VOLT" 

2050    HTAB  5:  PRINT  "COLLECTOR  POTENTIAL";  TAB(  27);CV;  TAB( 
36); "VOLT" 

2070    HTAB  5:  PRINT  "COLLECTOR  CURRENT";  TAB(  27);CC;  TAB( 
36);"NAHP" 

2080    PRINT  "CELL  POTENTIALS:  UPPER";  TAB(  27);UV;  TAB( 
36); "VOLT" 

2090    HTAB  18:  PRINT  "LOWER";  TAB(  27);LV;  TAB(  36); "VOLT" 
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2100  HTAE  18:  PRINT  "SIDE(F)";  TA6(  27);SF;  TAB(  36);"VOLT" 

2105  HTAB  18:  PRINT  "SIDE(C)";  TAB(  27);SC;  TAB(  36); "VOLT" 

2110  HTAB  18:  PRINT  "END1";  TAB(  27);E1;  TAB(  36); "VOLT" 

2115  HTAB  18:  PRINT  "END2";  TAB(  27);E2;  TAB(  36); "VOLT" 

2120  PRINT  "MAGNETIC  FIELD  STRENGTH";  TAB(  27);MF;  TAB( 
36); "KG" 

2130  PRINT  "TEMPERATURE";  TAB(  27); TO;  TAB(  36);"C" 

2140  PRINT  "R.F.  LEVEL";  TAB(  27);EP;  TAB(  36);"MV" 

2150  PRINT  "W1  FREQUENCY";  TAB(  27);W1;  TAB(  36);"KHZ" 

2160  PRINT  "W2  FREQUENCY";  TAB(  27);WF;  TAB(  36); "KHZ" 

2165  PRINT  "PRESSURE";  TAB(  27);PR;  TAB(  36);"TORR" 

2170  RETURN 

2200  POKE  33,9:  POKE  32,26:  POKE  34,3:  HOME 

2210  PRINT  GV:  PRINT  FV:  PRINT  CV:  PRINT  CC:  PRINT  UV:  PRINT 

LV:  PRINT  SF:  PRINT  SC:  PRINT  El:  PRINT  E2:  PRINT  MF:  PRINT 

TC:  PRINT  EP:  PRINT  W1 :  PRINT  WF:  PRINT  PR 

2230  RETURN 

2405  GOSUB  2000 

2410  GOSUB  2200 

2415  POKE  32,25:  VTAB  TV$:  HTAB  1 

2420  GET  A$:A  =    ASC  (A$) 

2430  IF  A  =  8  THEN  TV^  =  TVJ  -  1 :  GOTO  2470 

2440  IF  A  =  21  THEN  TV$  =  TV$  +  1 :  GOTO  2470 

2450  IF  A  =  27  THEN    TEXT  :  VTAB  20:  GOTO  8000 

2460  IF  A  =  13  THEN  2500 

2470  IF  TV$  <  4  THEN  1^%  =  19:  GOTO  2415 

2480  IF  TV$  >  19  THEN  TV$  =  4:  GOTO  2415 

2490  GOTO  2415 

2500  POKE  32,26:  VTAB  TVJ:  HTAB  1:  GOSUB  800 

2510  IF  ER  =  1  GOTO  2410 

2520  ON  TV5t  -  3  GOTO  2540  ,2550 ,2560 , 2580 , 2590 , 2600 ,261 0 ,261 5 , 

2620,2625,2630,2640,2650,2660,2670,2680 

2530  GOTO  2410 


2540 

GV 

PAR 

:  GOTO 

2410 

2550 

FV 

PAR 

:  GOTO 

2410 

2560 

CV 

PAR 

:  GOTO 

2410 

2580 

CC 

PAR 

:  GOTO 

2410 

2590 

UV 

PAR 

:  GOTO 

2410 

2600 

LV 

PAR 

:  GOTO 

2410 

2610 

SF 

PAR 

:  GOTO 

2410 

2615 

SC 

PAR 

:  GOTO 

2410 

2620 

El 

PAR 

:  GOTO 

2410 

2625 

E2 

PAR 

GOTO 

2410 

2630 

MF 

PAR 

GOTO 

2410 

2640 

TC 

PAR 

GOTO 

2410 

2650 

EP 

PAR 

GOTO 

2410 

2660 

W1 

PAR 

GOTO 

2410 

2670 

WF 

PAR: 

GOTO 

2410 

2680 

PR 

PAR: 

GOTO 

2410 

4000    ON  PVJ  -  2  GOTO  4020 , 4030 , 4040 , 4050 , 4060 ,4070 , 4080 , 4090 , 

4100,4140,4150 
4010  RETURN 

4020  UL  =  25:LL  =  0:S  r  SD  »  .5:V  =  GW:  GOSUB  4l60:GW  =  V: 
GOSUB  30:  RETURN 


■    •  t'*  ■  >     '  '■-         •■  J 

^    •^         '      -   ;r     •   ■  „■  . 
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M030  UL  =  QD  -  DW:LL  =  0  :S  =  SD  »  DS :  V  =  DD:  GOSUB  m60:DD  = 

V:TD  =  DD:  GOSUB  10:  RETURN 
4040  UL  =  100:LL  r  0:S  =  SD  «  5:V  =  DS:  GOSUB  4l60:DS  =  V: 

RETURN 

4050  UL  =  25:LL  =  1.5:S  =  SD  »  .5:V  =  DW:  GOSUB  4l60:DW  =  V: 

GOSUB  10:  RETURN 
4060  UL  =  16000:LL  =  0:S  =  SD  »  10:V  =  WD:  GOSUB  4l60:WD  =  V: 

GOSUB  20:  RETURN 
4070  UL  =  16000:LL  =  0:S  =  SD  «  10:V  =  WW:  GOSUB  4l60:WW  = 

V:SW  =  WW:  GOSUB  20:  RETURN 
4080  UL  =  DC  -  QW:LL  =  0:S  =  SD  »  100:V  =  QD:  GOSUB  4l60:QD  = 

V:  GOSUB  30:  RETURN 
4090  UL  =  10:LL  =  0:S  =  SD  »  1:V  =  QW:  GOSUB  4l60:QW  =  V: 

GOSUB  30:  RETURN 
4100  UL  =  16000:LL  =  QD  +  QW:S  =  SD  «  100:V  =  DC:  GOSUB 
4160:DC  =  V:  GOSUB  30:  RETURN 

4140    IF  EX  =  0  THEN  UL  =  4:LL  =  1:S  =  SD  »  1 :V  =  MDJ:  GOSUB 
4160:MD$  =  V:  RETURN 
4145    PRINT  BE$:  RETURN 

4150  UL  =  256:LL  =  1:S  =  SD  •  1:V  =  ?k%:  GOSUB  4160:PA$  =  V: 
RETURN 

4160  V  =  V  +  S:  IF  V  >  UL  THEN  V  =  UL:  PRINT  BE$:  RETURN 
4170    IF  V  <  LL  THEN  V  =  LL:  PRINT  BE$ 
4180  RETURN 

5000    TEXT  :  PRINT  :  FLASH  :  PRINT  "INSERT  DATA  DISK":  NORMAL 
5010    INPUT  "ENTER  FILENAME  :  ";DN$ 
5020    ONERR    GOTO  5800 
5030    PRINT  D$;"VERIFY";DN$ 

5040    POKE  216,0:  PRINT  :  PRINT  "FILE  ALREADY  EXISTS! ! !";BE$: 

PRINT  "RE"; 
5045    GOTO  5010 
5050    ONERR    GOTO  5810 

5060  PRINT  D$; "OPEN";DN$:  PRINT  D$; "WRITE";DN$:  PRINT  TW  »  N: 
FOR  I  =  0  TO  N  -  1  :  ON  MD5t  GOTO  506l  ,5061 ,5062,5062 

5061  PRINT  I:  GOTO  5064 

5062  PRINT  I  »  DS:  GOTO  5064 

5064    PRINT  DA(I,0):  NEXT  :  PRINT  D$; "CLOSE" ;DN$:  PRINT 

D$;"VERIFY";DN$ 
5070    POKE  216,0 

5100    PRINT  :  PRINT  "PARAMETER  HARDCOPY  (Y/N)?  "; 
5110    INPUT  A$:  IF  A$  =  "N"  GOTO  8000 
5120    IF  A$  <    >  "Y"  GOTO  5100 

5122    GOSUB  5125:  PRINT    CHR$  ( 1 2) ;D$; "PR#0" :  GOTO  8000 
5125    PRINT  :  INPUT  "DESCRIPTION:  ";DS$ 

5130    PRINT  D$;"PR#1":  PRINT    CHR$  (29 ) ;D$; "IN#1 " :  PRINT  CHR$ 
(23);"CS":  PRINT  D$;"PR#0":  INPUT  T$:  PRINT  D$; "IN#0" :DA$  = 
MID$  (T$,3,8):T$  =    MID$  (T$,12,8) 
5135    PRINT  D$;"PR#1":  PRINT    CHR$  (9);"T" 
5140    PRINT  :  PRINT  "FILENAME:  ";DN$:  PRINT  :  PRINT 

"DESCRIPTION:  ";DS$:  PRINT  :  PRINT  "DATE:  ";DA$:  PRINT  : 
PRINT  "TIME:  ";T$:  PRINT  :  GOSUB  930:  PRINT  :  GOSUB  2010: 
RETURN 

5400    GOSUB  5405:  GOTO  8000 

5405    TEXT  :  HOME  :  INPUT  "ENTER  PARAMETER  FILENAME:  ";PN$:  IF 
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LEN  (PN$)  =  0  THEN  PN$  =  "EXP  PARAMETERS" 
5410    ONERR    GOTO  5150 

5420    PRINT  D$ ; "OPEN" ;PN$;",D1";D$; "READ" ;PN$ 

5430    INPUT  GW:  INPUT  DD:  INPUT  DS:  INPUT  DW:  INPUT  WD:  INPUT 

WW:  INPUT  QD:  INPUT  QW:  INPUT  DC:  INPUT  HD%:  INPUT  Ph% 
5440    INPUT  GV:  INPUT  FV:  INPUT  CV:  INPUT  CC:  INPUT  UV:  INPUT 

LV:  INPUT  SF:  INPUT  SC:  INPUT  El:  INPUT  E2:  INPUT  MF:  INPUT 

PR:  INPUT  EP:  INPUT  W1:  INPUT  WF:  INPUT  TC 
5450    PRINT  D$; "CLOSE" ;PN$:  POKE  216,0:  RETURN 
5500    GOSUB  5505:  GOTO  8000 

5505    TEXT  :  HOME  :  INPUT  "ENTER  PARAMETER  FILENAME:  ";PN$:  IF 

LEN  (PN$)  =  0  THEN  PN$  =  "EXP  PARAMETERS" 
5510    ONERR    GOTO  5150 

5520    PRINT  D$;"0PEN";PN$;",D1";D$;"WRITE";PN$ 

5530    PRINT  GW:  PRINT  DD:  PRINT  DS:  PRINT  DW:  PRINT  WD:  PRINT 

WW:  PRINT  QD:  PRINT  QW:  PRINT  DC:  PRINT  MD*:  PRINT  Pk% 
5540    PRINT  GV:  PRINT  FV:  PRINT  CV:  PRINT  CC:  PRINT  UV:  PRINT 

LV:  PRINT  SF:  PRINT  SC:  PRINT  El:  PRINT  E2:  PRINT  MF:  PRINT 

PR:  PRINT  EP:  PRINT  W1 :  PRINT  WF:  PRINT  TC 
5550    PRINT  D$; "CLOSE" ;PN$:  POKE  216,0:  RETURN 
5800    IF    PEEK  (222)  =  6  GOTO  5050 

5810    PRINT  :  PRINT  "DISK  ERROR!  ERROR  NUMBER  =  ";  PEEK 
(222);BE$ 

5820    PRINT  :  PRINT  "PRESS  ANY  KEY  TO  CONTINUE":  GET  A$ 
5830    POKE  216,0:  GOTO  8000 

5900    TEXT  :  PRINT  D$; "CATALOG":  PRINT  :  PRINT  "PRESS  km  KEY 

TO  CONTINUE"; 
5910    GET  A$:  GOTO  8000 
6000    TEXT  :  HOME 

6005    PRINT  "ENTER  A  ZERO  FOR  THE  UNKNOWN  PARAMETER":  PRINT 
6010    INPUT  "ION  MASS  (G/M)  ?  ";MI 
6020    INPUT  "MAGNETIC  FIELD  STRENGTH  (TESLA)  ?  ";M 
6030    INPUT  "CYCLOTRON  FREQUENCY  (KHZ)  ?  ";CF 
6050    IF  MI  <    =0  THEN  MI  =  MF  *  1 .53575E7  /  (CF  *  1000): 
GOTO  6100 

6060    IF  M  <    =0  THEN  M  =  CF  »  1000  *  MI  /  1.5357E7:  GOTO 
6100 

6070  CF  =  M  »  1.5357E7  /  (MI  »  1000) 

6100    PRINT  :  PRINT  "ION  MASS  (G/M)  ";MI 

6105    PRINT  "MAGNETIC  FIELD  (TESLA)  ";M 

6110    PRINT  "CYCLOTRON  FREQ.   (KHZ)  ";CF 

6150    PRINT  :  PRINT  "PRESS  ANY  KEY  TO  CONTINUE":  GET  A$ 

6200    GOTO  8000 

6500    TEXT  :  HOME  :  IF  MD$  =  4  OR  MD%  =  3  GOTO  6507 
6502    PRINT  :  PRINT  "SCAN  MODE  ERROR" ;BE$:  GOTO  6590 

6507  INPUT  "FIRST, LAST  POINTS:  ";FP,LP 

6508  INPUT  "OFFSET:  ";KS 

6510    IF  FP  <  1  OR  LP  >  N  THEN    PRINT  "LIMIT  ERROR" ;BE$:  GOTO 
6590 

6520    IF  FP  >    =  LP  THEN    PRINT  "RANGE  ERROR" ;BE$:  GOTO  6590 
6530    PRINT  :TK  =  0:NK  =  0:  FOR  I  =  FP  TO  LP  -  1:  PRINT  I;"  , 

";I  +  1;  TAB(  15);"I/I0  =  "; 
6533    IF  DA(I  -  1,0)  =0  THEN  R  =  0:  GOTO  6537 
6535  R  =  (DA(I,0)  -  KS)  /  (DA(I  -  1,0)  -  KS):TK  =  TK  +  R 
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6537  PRINT  R:NK  =  NK  +  1 :  NEXT 

6538  AK  =  TK  /  NK:  PRINT  :  PRINT  "AVERAGE  I/IO  =  ";AK 

6540  PK  =    -    LOG  (AK)  »  (TO  +  273)  »  1000  /  (DS  *  3.53539E16 

»  273):  PRINT  :  PRINT  "PRELIMINARY  K  =  ";PK 
6590    PRINT  :  INPUT  "SATISFACTORY  (Y:N)  ?  ";A$:  IF  A$  r  "Y" 

THEN    GOSUB  920:  RETURN 
6595    IF  A$  =  "N"  GOTO  6507 
6597    GOTO  6590 

7000  DN$  =  "NONE":  GOSUB  5125:  PRINT  :  ON  MDJ  GOTO 

7010,7010,7020,7020,7030,7030,7030 
7010    PRINT  "POINT";  TAB(  1 0) ; "INTENSITY" :  FOR  I  =  0  TO  N  -  1 : 

PRINT    TAB(  2);I  +  1;  TAB(  10);DA(I,0):  NEXT  :  GOTO  7040 
7020    PRINT  "POINT";  TAB(  10); "TIME";  TAB(  20) ; "INTENSITY": 

FOR  I  =  0  TO  N  -  1  :  PRINT    TAB(  2);I  + '^  ;  TAB(  10)  ;I  »  DS ; 

TAB(  20);DA(I,0):  NEXT 
7040    PRINT    CHR$  (1 2) ;D$; "PR#0" :  GOTO  8OOO 
7400    POKE  49360,0:  POKE  49362,0:  PRINT  D$;"BLOAD  KIM  MO 

CODE":  POKE  49361 ,0 
7410    FLASH  :  PRINT  :  PRINT  "KIM  INITIALIZATION  INSTRUCTIONS 

:":  NORMAL 

7420  PRINT  :  PRINT  "ENTER  KIM  KEYBOARD  INSTRUCTIONS  :":  PRINT 
:  PRINT  "  RS":  PRINT  "  AD  00F1":  PRINT  "  DA  04":  PRINT 
"    AD  17FA":  PRINT  "      DA  BE":  PRINT  "      +  A2" 

7430    PRINT  :  PRINT  "PRESS  ANY  KEY  TO  CONTINUE":  GET  A$:A  = 
PEEK  (49362)  +    PEEK  (49363):  GOTO  8OOO 

7500    GOSUB  7800:  POKE    -  16304,0 

7510  CP  =    INT  (  PDL  (1)  »  4  +    PDL  (0)  /  8):A  =    PEEK  (  - 
16384):  IF  A  >  127  THEN    POKE    -  16368,0:  POKE    -  16303,0: 
RETURN 

7525  NX  =    INT  (CP  »  (N  -  1  )  /  IO3I  +  .5) 

7528  IF  NL  =  NX  GOTO  7510 

7529  IF  ME  <    >  0  THEN    XDRAW  1  AT  XL,YL 

7530  ON  MD%  GOTO  7540,7540,7560,7570 

7540  XV  =  NX:YV  =  DA(NX,0):DV  =  NX:CI  =  YV:  GOTO  7600 

7560  XV  =    INT  (NX  »  DS  »  XS  +  .5):YV  =  DA(NX,0):CI  =  YV:DV  = 

NX  »  DS:  GOTO  7600 
7570  XV  =    INT  (NX  »  DS  *  XS  +  .5):YV  =  DA(NX,0)  +  (AA  +  B  » 

NX): CI  =  DA(NX,0):DV  =  NX  «  DS 
7600    GOSUB  50:  XDRAW  1  AT  XV,YV:ME  =  1:XL  =  XV:YL  =  YV 
7700    GOSUB  7900:NL  =  NX:  GOTO  7510 

7800    POKE  32,0:  POKE  33,40:  POKE  34,20:  HOME  :  PRINT 

"POINT:":  PRINT  "X  VALUE:":  PRINT  "INTENSITY:  ":  GOSUB  7900: 
RETURN 

7900    POKE  32,11:  POKE  33,20:  POKE  34,20:  HOME  :  PRINT  NX  +  1 : 

PRINT  DV:  PRINT  CI:  RETURN 
8000  F%  =    FRE  (0):EX  =  0:  POKE    -  16254,0::  POKE    -  16297,0 
8010    TEXT  :  HOME  :  HTAB  15:  FLASH  :  PRINT  "MO  PROGRAM": 

NORMAL 
8020    PRINT  :  PRINT 
8030    PRINT  "COMMANDS  :" 
8040    PRINT  "1 ., .EXIT  PROGRAM" 
8050    PRINT  "2... STORE  DATA  ON  DISK" 
8060    PRINT  "3... CATALOG  DISK" 

8070    PRINT  "4... UPDATE  EXPERIMENTAL  CONDITIONS" 
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8080    PRINT  "5... CHANGE  PARAMETERS" 
8090    PRINT  "6... MASS  CALCULATION" 
8100    PRINT  "7... LOAD  PARAMETER  FILE" 
8110    PRINT  "8... SAVE  PARAMETER  FILE" 
8120    PRINT  "9... PRINT  DATA" 
8130    PRINT  "10. .LOAD  KIM  PROGRAM" 
8310    PRINT  :  INPUT  "COMMAND  ?  ";CMD$ 

8320    ON  CMD$  GOTO  9000,5000,5900,2405,900,6000,5400,5500, 
7000,7400 

8330    PRINT  :  PRINT  "ILLEGAL  COMMAND  !";BE$ 
8340    GOTO  8310 
9000    TEXT  :  END 

10000  D%  =  0:PA$  =  0:KV$  =  0:KE$  =  00 
10010    DIM  PP$(8,1),DA(279,1) 

10020  TW  =  2:WN  =  ^ =  3:TV*  =  4:FE  =  .625:EX  =  0:  R0T=  0 
SCALE=  2:ME  =  0:BE$  r  CHR$  (7)  +  CHR$  (7)  +  CKR$  (7):D$ 
CHR$  (13)  +    CHR$  (4) 

10030  PN$  =  "EXP  PARAMETERS":  GOSUB  5410:TD  =  DD:  GOSUB  10: 
GOSUB  20:  GOSUB  30 

10040    HGR  :  HC0L0R=  3:  TEXT  :  GOTO  7400 


APPENDIX  B 
ASSEMBLY  LANGUAGE  PROGRAMS 


The  three  assembly  language  programs  used  to  run  the  icr  experiment 
are  detailed  in  this  section.    All  programs  were  prepared  as  text  files 
with  a  text  editor,  PIE,  then  converted  to  machine  code  by  an  inter- 
active assembler,  LISA  [122].    The  manuals  should  be  consulted  for 
language  syntax  and  operation  of  the  assembler  and  editor.    Source  text 
listings  for  the  three  programs  end  with  the  word  "TEXT"  in  their 
filenames,  while  the  corresponding  assembled  code  are  stored  as  binary 
files  ending  in  "CODE". 

Apple  Program 

The  program  execution  sequence  for  the  Apple  machine  code  is  shown 
in  Figure  B.I.    As  explained  in  the  previous  appendix,  the  program 
expects  the  first 'BASIC  variable  to  be  the  data,  the  second  the  pass 
count,  the  third  the  program  vector,  and  the  fourth  the  KIM  availability 
status.    Again,  the  pulse  parameter  array  must  be  the  first  one 
declared.    The  program  code  begins  at  800  hex  or  2048  in  decimal,  thus, 
the  program  is  started  by  the  BASIC  statement  "CALL  2048"  in  line  number 
100  of  the  program  listed  in  Appendix  A. 

Upon  entering  the  routine,  the  system  registers  are  saved  and  the 
data  not  ready  flag  is  set.    Monitor  routines  in  ROM,  lOSAV  and  lOREST, 
save  and  restore  the  Apple  registers,  respectively.    This  allows  the 
program  to  be  called  at  any  time  without  disturbing  BASIC  or  any  other 
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Figure  B.l.     Flowchart  of  the  assembly  language  program  APPLE  ICR-MO 
showing  the  main  program  execution  sequence  (part  A)  and 
the  data  transfer  subroutine  (part  B)  which  passes  data 
and  parameters  between  the  Apple  and  KIM. 
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routine.    Next  the  KIM's  availability  status  is  examined,  and  if  the  KIM 
is  not  available  the  Apple  keyboard  status  port  is  checked  for  the 
pressed  key  flag  (bit  7).    If  a  key  has  been  pressed,  the  system 
registers  are  restored  and  program  execution  returns  to  BASIC  with  the 
data  not  ready  flag  set. 

When  the  KIM  is  available,  the  interface  is  configured  for  the 
S-100  memory  "on"  and  the  S-100  bus  requested.    The  data  and  KIM  program 
parameters  are  exchanged  and  then  the  interface  is  configured  "off"  by 
releasing  the  S-100  bus  and  switching  the  Apple  memory  "on".    Next  a  NMI 
request  is  sent  to  the  KIM  to  start  program  execution.    Now  the  data  not 
ready  flag  is  cleared,  registers  are  restored,  and  the  program  returns 
to  BASIC  control.    Since  data  and  program  parameters  are  passed  at  the 
same  time,  each  run  through  the  code  supplies  the  Apple  with  data  from 
the  last  experiment  vrtiile  sending  the  KIM  pulse  parameters  for  the 
next. 

Dual  microcomputer  timing  coordination  is  maintained  through  the 
KIM's  status  port.    Two  situations  are  possible  and  communicated  through 
bit  7  of  the  port.    A  one  indicates  the  KIM  is  busy  with  a  pulse 
sequence  and  should  not  be  disturbed.    A  zero  signifies  the  KIM  has 
completed  a  task  and  is  available  for  service  without  interfering  with 
an  experiment . 
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The  APPLE  ICR-MO  TEXT  Listing  follows.    When  the  program  is 
assembled,  it  occupies  memory  locations  800  to  869  hex  made  available 
the  initialization  procedure  described  in  Chapter  4. 


INS 


DEFINE  VARIABLES  AND  ADDRESSES 


im 

EPZ 

$69 

BASIC  INTEGER  VARIABLE  ADDRESS 

ARVEC 

EPZ 

$6B 

APPLE  ARRAY  ADDRESS 

TIMON 

EQU 

$A000 

PULSE  PARAMETER  ADDRESS 

PRGVEC 

EQU 

$A030 

KIM  PROGRAM  VECTOR 

MP  A  *N*N 

Ron 

KIM  NUMBER  OF  PASSES 

DATA 

EQU 

$A033 

KIM  DATA  ADDRESS 

KYBD 

EQU 

$C000 

APPLE  KEYBOARD  STATUS 

SI  BUS 

EQU 

$CODO 

S-100  BUS  REQUEST  ADDRESS 

S1MEM 

EQU 

$C0D2 

S-100  MEMORY  ENABLE  ADDRESS 

APMEM 

EQU 

$C0D6 

; APPLE  MEMORY  ENABLE  ADDRESS 

KSTAT 

EQU 

$C0D1 

KIM  STATUS  PORT  READ  ADDRESS 

KRES 

EQU 

$CQD2 

KIM  RESET  ADDRESS 

KNMI 

EQU 

$C0D3 

KIM  NMI  ADDRESS 

KSPORT 

EQU 

$C51F 

;KIM  STATUS  PORT  WRITE  ADDRESS 

lOSAV 

EQU 

$FF4A 

; APPLE  REGISTER  SAVE  PROGRAM 

lOREST 

EQU 

$FF3F 

; APPLE  REGISTER  RESTORE  PROGRAM 

ORG 

$800 

,ICR  EXP.  ENTRY  POINT 

THIS  ROUTINE  WILL  SET  THE  INITIAL  PARAMETERS 
REQUIRED  TO  OPERATE  THE  MO  PULSE  ROUTINE 
AND  COMMUNICATE  WITH  THE  KIM. 
THE  PROGRAM  IS  TOTALLY  RELOCATABLE. 
ALL  REGISTERS  ARE  RESTORED. 


ICRMO 


KWAIT 


KRDY 


JSR 

lOSAV 

;SAVE  ALL  REGISTERS 

LDY 

#$18 

;SET  DATA  NOT  READY  FLAG 

LDA 

#$01 

STA 

(INT) ,Y 

;PASS  TO  APPLESOFT 

BIT 

KSTAT 

; CHECK  KIM  STATUS 

BPL 

KRDY 

;BRANCH  IF  KIM  READY 

BIT 

KYBD 

; CHECK  KEYBOARD  STATUS 

BPL 

KWAIT 

;KEEP  CHECKING 

CLC 

;EXIT  FOR  ABORT 

BCC 

EXIT 

STA 

SI  BUS 

; REQUEST  S-100  BUS 

STA 

S1MEM 

; ENABLE  S-100  MEMORY 

JSR 

COMM 

; TRANSFER  DATA  WITH  KIM 

STA 

S1BUS+1 

; RELEASE  S-100  BUS 

STA 

APMEM 

; ENABLE  APPLE  MEMORY 
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LDA  KNMI 


I  START  KIM 


LDY  #$18 
LDA  #$00 
STA  (INT),Y 

EHT      JSR  lOREST 
RTS 


; CLEAR  DATA  NOT  READY  FLAG 

;PASS  TO  APPLESOFT 

; RESTORE  ALL  REGISTERS 
; RETURN 


THIS  ROUTINE  TRANSFERS  THE  NECESSARY  PARAMETERS 
BETWEEN  THE  KIM  AND  APPLE. 


COMM      LDA  DATA 

LDY  #$03 

STA  (INT),Y 

LDA  DATA+1 
DEY 

STA  (INT),Y 

LDA  #$80 

STA  KSPORT 

JSR  XFER 

LDY  #$0A 

LDA  (INT),Y 

STA  NPASS 

LDY  #$11 

LDA  (INT),Y 

STA  PRGVEC 
DEY 

LDA  (INT),Y 

STA  PRGVEC+1 
RTS 


;GET  LOW  BYTE 

;PASS  DATA  TO  APPLESOFT 

;GET  HIGH  BYTE 

;SET  TO  STORE  HIGH  BYTE 


;SET  KIM  NOT  AVAILABLE 


TRANSFER  PULSE  PARAMETERS 

SET  FOR  PASS  COUNT 

GET  PASS  COUNT 

SEND  TO  KIM 

SET  FOR  PROGRAM  VECTOR 

GET  LOW  BYTE 

SEND  TO  KIM 

SET  FOR  HIGH  BYTE 

GET  HIGH  BYTE 

SEND  TO  KIM 

RETURN 


THIS  ROUTINE  TRANSFERS  17  DOUBLE  PRECISION  VALUES 
FROM  THE  FIRST  BASIC  ARRAY  TO  THE  S-100  MEMORY. 


XFER 

LDX 

#!33 

; INITIALIZE  COUNTER 

LDY 

#!42 

; INITIALIZE  POINTER 

NEXT 

LDA 

(ARVEC),Y 

;GET  ARRAY  BYTE 

STA 

TIMON 

;SEND  BYTE 

DEY 

; DECREMENT  COUNTERS 

DEX 

BPL 

NEXT 

; REPEAT  TILL  DONE 

RTS 

; RETURN 

FINADD  EQU 

;END  OF  PROGRAM  ADDRESS 

> 

PAG 

END 
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KIM  Programs 

The  operation  of  the  KIM  requires  two  assembly  language  programs, 
the  slave  controller  and  pulse  generator.    The  two  routines  are 
assembled  separately  then  combined,  along  with  default  timing  and 
control  parameters,  to  form  the  load  time  package  of  KIM  MO  CODE.  For 
instruction  timing  considerations  and  KIM  operating  procedures,  one 
should  consult  the  manuals  [123]. 

Slave  Control 

The  program  KIM  ICR-MO  performs  slave  control  tasks  which  run  the 
pulsed  icr  experiment,  collect  and  average  the  data,  and  light  up  the 
KIM'S  display.    The  program  execution  sequence  can  be  traced  in  Figure 
B.2.    The  routine  is  entered  when  a  NMI  request  is  serviced  by  the  KIM's 
CPU.    After  receiving  the  NMI  the  computer  is  initialized:    the  .  , 

processor's  interrupt  request  line  is  disabled,  the  decimal  mode 
cleared,  and  the  KIM's  status  set  not  available.    An  indirect  jump 
command,  directed  by  the  KIM's  program  vector,  places  the  KIM  in  either 
the  continuous  run  or  signal  averaging  mode. 

The  continuous  run  mode  of  operation,  detailed  in  the  upper  half  of 
Figure  B.2,  first  initializes  the  display  parameters,  zeroes  the  data, 
initializes  the  Bianation  and  pulse  port  hardware,  then  transfers  the 
duty  cycle,  pulse  "on",  and  "off"  times  from  the  common  memory  to  page 
zero  to  satisfy  the  8-bit  pulse  routine.    After  all  initialization 
procedures  have  been  completed,  the  experiment  is  performed  with  the 
specified  timing  parameters.    The  data  from  each  pass  is  fetched  from 
the  Biomation  and  placed  in  the  low  byte  of  the  double  precision  data 
address.    The  program  repeats  until  another  tlMI  command  is  received  by 
the  KIM's  CPU.    Periodic  flickering  of  the  KIM's  display  each  pass 
indicates  the  code  is  functioning  properly. 
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Figure  B.2.  Flowchart  of  the  assembly  language  program  KIM  ICR-MO 
showing  the  continuous  run  and  signal  averaging  modes 
of  operation. 
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Signal  averaging  is  accomplished  through  the  sequence  shown  in  the 
lower  half  of  Figure  B.2.    During  initialization  the  number  of  passes  is 
set  and  execution  continues  similar  to  the  continuous  run  mode.  Next, 
the  KIM's  status  is  set  not  available  and  the  pulse  sequence  is  repeated 
the  number  of  times  specified  in  NPASS.    The  data  from  each  pass  is 
added  to  the  previous  sum  and  stored.    When  all  passes  have  been 
completed,  the  KIM's  status  is  set  available  and  the  program  enters  a 
display  loop  while  it  waits  for  service.    The  resulting  data  is  a  total; 
therefore,  it  must  be  divided  by  the  number  of  passes  for  the  average. 
The  program  remains  in  the  display  loop  until  the  next  NMI  request  is 
received. 

The  following  is  a  listing  of  the  assembly  language  program  KIM 
ICR-MO  TEXT.    When  assembled,  the  program  occupies  memory  locations  A28E 
to  A38E  hex  in  the  S-100  memory. 

INS 

;  DEFINE  VARIABLES  AND  ADDRESSES 
5 


DISD 

EPZ 

$F9 

,KIM  DISPLAY  DATA  BYTE 

DISLA 

EPZ 

$FA 

;KIM  DISPLAY  LOW  ADDRESS 

DISHA 

EPZ 

$FB 

HM  DISPLAY  HIGH  ADDRESS 

PAD 

EQU 

$1700 

KIM  PORT  A  DATA 

PADD 

EQU 

$1701 

KIM  PORT  A  DATA  DIRECTION 

PBD 

EQU 

$1702 

KIM  PORT  B  DATA 

PBDD 

EQU 

$1703 

KIM  PORT  B  DATA  DIRECTION 

PULSER 

EQU 

$A100 

8  BIT  PXSE  ROUTINE  ADDRESS 

SPAD 

EQU 

$F004 

S-100  PARALLEL  PORT  A  DATA 

SPBD 

EQU 

$F005 

S-100  PARALLEL  PORT  B  DATA 

KSTAT 

EQU 

$F01F 

KIM  STATUS  ADDRESS 

EXPB 

EQU 

$F026 

EXP.  PARALLEL  I/O  PORT  B 

KIMDP 

EQU 

$1F1F 

KIM  KEYBOARD  DISPLAY  ROUTINE 

! 

ORG 

$A000  ; 

SET  PULSE  PARAMETERS 

OBJ 

•-$9800 

TIMON 

EQU 

» 

PULSE  BIT  ON  TABLE 

DUTCYL 

EQU 

»+$10  \ 

DUTY  CYCLE  DURATION  TIME 

TIMOFF 

EQU 

»+$1 2  ; 

PULSE  BIT  OFF  TABLE 

! 

ORG 

$A030  ; 

SET  PARAMETER  TABLE 
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OBJ  «-$9800 
PRGVEC  ADR  CHUN 
NPASS    BYT  $01 
DATA      HEX  0000 
RDATA    BYT  $00 

PAG 


; PROGRAM  EXECUTION  VECTOR 
; NUMBER  OF  PASSES 
; DOUBLE  PRECISION  DATA 
;RAW  DATA  STORAGE 


ORG  $A28E 
OBJ  *-$9800 


; PROGRAM  BEGINNING 


THIS  ROUTINE  IS  ENTERED  BY  A  NON-MASKABLE 
INTERRUPT  REQUEST,  INITIALIZES  THE  KIM  AND  INTERFACE, 
THEN  JUMPS  TO  THE  DESIGNATED  ROUTINE. 


SEX 
CLD 

LDA  #$80 
STA  KSTAT 
JMP  (PRGVEC) 


; DISABLE  INTERRUPT 
; CLEAR  DECIMAL  MODE 
;KIM  STATUS:  BUSY 


THIS  ROUTINE  SETS  THE  KIM  FOR  CONTINUOUS  RUN 
OF  THE  ICR-MO  PULSE  ROUTINE. 


CRUN 


CONT 


T  nfl 

•  QiTT  riDiTDfiTTnM  Mnnc 
,oCii  UrCiriftlXUN  nUUCi 

STA 

DISHA 

;AND  UPDATE  DISPLAY 

LDA 

RDATA 

; DISPLAY  LAST  POINT 

STA 

DISD 

LDA 

#$00 

;ZERO  THE  DATA 

STA 

DATA 

STA 

DATA+1 

JSR 

INIT 

; INITIALIZE  THE  PARAMETERS 

JSR 

SETPUL 

;SET  PULSE  PARAMETERS 

LDA 

#$80 

;KIM  STATUS:  BUSY 

STA 

KSTAT 

JSR 

KIMDP 

; DELAY  AND  DISPLAY 

JSR 

BIOARM 

;ARM  THE  BIOMATION 

JSR 

PULSER 

;RUN  EXPERIMENT 

JSR 

BIODAT 

;GET  DATA  FROM  BIOMATION 

STA 

DATA 

;SAVE  DATA 

STA 

DISD 

;DISPLAY  DATA 

LDA 

#$00 

;KIM  STATUS:  AVAILABLE 

STA 

KSTAT 

INC 

DISLA 

;INCREMENT  COUNTER 

JSR 

KIMDP 

;DEALY  AND  DISPLAY 

CLC 

; CONTINUE  RUNNING 

BCC 

CONT 

PAG 

THIS  ROUTINE  SETS  THE  KIM  UP  FOR  DATA  COLLECTION 
SIGNAL  AVERAGING.   WHEN  THE  SPECIFIED  NUMBER  OF  PASSES 
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HAVE  BEEN  TAKEN  THE  KIM  EXECUTES  A  DISPLAY  LOOP, 


DATAVG  LDA  NPASS 
STA  DISLA 
LDA  #$02 
STA  DISHA 
LDA  #$00 
STA  DATA 
STA  DATA+1 
JSR  INIT 
JSR  SETPUL 
LDA  #$80 
STA  KSTAT 


;SET  THE  NO.  OF  PASSES 
;SAVE  WORKING  PASS  COUNT 
;SET  OPERATION  MODE 
;DISPLAY  THE  MODE 
;ZERO  THE  DATA 


; INITIALIZE  THE  PARAMETERS 
;SET  PULSE  PARAMETERS 
;KIM  STATUS:  BUSY 


NEXTP 


JSR  KIMDP 
JSR  BIOARM 
JSR  PULSER 
JSR  BIODAT 
STA  DISD 


; DISPLAY 

;ARM  THE  BIOMATION 

; EXECUTE  A  PASS 

;GET  DATA  FROM  BIOMATION 

; DISPLAY  RAW  DATA 


ac 

ADC  DATA 
STA  DATA 
LDA  DATA+1 
ADC  #$00 
STA  DATA+1 


; CLEAR  FOR  SIGNAL  AVERAGING 
; AVERAGE  SIGNAL 
; STORE  IT  ' 
;HIGH  BYTE 


DEC  DISLA 
BNE  NEXTP 
LDA  #$00 
STA  KSTAT 
LDA  #$00 
STA  DISHA 
WAITLP  JSR  KIMDP 
JMP  WAITLP 


; DECREMENT  PASS  COUNT 
; CONTINUE  AVERAGING 
;KIM  STATUS:  AVAILABLE 

;SET  OPERATION  MODE 
;DISPLAY  IT 
; DISPLAY  VITAL  SIGNS 
; CONTINUE  DISPLAYING 


THIS  SUBROUTINE  WILL  SET  THE  INITIAL  PARAMETERS 
REQUIRED  TO  OPERATE  THE  ICR-MO  PULSE  ROUTINE, 
BIOMATION,  AND  PULSE  OUTPUT  PORT. 
THE  PROGRAM  IS  TOTALLY  RELOCATABLE. 
ALL  REGISTERS  ARE  RESTORED. 


INIT 


PHA 

;SAVE  THE  ACCUMULATOR 

LDA 

#$FF 

;SET  PARALLEL  PORT  A  TO  ALL 

STA 

PADD 

LDA 

#$07 

;KIM  PB0,1,2  ARE  OUTPUTS 

STA 

PBDD 

LDA 

#$87 

; INITIALIZE  THE  BIOMATION 

STA 

SPBD 

LDA 

#$00 

; ACCESS  PIA  DATA  DIRECTION 

STA 

EXPB+1 

;EXP.  PARALELL  PORT  B 

LDA 

#$FF 

;SET  FOR  8  OUTPUT  LINES 

STA 

EXPB 
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LDA  #$04  ; ACCESS  DATA  PORT 
STA  EXPB+1 

FLA  ; RESTORE  ACCUMULATOR 

RTS  ; RETURN 

PAG 

SUBROUTINE  SETPUL  TRANSFERS  THE  PULSE  PARAMETERS  IN 
THE  COMMON  MEMORY  TO  PAGE  ZERO  FOR  USE  BY  THE  KIM. 


SETPUL  LDX  #$21 
NEXTV    LDA  $AOOO,X 

STA  $0000, X 

DEX 

BPL  NEXTV 

RTS 

PAG 


INITIALIZE  COUNTER 
GET  PULSE  PARAMETER 
TRANSFER  TO  ZERO  PAGE 
DECREMENT  COUNTER 
REPEAT  TILL  DONE 
RETURN 


SUBROUTINE  BREQ  SENDS  AN  OUTPUT  COMMAND  TO  THE 
BIOMATION  THEN  WAITS  FOR  THE  OPT  FLAG. 
ALL  REGISTERS  ARE  RESTORED. 

BREQ      PHA  ;SAVE  ACC 

LDA  #$85  JSEND  OPT  COMMAND 

STA  SPAD 

WFLG      LDA  #$01  ; CHECK  OPT  FLAG 

AND  SPBD  ;READ  STATUS  PORT 

BEQ  WFLG  ;WAIT  FOR  OPT  FLG 

PLA  ; RESTORE  ACC 

RTS  ;  RETURN  ■  '"^  ' 

SUBROUTINE  GETB  READS  A  DATA  POINT  FROM  THE  BIOMATION 
THEN  TOGGLES  WDC  FOR  THE  NEXT  POINT. 
THE  DATA  POINT  IS  IN  THE  ACCUMUUTOR  UPON  '     '  """ 
EHT  FROM  THE  ROUTINE. 

GETB      LDA  SPAD  ;GET  DATA  '   '      . ' 

PHA  ;SAVE  IT  .  '  ' 

LDA  #$84  ; TOGGLE  WDC  LOW 

STA  SPAD  ;FOR  NEXT  POINT 
LDA  #$85 
STA  SPAD 

PLA  ; RESTORE  ACC 

RTS  ; RETURN 

THIS  SUBROUTINE  INTERROGATES  THE  BIOMATION  FOR  DATA 
AND  STORES  THE  RESULT  IN  RDATA. 

THE  ACCUMULATOR  CONTAINS  THE  DATA  UPON  EXIT  FROM 
THE  ROUTINE. 

BIODAT  BIT  SPBD  ;WAIT  FOR  RECORD  DONE  FUG 

BMI  BIODAT 


JSR  BREQ  ; REQUEST  BIOMATION  OUTPUT 
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JSR  GETB 
JSR  GETB 
STA  RDATA 
PHA 

LDA  #$87 
STA  SPAD 
PLA 
RTS 


;SKIP  FIRST  POINT 
;GET  DATA 

; STORE  AS  RAW  DATA 

;SAVE  THE  DATA 

; RELEASE  BIOMATION 

; RESTORE  THE  DATA 
; RETURN 


THIS  ROUTINE  ARMS  THE  BIOMATION. 
ALL  REGISTERS  ARE  RESTORED. 


BIOARM  PHA 

LDA  #$07 
STA  SPAD 
LDA  #$87 
STA  SPAD 
PLA 
RTS 


;SAVE  THE  ACCUMULATOR 
;ARM  BIOMATION 


; RESTORE  THE  ACCUMULATOR 
; RETURN 


FINADD  EQU  » 

PAG 
END 


;END  OF  PROGRAM  ADDRESS 
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Pulse  generator 

Pulse  timing  is  controlled  by  executing  a  series  of  computer 
instructions  of  known  length  such  that  the  total  execution  time  is 
constant.    Since  each  instruction  takes  from  two  to  four  computer  clock 
cycles  to  execute  and  the  computer  clock  runs  at  a  constant  frequency, 
the  execution  time  of  each  instruction  is  known.    Even  though 
conditional  branch  statements  vary  the  code  executed  each  time  the 
program  runs,  the  number  of  clock  cycles  is  constant  for  every  possible 
execution  path. 

The  program  execution  sequence  is  outlined  in  Figure  B.3.  Upon 
entering  the  routine,  the  program  expects  the  duty  cycle  time,  "on",  and 
"off"  times  to  be  in  page  zero.    The  program  first  zeroes  the  pulse 
master  clock,  next  the  pulse  port  value  is  set  to  zero.  Now  each  bit's 
"on"  and  "off"  time  is  compared  to  the  pulse  master  clock  by  double 
precision  subtraction.    If  the  clock  is  greater  than  or  equal  to  the 
"on"/"off"  time,  the  carry  is  set  and  the  appropriate  bit  set  or  cleared 
in  the  pulse  port  value.    When  all  eight  bits  are  configured,  the  pulse 
port  value  is  sent  to  the  pulse  port  hardware.    It  is  at  this  time  that 
all  bits  seen  by  the  detection  electronics  change  synchronously.  The 
pulse  master  clock  is  then  compared  to  the  duty  cycle  time.    If  they  are 
equal,  program  execution  is  returned  back  to  the  calling  program;  if 
not,  the  pulse  master  clock  is  incremented  and  the  program  is  delayed 
before  repeating  the  code. 

This  code  can  be  used  on  either  the  Apple  or  KIM  microcomputers 
with  the  appropriate  modification.    The  CPU  phase  zero  clock  frequency 
is  1.00000  MHz  for  the  KIM;  consequently,  each  instruction  cycle  lasts 
one  microsecond.    Thus,  500  instruction  cycles  are  needed  for  a  0.5  msec 
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Figure  B.3.     Flowchart  of  the  8-bit  pulse  generator  subroutine. 
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program  cycle  time .    This  dictates  the  40  CPU  cycle  delay  after 
incrementing  the  pulse  master  clock.    However,  the  Apple  phase  zero 
clock  frequency  is  1.02050  MHz;  thus,  510  CPU  clock  cycles  approximate 
the  510.25  clock  cycles  needed  for  a  0.5  msec  program  cycle  time.  The 
50  CPU  clock  cycle  delay  at  the  end  of  the  routine  introduces  a  timing 
error  of  only  0.05$. 

The  code  used  to  generate  the  pulse  sequences,  KIM  PULSE  8  TEXT, 
follows.    When  assembled,  it  occupies  memory  locations  A100  to  A28D  hex 
in  the  S-100  memory. 


INS 


DEFINE  VARIABLES  AND  ADDRESSES. 

ALL  THE  FOLLOWING  VARIABLES  ARE  16  BITS  LONG. 
HIGH  ORDER  BYTES  ARE  FIRST. 


TIMON    EPZ  $00 

DUTCYL  EPZ  $10 

TIMOFF  EPZ  $12 

PaK      EPZ  $22 


PULSE  BIT  ON  TABLE 
DUTY  CYCLE  TIME 
PULSE  BIT  OFF  TABLE 
PULSE  MASTER  CLOCK 


THE  FOLLOWING  LOCATIONS  ARE  8  BITS  IN  LENGTH. 


BASE2    EQU  $F026 
PLSVAL  EQU  $2FF 

ORG  $A100 
OBJ  »-$9900 

8  BIT  PULSE  ROUTINE 


PULSER  LDA  #$00 

STA  PCLK+1 
STA  PCLK 

f 

PULSE    LDA  #$00 

STA  PLSVAL 


BO 


SEC 

LDA  PCLK+1 
SBC  TIMON+1 
LDA  PCLK 
SBC  TIMON 
BCC  BODELA 


;PIA  B  ADDRESS 
; PULSE  PORT  VALUE 

;8  BIT  PULSE  ROUTINE  ENTRY 
; STORE  AT  $800 


;ZER0  PULSE  MASTER  CLOCK 
;  LOW  BYTE 
;  HIGH  BYTE 

;ZERO  PULSE  PORT  VALUE 


; CHECK  FOR  BIT  0  ON 
;  LOW  BYTE 

;  HIGH  BYTE 

; BRANCH  IF  NOT  ON 
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LDX 

#$01 

JMP 

BOOFF 

BODELA 

LDX 

#$00 

NOP 

BOOFF 

SEC 

LDA 

PCLK+1 

SBC 

TIMOFF+1 

LDA 

PCLK 

SBC 

TIMOFF 

BCC 

BODELB 

LDX 

#$00 

JMP 

BOSET 

BODELB 

NOP 

BOSET 

TXA 

ORA 

PLSVAL 

• 

STA 

PLSVAL 

> 

B1 

SEC 

LDA 

PCLK+1 

SBC 

TIMON+3 

LDA 

PCLK 

SBC 

TIMON+2 

BCC 

B1DELA 

LDX 

#$02 

JMP 

B10FF 

B1DELA 

LDX 

#$00 

NOP 

B10FF 

SEC 

LDA 

PCLK+1 

SBC 

TIMOFF+3 

LDA 

PCLK 

SBC 

TIMOFF+2 

BCC 

B1DELB 

LDX 

#$00 

JMP 

B1SET 

B1DELB 

NOP 

Mno 
NUr 

B1SET 

TXA 

ORA 

PLSVAL 

STA 

PLSVAL 

5 

B2 

SEC 

LDA 

PCLK+1 

SBC 

TIMON+5 

LDA 

PCLK 

SBC 

TIMON+4 

BCC 

B2DELA 

LDX 

#$04 

JMP 

B20FF 

B2DELA 

LDX 

#$00 

NOP 

B20FF 

SEC 

LDA 

PCLK+1 

;SET  BIT  0  ON 

;SET  BIT  0  NOT  ON 
;  DELAY 

; CHECK  FOR  BIT  0  OFF 
;  LOW  BYTE 

;  HIGH  BYTE 

; BRANCH  IF  BIT  0  NOT  OFF 
; CLEAR  BIT  0 

;  DELAY 

;SET  PULSE  PORT  VALUE 
; STORE  IT 

; CHECK  FOR  BIT  1  ON 
;  LOW  BYTE 

;  HIGH  BYTE 

; BRANCH  IF  NOT  ON 
;SET  BIT  1  ON 

;SET  BIT  1  NOT  ON 
; DELAY 

1  CHECK  FOR  BIT  1  OFF 
:  LOW  BYTE 

:  HIGH  BYTE 

;  BRANCH  IF  BIT  1  NOT  OFF 
: CLEAR  BIT  1 

DELAY 

SET  PULSE  PORT  VALUE 
STORE  IT 

CHECK  FOR  BIT  2  ON 
LOW  BYTE 

HIGH  BYTE 

BRANCH  IF  NOT  ON 
SET  BIT  2  ON 

SET  BIT  2  NOT  ON 
DELAY 

CHECK  FOR  BIT  2  OFF 
LOW  BYTE 
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SBC 

TIMOFF+5 

LDA 

PCLK 

SBC 

TIMOFF+4 

BCC 

B2DELB 

LDX 

#$00 

JMP 

B2SET 

B2DELB 

NOP 

MOP 

B2SET 

TXA 

ORA 

PLSVAL 

• 

STA 

PLSVAL 

f 

B3 

SEC 

LDA 

PCLK+1 

SBC 

TIMON+7 

LDA 

PCLK 

SBC 

TIMON+6 

BCC 

B3DELA 

LDX 

#$08 

JMP 

B30FF 

B3DELA 

LDX 

#$00 

NOP 

B30FF 

SEC 

LDA 

PCLK+1 

SBC 

TIMOFF+7 

LDA 

PCLK 

SBC 

TIMOFF+6 

BCC 

B3DELB 

LDX 

#$00 

JMP 

B3SET 

B3DELB 

NOP 

NOP 

B3SET 

TXA 

ORA 

PLSVAL 

STA 

PLSVAL 

B4 


SEC 

; CHECK  FOR  BIT 

4  ON 

LDA 

PCLK+1 

;  LOW  BYTE 

SBC 

TIM0N+I9 

LDA 

PCLK 

;  HIGH  BYTE 

SBC 

TIM0N+I8 

BCC 

B4DELA 

; BRANCH  IF  NOT 

ON 

LDX 

#$10 

;SET  BIT  4  ON 

JMP 

B40FF 

LDX 

#$00 

;SET  BIT  4  NOT 

ON 

NOP 

; DELAY 

SEC 

; CHECK  FOR  BIT 

4  OFF 

LDA 

PaK+1 

;  LOW  BYTE 

SBC 

TIMOFF+19 

LDA 

PCLK 

;  HIGH  BYTE 

SBC 

TIMOFF+IB 

BCC 

B4DELB 

; BRANCH  IF  BIT 

4  NOT 

LDX 

#$00 

; CLEAR  BIT  4 

JMP 

B4SET 

HIGH  BYTE 

BRANCH  IF  BIT  2  NOT  OFF 
CLEAR  BIT  2 

DELAY 

SET  PULSE  PORT  VALUE 
STORE  IT 

CHECK  FOR  BIT  3  ON 
LOW  BYTE 

HIGH  BYTE 

BRANCH  IF  NOT  ON 
SET  BIT  3  ON 

SET  BIT  3  NOT  ON 
DELAY 

CHECK  FOR  BIT  3  OFF 
LOW  BYTE 

HIGH  BYTE 

BRANCH  IF  BIT  3  NOT  OFF 
CLEAR  BIT  3 

DELAY 

SET  PULSE  PORT  VALUE 
STORE  IT 
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B4DELB  NOP 
NOP 

B4SET  TXA 

ORA  PLSVAL 
STA  PLSVAL 

B5  SEC 

LDA  PCLK+1 
SBC  TIM0N+!11 
LDA  PCLK 
SBC  TIMON+I10 
BCC  B5DELA 
LDX  #$20 
JMP  B50FF 

B5DELA  LDX  #$00 
NOP 

B50FF  SEC 

LDA  PCLK+1 
SBC  TIM0FF+I11 
LDA  PCLK 
SBC  TIM0FF+!10 
BCC  B5DELB 
LDX  #$00 
JMP  B5SET 

B5DELB  NOP 
NOP 

B5SET  TXA 

ORA  PLSVAL 
STA  PLSVAL 

B6  SEC 

LDA  PCLK+1 
SBC  TIM0N+I13 
LDA  PCLK 
SBC  TIMON+112 
BCC  B6DELA 
LDX  #$40 
JMP  B60FF 

B6DELA  LDX  #$00 
NOP 

B60FF  SEC 

LDA  PCLK+1 
SBC  TIM0FF+!13 
LDA  PCLK 
SBC  TIM0FF+I12 
BCC  B6DELB 
LDX  #$00 
JMP  B6SET 

B6DELB  NOP 
NOP 

B6SET  TXA 

ORA  PLSVAL 
STA  PLSVAL 


;  DELAY 

;SET  PULSE  PORT  VALUE 
:  STORE  IT 

;  CHECK  FOR  BIT  5  ON 
LOW  BYTE 

:  HIGH  BYTE 

; BRANCH  IF  NOT  ON 
;SET  BIT  5  ON 

;SET  BIT  5  NOT 
: DELAY 

:  CHECK  FOR  BIT  5  OFF 
LOW  BYTE 

HIGH  BYTE 

;  BRANCH  IF  BIT  5  NOT  OFF 
: CLEAR  BIT  5 

DELAY 

;SET  PULSE  PORT  VALUE 
:  STORE  IT 

CHECK  FOR  BIT  6  ON 
LOW  BYTE 

HIGH  BYTE 

BRANCH  IF  NOT  ON 
SET  BIT  6  ON 

SET  BIT  6  NOT  ON 
DELAY 

CHECK  FOR  BIT  6  OFF 
LOW  BYTE 

HIGH  BYTE 

BRANCH  IF  BIT  6  NOT  OFF 
CLEAR  BIT  6 

DELAY 

SET  PULSE  PORT  VALUE 
STORE  IT 


R7 

SEC 

LDA 

PCLK+1 

SRC 

TIMON+ ! 1 5 

LDA 

PCLK 

SBC 

TIMON+ ' 1 4 

RCC 

R7DELA 

U  1  X/XjXjA 

LDX 

#$80 

T  ny 

MOP 

I  DA 

PCI  K+l 

SRC 

TTMOFF+ ! 1 S 

1  DA 

PPT  K 

TTMnFF+ 1 1 4 

RfC 

R7DELR 

.IMP 

R7SET 

B7DELB 

NOP 

NOP 

D  t  OCi  1 

TX  A 

ORA 

PLSVAL 

X  XjO  V  AXj 

• 

STA 

HASE2 

> 

SEC 

LDA 

PCLK+1 

X  OXj  iv~  1 

SBC 

DUTCYL+1 

LDA 

PCLK 

SBC 

DUTCYL 

X^  U  X  w  X  u 

RfC 

• 

RTS 

t 

I  DA 

Xj  x/a 

PCI  K+1 

ADC 

#$01 

STA 

PCLK+1 

I  DA 

PPT  Ff 

ADr 

• 

STA 

PPf  V 

* 

LDX 

#!7 

DEL 

DEX 

BNE 

DEL 

NOP 

• 

MOP 

t 
• 

TMD 

tinr 

DITT 

rULoCi 

FINADD 

EQU 

» 

PAG 

END 
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CHECK  FOR  BIT  7  ON 
LOW  BYTE 

HIGH  BYTE 

BRANCH  IF  NOT  ON 
SET  BIT  7  ON 

SET  BIT  7  NOT  ON 
DELAY 

CHECK  FOR  BIT  7  OFF 
LOW  BYTE 

HIGH  BYTE 

BRANCH  IF  BIT  7  NOT  OFF 
CLEAR  BIT  7 

DELAY 

SET  PULSE  PORT  VALUE 

SEND  PULSE  VALUE  OUT  HARDWARE 

CHECK  FOR  DUTY  CYCLE  END 
LOW  BYTE 

HIGH  BYTE 

BRANCH  IF  NO  DONE 
EHT,  END  OF  DUTY  CYCLE 

INCREMENT  PULSE  MASTER  CLOCK 
LOW  BYTE 

SAVE  IT 
HIGH  BYTE 

SAVE  IT 

START  40  CYCLE  DELAY 
DELAY  LOOP 


;  NEXT  PASS 

;END  OF  PROGRAM  ADDRESS 
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